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Chapitre I Physique et probabilités 


1 Probabilités 


a Exemple : jeux aléatoires 
a Le jeu de pile ou face 

Si on lance une pièce de monnaie N fois et qu’elle tombe N p n e fois sur pile, la 
proportion P p n e = N p n e /N mesure la probabilité de cet évènement. Si on fait tendre N 
tend vers l’infini, on obtient une mesure expérimentale de plus en plus précise. 

Si on liste tous les évènements possibles, dont le nombre 0 vaut ici fl = 2, on peut 
considérer qu’ils doivent être équiprobables (rien ne favorise une face plutôt qu’une 
autre). Comme ce sont des évènements élémentaires, ils vérifient (cf. annexe A) 

Eface S - -Epile = 1 e t Ef ace = Epj] e 
d OÙ Eface = EpH e = 1/2 . 

P Le jeu de dé 

Soit l’expérience qui consiste à lancer N fois un dé à six faces. Chacune des valeurs 
1, 2, • • • ,6 est une issue. Puisqu’on répète N fois le lancer de dé, disons qu’on obtient 
N\ fois le 1, ..., Nq fois le 6, on a: 


Ni + ■ ■ ■ + Nq = N . 


Soit l’événement : “le résultat est i”. On lui associe sa probabilité de réalisation Pi définie 
par le rapport du nombre d’issues favorables au nombre total de cas, soit : 



Comme 0 < iVj < N, on vérifie que la probabilité d’un évènement est un nombre compris 
entre 0 et 1 et on retrouve également la propriété de normalisation P x = 1, qui 

traduit le fait qu’obtenir une des valeurs 1 à 6 est une certitude. 

Considérons maintenant l’événement : “le résultat est pair”. Les issues favorables 
étant 2, 4 ou 6, on en déduit que 


P(“Obtenir un 2 ou un 4 ou un 6”) = 


n 2 + n a + n 6 

N 


N 2 n 4 n 6 

7r + 7r + 7r = ft+J, ‘ + fl > 


Dans cet exemple, les issues favorables constituent des événements incompatibles, 
c’est-à-dire qui ne peuvent se réaliser simultanément (le tirage ne peut pas donner un 2 
et un 4, par exemple). On va considérer maintenant des situations plus compliquées. 
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7 Jeu de cartes 

Considérons par exemple, la probabilité associée à l’événement : “tirer soit un 7, soit 
un carreau dans un jeu de 52 cartes”. Les issues favorables ici sont les treize carreaux 
(4x13 = 52) et trois 7 (car le 7 de carreau a déjà été compté comme issue favorable 
parmi les carreaux et ne doit pas être compté en double). Autrement dit, tirer un 7 
et tirer un carreau ne sont pas deux événements incompatibles puisque le 7 de carreau 
existe dans le jeu. Il y a 13+3=16 issues favorables, la probabilité associée vaut donc: 

16 4 

P (“Tirer soit un 7, soit un carreau dans un jeu de 52 cartes”) = — = — . 

v ’ J ' 52 13 

On retrouve ce résultat en appliquant (A 1) à cet exemple : 

. 4 13 1 4 

.P (Tirer soit un 7, soit un carreau dans un jeu de 52 cartes”) =-1-= — . 

v J ; 52 52 52 13 


b Dénombrement 


Un arrangement à p éléments d’un ensemble à n éléments est une sélection de p 
éléments parmi les n, dans un ordre déterminé (on parle également de permutation de p 
éléments parmi n). Comme il y a n façons de choisir le premier élément de la sélection, 
puis n — 1 de choisir le 2ènre, ... etc., le nombre d’arrangements An vérifie : 

77 t 

^ = n(n-l)-(n-(p-l)) = ^-^ î . 

On notera en particulier que le nombre de permutations de n objets (parmi n) vaut n! 
(par convention 0! = 1). 

Une combinaison à p éléments d’un ensemble à n éléments est une sélection de p 
éléments parmi les n, sans ordre déterminé. Comme le nombre de permutations de p 
objets est pl, le nombre de combinaisons, noté (”) ou Cn , vérifie: 

/n\ _ _ ni 

\p J pl pl{n — p)l 

Par exemple, il y a 3 = Q) combinaisons des lettres a, b et c prises 2 à 2 (ab, ac, 
bc), mais 6= arrangements (ab, ba, ac, ca, bc, cb). 


2 Systèmes physiques à l’équilibre 


a Équilibre thermodynamique 

Par système à l’équilibre, on entend un système dont les propriétés physiques n’évo¬ 
luent plus dans le temps à l’échelle macroscopique. A cette échelle, les états d’équilibre 
sont caractérisés par un petit nombre de variables appelées variables d’état. Ces variables 
ne sont pas toutes indépendantes, mais reliées entre elles par des relations fonctionnelles. 
La célèbre équation d’état du gaz parfait : 


N 

p=v RT 
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2. Systèmes physiques à Véquilibre 


est un exemple d’une telle relation qui associe les variables d’état, pression, densité 
et température. Il est remarquable que la thermodynamique physique macroscopique 
ait été capable d’établir phénoménologiquement de telles relations, à une époque où la 
structure atomique de la matière était inconnue. Nous sommes désormais convaincus que 
toutes les propriétés physiques observées à notre échelle dépendent en dernière analyse 
de la dynamique microscopique des constituants élémentaires de la matière. L’objet de 
la physique statistique est d’établir les relations existant entre les variables 
macroscopiques, à partir des équations qui régissent le comportement de 
la matière à l’échelle microscopique. Les états d’équilibre sont des états asymp¬ 
totiques. On entend par là qu’un système initialement perturbé, qu’on laisse ensuite 
évoluer librement, finit toujours par revenir à un de ses états d’équilibre sur des échelles 
de temps plus ou moins longues. L’étude de l’origine microscopique et la dynamique 
des mécanismes de retour à l’état d’équilibre est l’objet de la physique statistique des 
systèmes hors de l’équilibre. 

b Autres grandeurs thermodynamiques 

Les variables p, V, T et N sont des grandeurs thermodynamiques majeures, mais 
d’autres variables macroscopiques interviennent, qui sont également essentielles. 

Faisons maintenant l’expérience suivante : on mesure la capacité calorifique du gaz en 
fonction de la température en apportant une quantité de chaleur ôQ à un système isolé 
thermiquement, à volume constant, et en regardant de combien a varié sa température 
(on note la variation ST) autour d’une température de consigne Tq : on en tire la capacité 
calorifique à Tq, C v (Tq ) = SQ/ST. La courbe C v (Tq) montre que celle-ci dépend de la 
température. La capacité calorifique est une grandeur thermodynamique très importante 
qui intervient dans les bilans thermodynamiques. Or, la thermodynamique seule ne 
permet pas d’expliquer pourquoi C v (Tq) a ce comportement particulier. C’est le rôle de 
la physique statistique de le comprendre et de le modéliser. 

c Ordres de grandeur 

Considérons une certaine masse m de molécules de diazote (constituant essentiel de 
l’air) dans un ballon de volume V, dont on fixe la température par le bain ambiant. On 
pourra en mesurer la pression et vérifier la loi des gaz parfaits. Examinons les ordres de 
grandeurs des quantités en jeu. 

* La masse d’une molécule de diazote peut être évaluée à partir de la masse molaire 
M = 28 g.mol -1 et on trouve m ~ 5 ICO 26 kg. La taille d’une molécule est de 
l’ordre de a 10 10 m. 

’ Combien il y a-t-il de molécules ? Le nombre qui fait le lien entre le microscopique 
et le macroscopique est le nombre d’Avogadro : 

A/a = 6,02210 23 ; 

c’est le nombre d’atomes qu’il y a dans 12 g de carbone. C’est un nombre ex¬ 
traordinairement grand. Dans les conditions standard, une mole de gaz parfait 
occupe un volume de 22,4 L ce qui correspond à une densité 
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De la densité, on peut évaluer la distance typique entre molécules 

d = ro” 1//3 ~ 33 a . 

* Quelle est la vitesse typique des molécules dans le gaz ? Pour cela, on peut choisir 
une approche expérimentale et percer un petit trou dans le ballon pour mesurer 
la vitesse des particules qui sortent. La formule (19) appliquée à un système 
isotrope donne une distribution des vitesses de la forme : 

V(v) = Av 2 e~ Bv2 . 

En faisant varier les paramètres, on observe que A ~ (w) 3/2 et B = mr ■ On 
trouve donc, quand on normalise, la distribution : 

, f M \ 3/2 A 2 mA 

V{v)= \2^f) 

Le lien entre cette loi microscopique et celle des gaz parfaits, macroscopique, n’a 
rien d’évident, mais il est remarquable que toutes deux sont paramétrées par le 
même rapport RT /M. La distribution des vitesses est telle que l’écart type a est 
de l’ordre de la vitesse moyenne (v ), qui est définie par (v) = J^°vV(v)dv. On 
en déduit l’ordre de grandeur typique de la vitesse dans le gaz 

(v) ~ 500 ms -1 . 

On aurait pu calculer 

VW) = 

qui est toujours du même ordre de grandeur, 1 et donne l’ordre de grandeur de 
l’énergie cinétique des molécules : 

e = -mv 2 — 6,24 1CT 21 J. 

2 

Cet ordre de grandeur est microscopique. Cependant, ramené à A/a molécules, 
le gaz contient une énergie cinétique totale de l’ordre de 3500 J, qui est du 
bon ordre de grandeur des énergies à l’échelle macroscopique. On pourra retenir 
que les énergies à l’échelle atomique sont souvent de l’ordre de grandeur d’une 
fraction d’eV (ici 6,24 1CD 21 J ~ 0,039 eV). 

d Collisions et libre parcours moyen 

En réalité, la présence d’interactions à courte portée entraîne des collisions élastiques 
entre molécules. Pour simplifier, on peut imaginer que les molécules se comportent 
comme des boules dont le diamètre est environ a. Evaluons le nombre moyen de collisions 
que subit une molécule. La section efficace de collision de la molécule (surface apparente 

1. Pour la distribution choisie, quand on normalise correctement P(v), on trouve que la seconde est 
égale à ^3^/8 « 1,085 fois la première. 
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3. Bilan et objectifs du cours 


de la cible) est a ~ ira 2 . Entre deux collisions, elle suit une trajectoire balistique 2 et 
parcourt une distance vt pendant le temps typique r tel que navr = 1 où n est la 
densité de molécules. Cela donne un temps typique de parcours libre (sans collision) de 
l’ordre de : 

r ~-~ 10~ 9 s. 

nav 

On peut en déduire la distance typique parcourue entre deux collisions, appelée libre 
parcours moyen : i = vt ~ ICC 6 m (on notera que l 3> a si bien que notre approche 
est raisonnable). On s’attend donc à un nombre gigantesque de collisions (au moins un 
milliard) aux échelles de temps auxquelles on fait les mesures (au moins une seconde). 

Au niveau microscopique, le mouvement d’une molécule est très chaotique, elle suit 
une marche aléatoire à travers le récipient. 

3 Bilan et objectifs du cours 

On a dégagé les idées suivantes : 

• il y a un grand nombre de particules, il est inconcevable de résoudre les équations 
de Newton pour les N molécules ; 

• le mouvement chaotique des molécules fournit une source d’aléas qui favorise une 
description probabiliste, un peu comme dans le jeu de pile ou face. 

On va donc adopter une démarche probabiliste pour laquelle il faudra : 

• déterminer les évènements pertinents, c’est-à-dire donner une description micro¬ 
scopique des différents états du système ; 

• attribuer une certaine probabilité à ces évènements, compte tenu des contraintes 
globales (ou macroscopiques) sur la base d’arguments généraux et simples, un 
peu comme pour le jeu de pile ou face. 

Une fois ces deux ingrédients définis, on pourra calculer des valeurs moyennes et des 
écarts quadratiques et on verra que ceux-ci sont reliés au comportement thermodyna¬ 
mique du système. Comme on l’a vu, la méthode sera efficace grâce aux simplifications 
qui se produisent lorsqu’on travaille avec un grand nombre de particules. 


2. C’est à dire droite et à vitesse constante. 
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4- Exercices 


4 Exercices 

a Marche au hasard (approche par la combinatoire) 

Un marcheur ivre se déplace sur une droite en sautant aléatoirement à gauche ou 
à droite (de manière équiprobable) à une distance a de là où il se trouve. On prend 
comme origine des coordonnées l’endroit où il se trouve à l’instant t = 0. La durée qui 
s’écoule entre chacun de ses sauts est notée r. 

i Quelle est la probabilité P(n; N) que le marcheur ait, après N sauts, effectué n 
sauts vers la droite et N — n vers la gauche ? Quelle sera alors sa position x en 
fonction de N, n et a ? Quel temps t lui faut-il pour faire N sauts ? 

ii Quel est le nombre moyen de pas vers la droite? 

iii Exprimer la variance de la variable x, définie par (Ax) 2 = (x 2 ) —(x) 2 , en fonction 
de A et a, puis en fonction de t, r et a. 

iv On s’intéresse à présent au comportement de P{n\ N) lorsque n et N sont très 
grands devant 1, mais que Ç = n/N ~ 1/2, c’est à dire £ — 1/2 <C 1. En utilisant 
la formule de Stirling, montrer que P(n\N ) tend vers une loi gaussienne de la 
variable £. 

v Comment varie la distance typique atteinte en fonction du temps ? 

b Distribution des vitesses d’un gaz 

On considère un gaz dont les N particules ont une vitesse notée v = (v x ,v y ,v z ). On 
suppose que les vitesses des particules sont indépendantes entre elles et que, pour chaque 
particule, les variables aléatoires v x ,v y ,v z le sont également. On admet connu que la 
distribution de probabilité de vitesse s’écrit comme le produit de fonctions gaussiennes 
f(v x ,v y ,v z ) = g{v x )g{y y )g(v z ) avec g(x) = Ae~ Bx . 

i Calculer la valeur moyenne E c j n de l’énergie cinétique. En utilisant le résultat 
sur les gaz parfaits E c [ n = ‘-^Nk^T. récrire les paramètres de la distribution de 
probabilité de vitesse en fonction des paramètres physiques. 

ii Quel est l’ordre de grandeur typique de v pour une molécule d’air de la pièce ? 

iii Calculer la valeur moyenne de la norme de la vitesse. Déterminer la distribu¬ 
tion de probabilité y associée et la tracer qualitativement. Calculer son écart 
quadratique. Commenter. 

c Probabilités conditionnelles 

On tire sans remise deux boules successivement dans une urne contenant trois boules 
blanches et deux boules noires. Quelle est la probabilité qu’elles soient noires? 

d Événements indépendants 

On lance deux fois de suite un dé non pipé à six faces. Les événements suivants 

■ A : la somme des deux lancers est 7, 

■ B : la somme des deux lancers est 6, 

' C : le premier lancer donne 4, 

sont-ils indépendants ? 
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Chapitre II Description microscopique et 
entropie statistique 


1 Micro-états et macro-états 

Il est nécessaire de distinguer les états microscopiques (dits micro-états ) et les états 
macroscopiques (dits m.acro-états). 

► Un micro-état est la donnée de grandeurs physiques, fournies par la mécanique clas¬ 
sique ou quantique et définies à l’échelle microscopique. Un macro-état est la donnée de 
grandeurs physiques macroscopiques, qui suffisent à caractériser complètement le sys¬ 
tème étudié à l’échelle macroscopique. 

Intéressons-nous par exemple à un gaz constitué de N particules. Dans l’hypothèse 
où ces particules suivent les lois de la mécanique classique, un micro-état particulier 
est défini par l’ensemble des positions r*i, ■ ■ ■, rjv e t des impulsions p\, ■ ■ ■,piy de toutes 
les particules. Pour un volume V donné, il existe une infinité de tels micro-états. Un 
macro-état caractérisant ce même gaz consistera par exemple en la donnée du triplet 
(N, V, T) si le gaz est confiné dans un volume V fixé et que la température T est imposée 
par un thermostat. Un gaz dont la pression p serait maintenue constante serait quant à 
lui caractérisé par le macro-état ( N,p,T ). 

► Dans certains cas, l’utilisation de la mécanique classique conduit à des résultats trop 
approximatifs et il est plus pertinent alors d’utiliser la mécanique quantique : on doit ré¬ 
soudre l’équation de Schrôdinger, dont la solution (cf. le cours de mécanique quantique) 
fournit un ensemble de niveaux d’énergie accessibles aux particules. Chaque niveau 
d’énergie est occupé par un certain nombre de particules qui définit le tau d’occupation 
du niveau. 1 

Considérons l’exemple suivant : Soit un système quantique constitué de quatre par¬ 
ticules sans interaction, qu’on notera 1, 2, 3 et 4, qui peuvent se répartir sur deux 
niveaux d’énergie ei et e 2 . Considérons par exemple les deux états représentés sur la 
figure 1 . Ces deux états ne diffèrent que par la permutation des particules discernables 
3 et 4 entre les deux niveaux d’énergie accessibles : ils constituent un exemple de deux 
micro-états différents. Le détail de cette répartition n’est cependant pas accessible à un 
niveau macroscopique. Au contraire, on notera que ces deux micro-états possèdent la 
même énergie totale E = 3ei + 62 - Ils constituent un seul macro-état. Autrement dit, un 

1. Le tau d’occupation est défini par les physiciens dans un contexte quantique. Selon la définition 
la plus adéquate, c’est le nombre de particules dans un niveau d’énergie divisé par la dégénérescence de 
ce niveau. Les chimistes divisent plutôt par le nombre total d’états = le nombre de niveaux pondérés 
par leur dégénérescence. 
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FIGURE 1 - Exemple de deux micro-états correspondant au même macro-état. 


macro-état est entièrement déterminé par son énergie totale et les taux d’occupation des 
niveaux (ici, respectivement, 1 et 3). Le schéma de la Fig. 2) suivante répertorie tous les 
macro-états possibles pour quatre particules discernables distribuées sur deux niveaux 
d’énergie, ainsi que le nombre d’états microscopiques différents (les cercles représentent 
une des quatre particules possibles). 

— oooo - —eee- —ee- —e- - 


- -e- -ee- -eee- — oooo - 

Ci = 1 Ci = 4 Ci = 6 Ci = 4 Ci = 1 

FIGURE 2 - Macro-états de différentes énergies totales, ainsi que le nombre de micro¬ 
états correspondants. 

On notera qu’il existe en général plusieurs micro-états pour un même macro-état. 
Comme on le verra à la section suivante, quand on suppose que tous les micro-états 
ont la même probabilité, le macro-état le plus probable (celui qui est observé à l’échelle 
macroscopique) est celui qui correspond au plus grand nombre de micro-états. 

2 Espace des phases et densité d’états 

a Définition de l’espace des phases 

Il est connu que la dynamique d’une particule est parfaitement déterminée dans le 
cadre classique, dès lors que sa position et son impulsion (ou sa vitesse) sont connues à 
un instant particulier. Le passé ou le futur de la particule peuvent ensuite être obtenu 
-au moins en principe- par intégration des équations d’Hamilton (ou de Newton). 

Soit un système constitué d’une seule particule. Le nombre de paramètres indépen¬ 
dants nécessaires pour préciser la position de la particule à chaque instant constitue les 
degrés de liberté du système. Ce nombre est 6 dans le cas d’une particule libre et 
l’état du système à chaque instant est déterminé par un point défini dans un espace 
euclidien à 6 dimensions (3 pour la position et 3 pour l’impulsion) : l’espace des phases 
du système constitué de la particule libre. Dans le cas d’un système constitué de N par¬ 
ticules, la dimension de l’espace des phases sera donnée par le nombre de coordonnées 
généralisées qi et d’impulsions généralisées associées p^. Le nombre de degrés de liberté 
est 6N dans le cas où les N particules sont libres et sans interaction entre elles (< 6N 
dans le cas de contraintes sur le système). Dans cet espace des phases, aussi appelé 
ensemble de Gibbs ou espace T, introduit par Gibbs en 1901, l’état du système à un 
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2. Espace des phases et densité d’états 


temps t quelconque est complètement défini par la position d’un point dans l’espace 
des phases. Au cours du temps ce point représentatif du système décrit une trajectoire 
dans l’espace des phases. 2 

b Densité d’états dans l’espace des phases 

Considérons un gaz composé de N particules ponctuelles. La dimension de l’espace 
des phases est 6N. La distribution dans l’espace des phases est en général continue, 
on la décrit par une densité des états p(p,q,t) à chaque instant. On utilisera sou¬ 
vent la notation p = {pj} et q = {qj}- Le nombre dTl de points représentatifs (i.e. le 
nombre de micro-états) contenus dans le volume de l’espace des phases dT = d 3N qd 3N p = 
dxidyidzidx2 ■ ■ ■ dzjsrdp x idp y idp z idp x 2 ■ ■ ■ dp z jy au temps t et positionné autour du point 
(p, q ) est donné par 

dQ(T,t) = p(p,q,t) d 3N qd 3N p . 

L’espace des phases est complètement spécifié par la donnée de la densité p(p,q,t). La 
probabilité dV que le système occupe un point T = ( p,q ) = (pj, • • • ,Pn, <fi, • • •, (Jn) 
contenu dans le volume de l’espace des phases dT à l’instant t s’écrit : 

~,/ r .x _ p(p,q,t)d 3N qd 3N p 
’ f p(p',q',t)d 3N q'd 3N p' ’ 

et la densité de probabilité dans l’espace des phases : 

dV{T,t) = p(p, q, t) 
dT fp(p>,q>,t)dr ' 


Remarques 

• Les éléments d’un ensemble de Gibbs n’interagissent pas entre eux. 

• Il existe une infinité de points formant l’espace des phases d’un système classique. 

• Les niveaux d’énergie et les états propres d’un système décrivant une ou plusieurs 
particules quantiques (particule dans un puits, ...) sont discrets; dans ce cas, 
l’ensemble de Gibbs est discret et la loi de probabilité associée aussi. 

c Exemples d’espace des phases 

a Oscillateur harmonique à une dimension et à une particule 

On considère l’hamiltonien iL 0 hi à une dimension 

TT P " , 1 2 2 

Hohi = x- V -mu x , 

2m 2 

où m est la masse de la particule et oj la pulsation du potentiel. L’espace des phases est 
de dimension deux et les courbes d’énergie classique iL 0 hi(x,p) = £ sont des ellipses de 

petit et grand axes a = \J et b = \J2me, voir Fig. 3. 

2. On résout les équations canoniques dites de Hamilton : qj = pj = avec j — 1, • • •, M, 

dans un espace des phases de dimension 2 M. L’intégration des 2 M équations précédentes, connaissant 
la position {• • • ,p°, • • •, • • •} du point dans l’espace T à un instant to donné, détermine complètement 

l’ensemble {• • • ,Pj(t), • • •, qj (t), • • •} i.e. la trajectoire dans cet espace. 
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Figure 3 - Représentation de l’espace des phases de dimension deux (une dimension 
pour x et une pour p x ) représentant une particule dans un oscillateur harmonique 
unidimensionnel. 

(3 Particules classiques identiques et indépendantes à 3 dimensions 

Soient N particules classiques indépendantes, on note qi le vecteur position de la 
particules i (les composantes de (p sont ( Xi,yi,Zi )) et p t son vecteur impulsion (i.e. sa 
quantité de mouvement). L’hamiltonien total du système à N particules peut s’écrire 

N 

H(qi,pu ■ ■ ■ ,qn,pn) = ^) > 

i=l 

où l’hamiltonien à une particule h est le même pour toutes les particules. Dans ce cas, 
la dimension de l’espace est 6 N. 

► Un cas particulier que l’on rencontrera très souvent est celui du gaz parfait 

N 9 

ip ~ 

^{qiiPi,--- Ân,pn) = 2_j ^ > 


où pi = \\pi\\ et m est la masse des particules. 

► À l’inverse, si on inclut des interactions entre particules, 'H n’est plus réductible à une 
somme de termes indépendants mais la dimension de l’espace des phases est inchangée. 

d Surface d’énergie constante 

Les particules n’occupent pas tout l’espace des phases mais respectent au contraire 
l’équation maîtresse : 

H(qi,pi,--- ,cin,Pn) = E , 

où E est l’énergie totale du système, qui est conservée au cours du temps si le système 
est isolé. Cette équation définit une hypersurface de l’espace des phases. 

Dans l’exemple précédent des N particules dans un espace de dimension 3, la surface 
d’énergie constante a pour dimension 61V — 1. 

e Particules indiscernables 

Pour des particules indiscernables, il faut diviser par N\ qui compte le nombre de 
permutations des configurations microscopiques que l’on peut faire sans modifier l’état 
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2. Espace des phases et densité d’états 


physique. Ainsi, dans le cadre classique, les formules précédemment définies (fonctions 
de partition, moyennes, ..) devront être interprétées avec la règle de substitution: 



m 


d 3 q N d 3 p N 


(1) 


Au bilan, on utilisera les règles de dénombrement des micro-états suivantes : 


• Chaque point de l’espace des phases occupe un volume élémentaire h a où a 
est le nombre de couples bi-scalaires ( Çj,Pj ) dans l’espace des phases. 

• Si les particules sont classiques, identiques et indiscernables, on divise l’élé¬ 
ment de volume de l’espace des phases par Ai! 



ol Moyenne d’ensemble 

Soit G(p, q) une grandeur physique quelconque, sa moyenne d’ensemble à un instant 
quelconque t est définie par : 

L G(p,q ) p(p,q) dT 

= JV T wvu — ( 2 ) 

J Vr p(p, Q) dT 

où Vp représente le volume de l’espace des phases et p et q désignent {pi}, {%}• 

j3 Hypothèse ergodique 

L’approche de Gibbs repose sur l’hypothèse suivante : pour les ensembles conve¬ 
nablement choisis, la moyenne d’ensemble (2) est identique à la moyenne temporelle, 
celle-ci étant alors effectuée sur le système : 


(G) = G, 

_ 1 l'to+A.t 

où G = A lim Â7 / G(p(t),q(t))dt 
At —»oo Ai J to 

et les variables canoniques évoluent au cours du temps, p,q = p(t),q(t). Si le système 
est à équilibre, la limite ne dépend pas de tç,. 

Il est souvent plus facile de calculer la moyenne d’ensemble que la moyenne tempo¬ 
relle G(t), en particulier pour des systèmes hors équilibre. 
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g Variables extensives et intensives 

Parmi les variables macroscopiques introduites précédemment, certaines sont dites 
extensives, d’autres intensives. Pour définir une quantité extensive ou intensive, on peut 
faire l’expérience de pensée suivante : on double le système étudié, de façon à ce qu’il 
garde les mêmes propriétés physiques ; les quantités extensives sont doublées lors de 
l’expérience, les quantités intensives sont celles qui sont conservées. 

Sont intensives la température, la densité, le potentiel chimique, la vitesse de dérive, 
l’énergie d’agitation par particule, la pression... Les quantités extensives sont plus rares : 
l’énergie totale du système (souvent notée U ou E), V le volume qu’il occupe, N le 
nombre de particules total (Ni, N 2 , ..., les nombres de particules de type Ai, A 2 , ..., 
quand il y a différentes espèces de particules), la capacité calorifique. Tout rapport entre 
deux quantités extensives redonne une variable intensive, par exemple E/N l’énergie par 
particule, ou N/V la densité volumique. 


3 Entropie statistique 


a Information manquante 

Par rapport à la certitude absolue, une description probabiliste comporte intrinsè¬ 
quement un certain manque d’information. Ce défaut d’information peut être mesuré 
par une notion empruntée à la théorie de l’information et que l’on appelle l’entropie 
statistique associée à une loi de probabilité. 

La théorie de l’information établit comment mesurer quantitativement l’information 
a priori que l’on peut avoir sur un système donné à partir de l’entropie statistique. 

b Définition de l’entropie statistique 

Soit {e m }, {m G 1, • • •, P}, un ensemble discret d’événements possibles auquel est 
associé l’ensemble { P m } des probabilités respectives vérifiant les relations : 


0 < P m < 1 , V m = 1..P , 

a 

Y P ™ = 1 • 

m= 1 


(3a) 

(3b) 


• Si pour un événement e mQ on a P mQ = 1, cet événement certain exclut toutes 
les autres possibilités (P m . Q = 0, ra / m 0 ) : l’information est a priori complète 
sur l’état du système considéré (ou sur le résultat d’une expérience relative aux 
P événements considérés effectuée sur le système). 

• Toute autre distribution des probabilités {Pm} introduit un manque d’infor¬ 
mation par rapport au cas précédent. 

Ce manque d’information est mesuré par l’entropie statistique définie selon le cas 
d’une distribution discrète ou continue par : 
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• Cas discret : 

(4a) 


• Cas continu : 

(4b) 


où V(x) représente la densité de probabilité de la variable aléatoire continue x, 
telle que f R P(x ) dx = 1 . 

► Dans ces définitions, 3 k est une constante positive et a priori quelconque. Elle est 
fixée par les considérations spécifiques suivantes : 

• En Physique Statistique on choisit 


S(V) = —k / V(x)lnV{x)dx 



£7 

S(Pl,P2,'- 

' , Pü) = -k Ÿ2 Prn lu Pm 


m= 1 


k = k~B \k~B = 1,38 1 0 23 J/K] , 

où est la constante de Boltzmann: k-Q = rapport de la constante molaire 
des gaz parfaits classiques au nombre d’Avogadro. 

Ceci permet (voir plus loin dans le cours) d’identifier l’entropie statistique à 
l’entropie classique thermodynamique lorsque le système est à l’équilibre. 

• En théorie de l’information on choisit plutôt 


Par exemple, si on considère deux événements équiprobables (Pi = P -2 = 1/2), 
avec P = 2 , on trouve : 


S 



= 1 , 


l’unité étant le bit. 


Le manque d’information pour ce système à deux événements est maximum et 
égal à 1 bit dans le cas d’équiprobabilité. 

c Propriétés 

Les propriétés suivantes montrent que la définition donnée de S convient pour la 
mesure du manque d’information concernant un système. 

• S est définie positive: pour tout {P m } vérifiant les conditions (3a) et (3b), 


S > 0 . 


• S est totalement symétrique : toute permutation des P m laisse S inchangée : 


S (Pl, P 2 , • • •, Pçi) = S (P 4 , Pl, •••) = ••• 


3. La définition de l’entropie statistique (4a) est complétée par Pm ln P m = 0 si P m = 0 puisque 
xlnæ —> 0 lorsque x —» 0. 
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• S est additive : considérons deux ensembles indépendants d’événements pos¬ 
sibles {e m }, m = 1..M et m,' = 1 ..AI' auxquels correspondent les deux 

ensembles de probabilités {P m } et {P' m /}, vérifiant les conditions générales (3a) 
et (3b). Une expérience “double” dans laquelle se produisent un événement e m 
et un événement e 1 , a pour probabilité : 

Pm,m! = Pm ' Pm' (5a) 

en appliquant (A 2) car les deux séries d’évènements sont indépendants. D’après 
la définition (4a), l’entropie statistique associée à cette expérience double est 

M M 1 

s ({p m ,m r }) = —k Pm,m' In Pm ,ra' (5b) 

m= 1 m '=1 

et, en reportant (5a) dans (5b), on obtient : 

S{{P m ,m'}) = S ({P m }) + S ({P ml }) . (5c) 

Le manque d’information dans l’expérience double est ici la somme des manques 
d’information sur chacun des deux ensembles indépendants. 

• L’entropie est extensive : c’est lié à l’additivité de S et on le revérifiera ultérieu¬ 
rement sur des exemples physiques. 

• Minimum de S : l’entropie statistique s’annule lorsque l’une des probabilités vé¬ 
rifie P mQ = 1 : 

5(0, 0 , 0 , •••, 1 , 0 , 0 - - -, 0) = 0 . 

Cette valeur est un minimum absolu, ce qui signifie que l’information est complète 
(manque d’information nul). 

• Maximum de S : comme P m est une probabilité avec 0 < P rn < 1, l’entropie 
statistique est une fonction positive ou nulle, à concavité tournée vers le bas : 



Le maximum de S correspond à la configuration d’équilibre thermodynamique. 
► Quand il n’existe aucune contrainte physique imposées au système, notamment 
s’il est isolé, 4 ce minimum est réalisé pour 

p =1 

m n ’ 

où U est le nombre de configurations. On verra que cette situation correspond à 
l’ensemble d’équilibre microcanonique. 

4. Un système isolé n’échange ni matière ni énergie avec l’extérieur. 
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4- Postulat d’entropie maximale 


d De l’entropie statistique à l’entropie thermodynamique 

L’entropie définie en (4a), vérifiant les propriétés développées au paragraphe pré¬ 
cédent, possède de ce fait les propriétés requises par l’entropie thermodynamique 
classique à l’équilibre : 

i La relation (5c) établit que S est une fonction additive d’ensemble ; propriété 
partagée avec l’entropie thermodynamique classique. 

ii S présente un maximum dans le cas où tous les événements sont équiprobables 
(cas du système isolé). Ce sera l’objet de la physique statistique de vérifier que, 
dans ce cas, les événements équiprobables ont lieu à l’équilibre et que l’entropie 
est alors maximale. Ceci correspond entièrement au deuxième postulat de la 
thermodynamique. 

iii L’entropie statistique S, définie par la relation (4a), s’identifie à l’entropie ther¬ 
modynamique classique avec le choix k = ks pour les systèmes à l’équilibre. 
Attention, cependant, que les notions d’entropie statistique S ainsi que d’entropie 
thermodynamique peuvent être prolongées pour des systèmes hors équilibre. 
Dans ce cas, elles ne coïncident pas. 

4 Postulat d’entropie maximale 

Le principe même de la démarche probabiliste adoptée en physique statistique 
conduit à affecter à chaque état microscopique -qu’on notera génériquement m- une 
certaine probabilité P m de concourir à la réalisation d’un état macroscopique M par¬ 
ticulier. Le problème central de la physique statistique est de déterminer cette loi de 
probabilité. D’une façon générale, lorsque le système est à l’équilibre, il n’y a, a priori, 
aucune raison de penser qu’il y ait une quelconque différence de réalisations entre les 
différents micro-états accessibles, ce qui doit conduire à une entropie maximale puisque 
le défaut d’information est alors le plus grand. Cette conviction est érigée en un prin¬ 
cipe, appelé principe d’entropie maximale qui sert d’unique fondement à toute la 
physique statistique des systèmes à l’équilibre. Il s’énonce ainsi : 

Un état macroscopique est à l’équilibre lorsque son entropie statistique 
est maximale compte tenu des contraintes imposées au système. 

Dit d’une autre façon encore: un état hors d’équilibre évoluera jusqu’à ce qu’il ne 
contienne pas plus d’informations que celles associées aux contraintes qui s’exercent sur 
lui. La validité de la Physique Statistique repose sur le fait qu’à la limite thermodyna¬ 
mique : 

- tous les ensembles d’équilibre sont équivalents, 

- à chacun d’eux est associée une fonction macroscopique identifiable à un 
potentiel thermodynamique. 
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5 Exercices 

a Micro-états et macro-états 

Un système physique composé de deux particules compte trois niveaux d’énergie 

équidistants (eo = 0, e\ = e, 62 = 2e) sur lesquels peuvent se répartir ces deux particules. 

i On se place d’abord dans le cas de particules discernables et de niveaux d’énergie 
non dégénérés. 

• Combien il y a-t-il d’états macroscopiques d’énergie totale E = 2e? 

N.B. : on les distingue selon le tau d’occupation de leurs niveaux d’énergie. 

• Combien il y a-t-il de configurations microscopiques possibles pour chacun 
des états macroscopiques déterminés précédemment ? 

• En déduire l’état macroscopique le plus probable. 

ii Reprendre le problème dans le cas de particules discernables à répartir sur les 
niveaux d’énergie 60,61,62, de dégénérescences, respectivement, 1,2,1. 

iii Reprendre le problème dans le cas de particules indiscernables à répartir sur 
les niveaux d’énergie 60,61,62 d’abord non dégénérés, puis de dégénérescences, 
respectivement, 1,2,1. 

b Oscillateur harmonique à une dimension 

Un oscillateur harmonique classique unidimensionnel a pour hamiltonien 

Tj, i P 2 , 1 22 

H \ x ,P) = 7 t -t -mu x , 

2m . 2 

où uj représente la pulsation et m la masse de la particule qui vibre sur un axe. 

i Quelle est la dimension de l’espace des phases ? 

ii Donner la trajectoire du système ainsi que son équation pour une énergie e fixée. 

iii Connaissant la condition initiale où la position et l’impulsion vérifient 
x(t = 0) = xo, p(t = 0) = 0, exprimer l’énergie fixée e en fonction de xq. 

iv Calculer U le nombre de micro-états accessibles à l’oscillateur harmonique clas¬ 
sique d’énergie comprise entre e et e + 5e. 

v Comparer ce résultat au calcul quantique : les énergies d’un oscillateur harmo¬ 
nique quantique à une dimension sont données par e n = hio(n + 1/2). 

vi On fait l’approximation suivante : on calcule D non plus sur la couronne définie 
par e < H(x,p) < e + 5e mais sur toute l’aire H(x,p) < e. Recalculer U. 

vii Déduire de la dernière expression de D l’entropie S. 

viii On montrera au chapitre suivant la formule 1/T = dS/de\ calculer la tempéra¬ 
ture T et l’exprimer en fonction de l’énergie e. 
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chapitre ni L'ensemble microcanonique 


1 Postulat fondamental 

a Système isolé 



On considère un gaz, un fluide ou un solide, composé de particules. Si elles sont 
toutes identiques, on dit que le fluide ou le gaz est simple, ou que le solide est pur. Sinon, 
c’est un système à plusieurs composants. On étudie le système formé par ces particules 
à l’équilibre, considéré comme isolé, c’est à dire n’échangeant ni énergie, ni chaleur, 
ni particules avec l’extérieur. On supposera de plus qu’il est placé dans une enceinte 
rigide. On travaille alors dans l’ensemble microcanonique d’équilibre thermodynamique 
(ou plus simplement ensemble microcanonique). Les contraintes SU = 0 et SQ = 0 
imposent 1 SW = 0, autrement dit aucune force ne travaille. 

Les paramètres macroscopiques suivants sont donc fixés : E l’énergie totale du sys¬ 
tème, N le nombre de particule (plus précisément, s’il existe N\ particules de l’espèce 
A\, N 2 de l’espèce A 2 , et ainsi de suite, les nombres N t sont conservés composant par 
composant) et V le volume de l’enceinte. De plus, on montrera que l’entropie S est 
également constante dans l’ensemble microcanonique. 

On peut toutefois généraliser ces définitions et inclure le cas où des forces travaillent, 
par exemple un champ de force extérieur agissant sur les particules, à condition de 
compter les énergies d’interaction dans le bilan. On perd alors la contrainte SU = 0 

1. Voir la section de thermodynamique ci-après. 
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et on ne garde que ôQ = 0. On peut montrer dans ce cas que l’entropie S est le seul 
paramètre macroscopique, mais nous ne l’étudierons pas. 

b Application du principe d’entropie maximale 
a Multiplicateur de Lagrange 

On calcule l’entropie thermodynamique à partir de l’entropie statistique donnée par 
la formule (4a) avec k = k B . On va appliquer le principe d’entropie maximale pour la 
déterminer. Remarquons d’abord que l’entropie statistique est effectivement maximale 
pour la distribution uniforme P m =c.onst.= 1 /fl, Vm, avec fl le nombre de configura¬ 
tions, lorsque le système n’est soumis à aucune contrainte. Le problème se ramène à dé¬ 
terminer le maximum de la fonctionnelle s[P ] = s[P\, • • •, Pq\ = — k B X^m=i lu Pm , 
compte tenu de la contrainte de normalisation des probabilités y) . P m = 1. Ce 
problème d’extrémalisation avec contraintes peut être appréhendé par la méthode des 
multiplicateurs de Lagrange (voir annexe D). Il faut donc déterminer la solution de 
l’équation 2 

n / n 

— ks ^2 Pm ln Pm + Ao f ^ -Pm - 1 

m =1 \m= 1 

où Ao désigne le multiplicateur associé à la contrainte. Les équations d’extrémalisation 
de S s’écrivent 

dS_ 

dpi 

dS 
d\ o 

(3 Équiprobabilitê 

La relation (6b) qui exprime que S est extrémale par rapport à Ao redonne exacte¬ 
ment la relation (3b) qu’on avait écartée. La relation (6a) donne 

b _ P —i+-Wfc B 

1 m — c 5 

ce qui n’est guère parlant. Cependant, Ao est déterminé par la relation Y2L=i Pm = 1, 
qui implique la relation fl e -1+ ^°/ fc B = 1. La distribution d’entropie statistique maximale 
devient alors 

Pm= Q- (7) 

On constate que « les probabilités de chaque configuration sont égales quand on se place 
dans l’ensemble microcanonique. » L’ensemble microcanonique correspond à l’ensemble 
probabiliste d’information minimale et toutes les configurations y sont équivalentes. 

Il ne faut pas se méprendre sur l’interprétation de ce principe. Considérons par 
exemple une configuration totalement improbable où les particules sont confinées dans 
un coin du volume : 

2. En toute rigueur, la condition qui suit ne permet que de déterminer des extrema de S (maximum 
ou minimum). Il faudrait encore montrer que la distribution trouvée est un maximum , ce qui n’est pas 
vraiment nécessaire puisque l’entropie statistique est concave. 
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1. Postulat fondamental 



Cette configuration a autant de chance de se produire que la configuration plus réaliste 
représentée au début de la section ! ce qui semble contre-intuitif. Il faut considérer l’en¬ 
semble des configurations irréalistes et les comparer avec l’ensemble des configurations 
réalistes : le rapport entre le cardinal du premier et celui du second est un nombre infini¬ 
ment petit, de sorte que, globalement, les premières n’ont aucune chance d’être réalisées, 
comparées aux secondes. La statistique sur les nombres de configurations contient toute 
l’information pertinente et compense l’absence apparente d’information qui s’exprime 
par l’équiprobabilité de chaque configuration. 

7 Fonction de partition 

Le nombre de configuration est appelé fonction de partition de l’ensemble microca¬ 
nonique et noté 12 (ou parfois Z mc ). 

ô Formule de Boltzmann de l’entropie 

On va maintenant retrouver la formule de Boltzmann, qui permet d’exprimer l’entro¬ 
pie à l’aide de 17. On réinjecte la formule (7) dans la définition de l’entropie statistique, 
ce qui donne : 


S = -l 


kB Y, Pi lo ë(^) = kB l °ëi^)Y Pi = kB log (^) 


( 8 ) 


iex 


iei 


=i 


Cette dernière formule est la définition de l’entropie thermodynamique, proposée par 
Boltzmann. À l’équilibre thermodynamique, l’entropie statistique et l’entropie thermo¬ 
dynamique se confondent ; attention, ceci ne reste valable qu’au voisinage de l’équilibre. 

L’entropie mesure le désordre d’un système. Le cas où un système ne possède qu’un 
seul état (quand tous les paramètres sont gelés) correspond à S = 0, c’est l’ordre absolu. 
On peut remarquer que le multiplicateur de Lagrange s’écrit 


Ao = -kn log(fi) + k B , 


soit Aq = —S + k B - On vérifie que Àq est un paramètre physique macroscopique. 


23 





Chapitre III. L’ensemble microcanonique 


c Densité d’états 


ol Calcul de la fonction de partition 

On calcule la fonction de partition avec les variables conjuguées de position et d’im¬ 
pulsion, en passant au continu. Plaçons-nous, pour clarifier les idées, à trois dimensions 
d’espace. On note toujours dT = d 3 q\d 3 p\ ■ ■ ■ d 3 qjyd 3 pN et on obtient alors 


n = 


1 

hW 



(9) 


où on a précisé le domaine d’intégration Vp = {'H = E, V ) conformément aux contraintes 
dans l’ensemble microcanonique. 


(3 Densité d’états intrinsèque 

Par définition, la densité d’états est le facteur qui provient du passage d’une somme 
discrète à une intégrale. Si on choisit des variables conjuguées et qu’on travaille dans 
l’ensemble microcanonique, on voit, dans l’expression précédente, qu’il vaut 

1 

L’expression (10) de p doit parfois être reformulée dans des cas particuliers. Toute¬ 
fois, l’illustration du résultat fondamental de l’ensemble microcanonique demeure tou¬ 
jours : la densité d’états pour des variables conjuguées est constante. 



7 Facteur d’indiscernabilité 

La formule (9), adjointe à (8), met déjà en relief l’extensivité de l’entropie. Cepen¬ 
dant, il n’a pas été tenu compte de l’indiscernabilité des particules. Ce problème a été 
soulevé puis résolu par Gibbs à la fin du 19 cme siècle. On a vu au chapitre précédent 
qu’il faut diviser l’expression de O par N\ pour tenir compte de toutes les permutations 
possibles entre particules. 3 Précisons les choses. 

Quand on effectue une permutation, on passe d’un état à un autre état, ces deux 
états ont été comptés séparément dans fî. Pourtant, ils sont physiquement identiques. 
Le génie des physiciens a été de comprendre qu’il fallait ne pas les compter séparément. 

De façon plus générale, quand on fait agir le groupe des permutations sur les confi¬ 
gurations, on peut montrer que les classes d’équivalences ont toutes exactement Ai! 
éléments, qui doivent être identifiés, ce qui explique qu’il suffit exactement de diviser 
par ce facteur pour trouver le bon nombre de classes. Cet argument n’est valable que 
dans une approche classique continue. Quand des observables, notamment l’énergie, 
sont quantifiées, on vérifiera, dans la partie sur les calculs statistiques semi-quantiques, 
que ce facteur n’est correct que pour des fermions. 


2 Exemple de calcul de la fonction de partition 


a Le cas du gaz parfait 

Un gaz est considéré parfait lorsque les particules n’interagissent pas entre elles : 
l’énergie potentielle d’interaction entre les particules U (ri, rj ) =0 est nulle V?y, Tj. Dans 

3. Ceci ne change pas la densité (10) calculée précédemment car le facteur 1/N\ ne participe pas 
formellement à la densité mais il faut bien le prendre en compte. 
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2. Exemple de calcul de la fonction de partition 


ce cas, l’hamiltonien se réduit à la contribution cinétique de chacune des particules : 
fi = b • Le système des N particules libres de se mouvoir dans le volume V 
est à l’équilibre avec l’énergie E fixée à 5E près (grâce à des parois adiabatiques). 

La fonction de partition microcanonique de ce gaz parfait se calcule d’après (9) : 

N 

m.v.E ) = 1 ^ SW T]I dh.êpi 


N\(h) 

1 


U 


J E<H<E+ÔE,V 


N\(h ) 


3N 


d 3 ri ■ ■ ■ d 3 r N 


'V 


L 


d 3 pi ■ ■ ■ (ï 3 p N . 


E<H<E+SE 


^ — 


=V N 


L’équation ^L- \\pi\\ 2 = E (resp. ^ Yli \\Pi\\ 2 = E + ô) correspond à une hyperboule 

de rayon \J2rnE (resp. y/2m(E + SE)) en dimension 3N (les 3N degrés de liberté des 
impulsions puisque les 3 N degrés de liberté des positions ont déjà été intégrés pour 
donner le facteur V N ). Le volume de l’espace des impulsions accessibles au système 
est donc celui compris entre les deux hyperboules de rayon \J2mE et yj2m.{E + ÔE). 
Le volume d’une hyperboule de rayon R en dimension d étant donné par: Vd(R ) = 
R d (voir annexe C), on en déduit : 


r(f+i) 


n(N,V,E) = 


V 


N 


( 

7T 2 


N\h 3N r(^ + l) 

yN _ {n) 3N/2 (2 mE) 3N / 2 

mh? N r(^ + i) 1 j 

V N (2vr mE) 3N / 2 3 N ÔE 

N\h 3N rïWTï) — ~Ë ' 


m{E + ÔE) 3N - V2mE 3N 


1 + 


3 NÔE 
~2E 


- 1 


r(^ + i) 

b Calcul des grandeurs thermodynamiques 
a Variante simple du calcul de 17 

On a accès aux grandeurs physiques (température, pression et potentiel chimique 
microcanoniques) en calculant directement le nombre de micro-états T pour lesquels 
fi < E (cette expression est équivalente à la précédente et la remplacera dorénavant) 

V N 

Ü{N,V,E) = 


, 3N 

(27 xmE) 2 


N\h 3N r(^ + l) ' 

On utilise l’approximation de Stirling (voir annexe C) qui conduit à 

3 N 


ln 


r(f + D 


3iV 3 N 


pour obtenir le logarithme de la fonction de partition et, partant, l’entropie du système : 

/ \ 


ln Çî(V,E,N) = N 


5 . V 3, , AirmE 1 

— ln — ^ — ln (-—— ) -)- — 

2 N 2 V 3 h 2 N ’ N 


, 3 N , 5E\ 

ln-h ln- 

2 EJ 


termes provenant de ^expression » 
première de et négligeables / 
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Chapitre III. L’ensemble microcanonique 


f3 Grandeurs thermodynamiques 

En définissant la densité n = y et l’énergie par particule e = jj (deux grandeurs 
intensives) et en prenant la limite IV —>• oo, qui correspond à la limite thermodynamique, 
on obtient (les termes éventuels en 1/N , ln N/N ou ÔE/E sont négligeables) : 


S(N,E,V) 


k B ln Ü(N,E,V) 


Nk B 


5 ( 3 1 ( Annie \ 
2 + 2 Ü1 \3(/i 3 n ) 2 / 3 ) 


et on vérifie le caractère extensif de S. 


3 Lois thermodynamiques 

a Premier principe 
a Énonce historique 

Le premier principe de la thermodynamique stipule que la variation d’énergie, qui 
est une différentielle totale, s’écrit : 


dE = 



+ 


échange de chaleur 



travail des forces 


(lia) 


C’est la généralisation aux gaz, fluides et autres systèmes statistiques, du théorème 
de l’énergie cinétique de la mécanique ponctuelle. Le premier terme correspond aux 
transferts de température entre corps et on verra que la température exprime l’agitation 
cinétique des particules, appelée ainsi agitation thermique. 

Le second terme s’identifie directement, mais son calcul est très différent des calculs 
de travaux, tels qu’on les fait en mécanique du point. 

Ces deux termes ne sont pas des différentielles exactes, d’où leur notation. 4 

(3 Chaleur 

Si un système évolue de façon quasi-statique, l’évolution est réversible et on est 
autorisé à écrire : 

ôQ rev = TdS , (11b) 

tandis que, lors d’une réaction irréversible, on a 


5Q ur < TdS 


car dS se décompose en 


dS 


dS e ch + dS c 


SQ , rn 

-je + dS c 


(11c) 


Autrement dit, il existe une partie dite de création de la variation d’entropie, qui ne 
correspond à aucun échange de chaleur. 

4. Contrairement à f dE, l’intégrale f SW dépend du chemin choisi pour l’intégration. 


26 




3. Lois thermodynamiques 


7 Travail des forces extérieures 

Si le système est soumis à une force extérieure, comme en mécanique du point, il lui 
correspond un travail, qui s’écrit de façon analogue. On peut citer la force de gravitation 
(qu’on néglige souvent), un champ électrique ou magnétique. 

Certaines forces sont plus spécialement adaptées à la physique statistique. Par 
exemple, à la force de pression p exercée sur les parois correspond le terme 

ÆWpression = ~p5V , 


où V est le volume. 

S Potentiel chimique 

Le potentiel chimique p est l’énergie potentielle résultant de toutes les interactions 
du milieu que subit une particule, —p est l’énergie qu’il faut fournir pour extraire la 
particule. Spontanément, les particules vont là où leur potentiel chimique est le plus 
bas (là où elles sont le plus stables). On l’appelle, dans d’autres contextes, énergie de 
cohésion. 

Les N a particules d’une même espèce A a voient le même potentiel chimique p a ; à 
une variation 6N a de leur nombre leur correspond un terme de travail 

^fLchimi = P*SN a . 

e Énoncé complet 

On obtient finalement 

dE = TdS-pdV+ ^2p a dN a , (12) 

a 

équation qui pourra être utilisée tout le temps. En effet, quand TdS > ÔQ (cas irré¬ 
versible), il vient également des termes irréversibles dans le travail, ôWi rr > SW Tev . qui 
compensent exactement les termes manquants dans la chaleur. 

b Le second principe retrouvé 

Le principe d’entropie maximale est compatible avec le second principe de la ther¬ 
modynamique. Celui-ci, qui s’applique plus généralement aux systèmes hors équilibre, 
stipule que l’entropie croît toujours, lors d’un processus thermodynamique (elle reste 
constante si le processus est réversible). Ainsi, l’état d’équilibre peut être compris comme 
l’aboutissement d’un processus irréversible particulier de retour à l’équilibre. À l’équi¬ 
libre, l’état du système ne peut plus évoluer et l’entropie, de façon correspondante, ne 
peut plus croître ; elle donc est maximale, comme nous l’avons établi dans les sections 
précédentes. 

c Du microscopique au macroscopique 

On peut maintenant résumer les étapes essentielles de la physique statistique, com¬ 
prise comme le passage du microscopique au macroscopique : 
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Chapitre III. L’ensemble microcanonique 


description microscopique description probabiliste thermodynamique 

' H(qi,pi,...,q N ,PN ) <—> 

fl(E, V, N) ► S(E, V, N) 

postulat 

formule 

fondamental 

de Boltzmann 

+ contraintes 

S = ks ln fl 


En pratique, on essaiera de suivre cette démarche en calculant fl pour un hamil¬ 
tonien donné afin d’obtenir l’équation d’état et d’en déduire la thermodynamique. On 
en a vu un exemple avec le gaz parfait. Au-delà de l’aspect méthodologique et calcula- 
toire que propose cette démarche, la force de l’approche statistique est de donner une 
interprétation simple et une image microscopique de situations thermodynamiques non 
intuitives. On étudiera un exemple élémentaire de création d’entropie irréversible lors 
d’une détente de Joule Gay-Lussac. Mais avant cela, on va montrer comment on obtient 
les grandeurs thermodynamiques à partir de la fonction S(E,V, N). 

d Calcul des variables macroscopiques à l’équilibre thermodynamique 

Pour trouver les quantités thermodynamiques, on va identifier les coefficients de 
cette différentielle avec ceux de la différentielle formelle 



L’énergie est une fonction génératrice, qui permet de calculer toutes les autres quan¬ 
tités thermodynamiques, elle est pour cela appelée potentiel thermodynamique de l’en¬ 
semble microcanonique. 

a Température thermodynamique 

On définit la température thermodynamique à partir de (12) : 


V,N 

1 3 Pression thermodynamique 

De même, on définit la pression thermodynamique à partir de (12) : 


S,N 

7 Potentiel chimique thermodynamique 

De même, on définit le potentiel chimique thermodynamique à partir de (12) : 


s,v 


(13c) 



(13b) 



(13a) 
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3. Lois thermodynamiques 


e Propriétés de l’entropie 
a Extensivité 

L’entropie est extensive, comme déjà annoncé, ce qui s’écrit : 

S(E r + E 2 , Ld + V 2 ,Ni + N 2 ) = S(E U P, Ad) + S(E 2 , V 2 , N 2 ) ; 
S(E, V, N) = N s(e,l/n) . 


(3 Stabilité thermodynamique de l’équilibre 

De plus, elle vérifie 

d 2 S\ 

dËPJ, 


< 0 


Or, la capacité calorifique s’écrit ! 


1 


C v = 


V,N 


d 2 S 


T 2 \ dE 2 


V,N 1 


où on utilise la température T, qui a été calculée plus haut. En conséquence, C v est 
toujours négative. Il faut noter que C v est également une quantité extensive. 

f Équation d’état du gaz parfait microcanonique 

On peut aller plus loin et retrouver l’équation d’état du gaz parfait (on notera 
n = N/V la densité de particules du gaz et e = E/N l’énergie par particule du système). 


KN,E,V) 

T(N,E,V ) 

P(N,E,V) 
T(N,E,V ) 

1 

T(N,E,V) 


d S 3 ( 47 t me \ 

~dN^ N? ’ ) = 11 \3(n/i 3 ) 2 / 3 J 


^(N,E,V) = ^ 

dv y ' ' 1 v 


§ ( N,E,V) = 3 -m B (l 


pV = Nk B T 


E = -Nk B T 


IL#’ Remarques : 

• Les deux dernières relations montrent que le traitement du gaz parfait dans l’en¬ 
semble d’équilibre microcanonique redonne bien les équations d’état et d’équi- 
partition de l’énergie (d’un gaz parfait) à la limite thermodynamique. 

• On vérifie que, pour le gaz parfait, la capacité à volume constant est constante 
et s’écrit : 

C v = | Nk B . (14) 


5. Cette formule n’est jamais indiquée dans les livres, ce qui montre qu’on calcule rarement la 
capacité calorifique dans l’ensemble microcanonique. 
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Chapitre III. L’ensemble microcanonique 
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4- Exercices 


4 Exercices 

a Entropie de mélange ; paradoxe de Gibbs 

On considère un mélange idéal de deux gaz parfaits, c’est-à-dire sans interaction 
entre eux, dont les particules sont différentes et donc discernables entre elles (mais 
indiscernables au sein d’une même espèce). On les indicera par 1 et 2 et on notera N\ 
et N 2 le nombre de particules de chacune des espèces. Les deux gaz sont initialement 
(situation I ) dans deux enceintes de volume, respectivement, V\ et V 2 , séparées par 
une paroi diathermique (à l’équilibre thermique), et libres de se déplacer (à l’équilibre 
mécanique). Le tout est isolé. On retire ensuite la paroi et on laisse les gaz se mélanger 
(situation F) dans le volume V = V\ + V 2 . 

i Dans la situation /, on suppose provisoirement les énergies E\ et E 2 fixées. 
Montrer que le nombre de configuration de l’ensemble 1 U 2 à l’équilibre s’écrit 

Vi,V2,Ni,N 2 ) = T th (-Ei, E 2 , Ni, N 2 ) . 

ii En réalité, l’énergie E\ peut varier entre 0 et E = E\ + E 2 , de sorte que le 
nombre de configurations s’écrit, dans la situation I, 

pJVl y JV2 

Qi(E, Ei, V 2 , N 1} N 2 ) = ^ J T th (E, N U N 2 ) 
avec T th (E,N u N 2 ) = J 0 E f th (E 1 ,E-E 1 ,N 1 ,N 2 )^. 

Calculer les pressions des gaz en fonction de Ni, N 2 , V], V 2 et T (qu’on ne 
calculera pas). 

iii Dans la situation F, justifier que le nombre de micro-états accessibles s’écrit : 

yAà yN 2 

n F (E,V,N u N 2 ) = T t h{E, Ni,N 2 ) . 

iv Montrer que la pression totale s’écrit comme la somme des pressions partielles de 
chacun des gaz (la pression partielle est la pression d’une espèce si elle occupait 
seule l’enceinte). 

v Déterminer la variation d’entropie totale AS = S F — Sj en fonction de Ni, N 2 , 
V\ et V 2 . Interpréter physiquement l’irréversibilité de la transformation I —>• F. 

vi On suppose dorénavant que les gaz ont initialement des pressions égales, p\ = p 2 - 
Donner la relation entre Ni, N 2 , Vi et V 2 . 

vii Donner, dans ce cas, l’expression de AS en fonction des nombres de particules 
Ni et N 2 . 

viii Si les deux gaz sont constituées de particules du même type, qui sont donc 
toutes indiscernables, les pressions initiales étant supposées égales, le précédent 
calcul aboutit à une création d’entropie AS > 0 alors que l’expérience est, dans 
ce cas, réversible. 

Trouver l’expression de Q F (E, V, Ni, N 2 ) valable dans ce cas et calculer la varia¬ 
tion d’entropie correcte. 
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Chapitre III. L’ensemble microcanonique 


b Défauts dans les solides 

Un cristal est constitué par un réseau régulier de N sites. À température nulle, tous 
les sites sont occupés. Lorsque la température augmente, divers défauts apparaissent, 
parmi lesquels des lacunes (défaut d’occupation des sites) que nous étudions ici. 

Des lacunes se forment dans le réseau cristallin lorsque des atomes migrent de l’in¬ 
térieur du cristal vers sa surface. On notera e l’énergie de création d’une lacune et on 
considérera que les lacunes sont sans interactions au cours de leur formation, autrement 
dit, l’énergie (fixée) du cristal comprenant n lacunes vaut E = ne. 

i Quels sont les degrés de liberté qui permettent de caractériser les micro-états du 
cristal en présence de lacunes ? 

ii On se place désormais dans l’ensemble microcanonique, et on considère le cristal 
dans le macro-état à n lacunes, d’énergie E = ne fixée. 

Calculer la dégénérescence de ce macro-état, c’est-à-dire le nombre O de configu¬ 
rations microscopiques différentes (micro-états) correspondant à ce macro-état. 

iii En déduire l’expression de l’entropie. 

iv On introduit c, la concentration en lacunes, définie par c = -ft. 

En utilisant l’approximation de Stirling, montrer que l’entropie S du cristal 
s’écrit : 

S = — k-^N [c Inc + (1 — c) ln(l — c)] . 

v Utiliser le premier principe de la thermodynamique sous sa forme différentielle 
pour retrouver l’expression donnant la température microcanonique T à partir 
de l’entropie. 

vi Utiliser la définition de la température pour montrer que la concentration c 
vérifie : 

1 

où j3 est un paramètre que l’on précisera. 

vii On donne e = 1 eV. Calculer c à température ambiante. 

viii Calculer la température de fusion du cristal si l’on estime que 10% de lacunes 
suffisent à faire effondrer (et donc fondre) le cristal. 
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4- Exercices 


c Chaîne polymère et force entropique 



On décrit ici le rôle de l’entropie dans la force exercée sur une fibre de laine ou sur 
une chaîne polymérique (décrite sur la figure ci-dessus). Le brin est fixé à une extrémité 
et tracté par une force / à l’autre, de façon que sa longueur L soit maintenue constante. 
Les N liens entre monomères peuvent avoir deux longueurs possibles a ou b suivant leur 
état (a > b). Précisons que l’on ne s’intéresse ici qu’à l’entropie de configuration spatiale 
de la chaîne S(L, N), dont on suppose qu’elle ne dépend pas de l’énergie (un peu comme 
pour un gaz parfait). A priori, d’autres degrés de liberté contribuent à l’entropie avec 
une dépendance en énergie E mais pas en L, si bien que l’entropie totale peut s’écrire 
S tot (E, L, N) = S(L, N) + S'(E, N). On peut attribuer une température à la chaîne par 


l’equation — 


dS tot 

dE 


d S' 

1ËË' 


Ce calcul ne nous intéresse pas ici. 


i Exprimer la longueur de la chaîne L en fonction du nombre n a de liens dans 
l’état a, de N, a et b. On pourra également introduire les longueurs minimale et 
maximale, L m ; n et L max , que peut acquérir la chaîne. 


ii Déterminer le nombre de configurations accessibles pour L donné. En déduire 
l’expression de l’entropie S(L,N). 

iii En appliquant le premier principe à une transformation infinitésimale et réver¬ 
sible, montrer que 



iv Déterminer l’équation d’état f(L,T) puis L(f,T). En particulier, que vaut la 
longueur d’équilibre en l’absence de force, que l’on notera Lq ? Décrire les cas 
limites ainsi que la limite de petite force /. 

v On procède à l’expérience suivante : on pend un morceau de caoutchouc à un sup¬ 
port et on le tend en y accrochant une masse constante. On chauffe le caoutchouc 
à l’aide d’un sèche-cheveux. Que se passe-t-il? Est-ce surprenant? Interpréter 
physiquement le comportement. 
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Chapitre III. L’ensemble microcanonique 
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chapitre iv Ensemble canonique 


1 Distribution et fonction de partition canonique 

Une situation fréquente, dite canonique , est celle où un système constitué de N 
particules comprises dans un volume V peut échanger de l’énergie avec l’extérieur. 
C’est le cas, par exemple, lorsque le système est en relation avec un thermostat qui 
impose une température T. L’énergie du système n’est plus connue d’une façon absolue, 
mais fluctue autour d’une valeur moyenne (E) fixée (indirectement) par les conditions 
expér iment aies. 

Pour qualifier la situation canonique, on parle aussi de systèmes notés par (N, V, T) 
pour souligner quelles sont les variables non fluctuantes. On représente un ensemble à 
l’équilibre canonique par le schéma de la figure 4. 



Figure 4 - Système canonique à N, V fixés en équilibre avec un thermostat. 

a Expérience de deux systèmes en contact thermique 

Considérons un gaz confiné dans un volume, en contact avec un autre gaz, comme 
représenté sur le schéma 5. On suppose que le contact permet des échanges de chaleur 
mais pas de particules (de sorte qu’il est inutile de savoir si les gaz des deux systèmes 
sont différents ou non). L’ensemble est isolé. 

On notera E\ l’énergie du premier système, V\ son volume et N\ le nombre de 
particules qui le constituent, et £ 2 , V 2 , IV 2 les grandeurs statistiques analogues pour le 
second système. Les hypothèses physiques choisies au départ indiquent que V \, V 2 , N\ 
et N 2 sont conservées au cours du temps. 

De plus, E = E\ +E 2 , l’énergie totale de l’ensemble des systèmes, est conservée, car 
ce grand système est isolé ; on l’étudie dans l’ensemble microcanonique (V = V\ + V 2 et 
N = N 1 +N 2 sont également conservées, mais ce n’est pas pertinent et d’ailleurs N peut 
ne pas être une grandeur statistique, si les systèmes n’ont pas les mêmes constituants). 
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Figure 5 - Les systèmes 1 et 2 sont en contact le long de la paroi, qui ne laisse pas 
passer les particules mais permet les échanges thermiques. 


ol Entropie maximale à l’équilibre 

Comme le grand système est dans l’ensemble microcanonique, à l’équilibre son en¬ 
tropie doit être maximale. Comme les systèmes sont rigidement séparés, le nombre de 
micro-états s’écrit 1 fl = et S = k-Q log(fi) = S\ + S 2l où Si est l’entropie du 

système i. 

S doit être maximale par rapport à toute variable, en particulier par rapport à E\ 
(ou E‘ 2 - ça donnerait le même résultat). On écrit : 



à l'équilibre 


dSi dS 2 
dËh + ~cLËh 


dSi 

dËl 


+ 


Vl.iVl 

ds 2 


dS 2 dE — E\ 
dE 2 dE\ 


= -1 


1 

T[ ~ 3E 2 

1 1 

t v ~T 2 ' 


v 2 ,n 2 


(15) 


où on a utilisé le fait que E, V\, V 2 , Ni et N 2 sont conservés. T\ et T 2 sont les tempéra¬ 
tures thermodynamiques des deux systèmes selon la définition (13a). On en déduit qu’à 
l’équilibre T\ = T 2 : les deux systèmes ont la même température, dès lors qu’on autorise 
les échanges de chaleur. 

(3 Bain thermostatique 

Lorsqu’un des deux systèmes est plus important que l’autre (soit N 2 2> Ah, soit 
V 2 3> Vi, ces deux conditions sont équivalentes à E 2 2> E±), le couplage entre les deux 
ne peut perturber le plus grand. Sa température est fixée et le résultat qu’on vient de 
montrer prouve que le petit système adopte la température du grand. 

On dit que le petit système est thermostaté, ou que le grand est un bain thermosta¬ 
tique. D’autres configurations existent pour imposer la température d’un système. 


1. Il n’y a pas de facteur d’indiscernabilité entre les deux systèmes, puisqu’on peut toujours déter¬ 
miner dans quel système une particule se trouve. Les facteurs d’indiscernabilité sont l/Aù! et 1/N 2 \. 
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2. Entropie d ’équilibre 


b Application du principe d’entropie maximale 

L’application du principe d’entropie maximale (cf. II 4) conduit à l’équation 
d 


0P„ 


^2 P m l n Pm + Aq ( ^2 P m ~ 1 J + A e I ^2 Pm.Prn — (E) 


= 0 


avec des notations comprénehsibles. On obtient ainsi 

n -l+±Q-+^-E(m) 

P m = e fc B *b 

et, en utilisant encore la contrainte P m = 1, on obtient la distribution canonique 
(ou distribution de Gibbs-Boltzmann), 

e fc B 

P — 

1 m. — 


— _B(m) 


Z(X E ) 


où la constante de normalisation 


z(\ E ) = j2 cl 


( E(m) 


est connue sous le nom de fonction de partition canonique. Cette distribution n’est pas 
uniforme, mais pondérée exponentiellement par l’énergie du micro-état concerné. 

Le multiplicateur de Lagrange X E est, quant à lui, fixé par la deuxième contrainte: 


U = ^2 p(m)E(m) , 


où on a utilisé la notation traditionnelle U = (E) é pour désigner l’énergie interne 
des états d’équilibre en physique macroscopique. On montre que X E s’obtient de façon 
implicite par l’expression : 

d ln Z(X E ) TT 

kB ~â>^ = u - 

La distribution canonique se trouve ainsi parfaitement caractérisée. La signification 
physique de À^; est expliquée à la section suivante. 

2 Entropie d’équilibre 

a Statistique de Boltzmann 

Calculons l’entropie. Puisque p(m) = Z -1 e fc B un ca ^ cu j élémentaire donne 

5 = k B lnZ(X E ) -X E U . 

Mais comme U = k B dln Z/dX E , d’après la transformation de Legendre (cf. an¬ 
nexe E), il apparaît que S (U) est la transformée de Legendre de k B \nZ{X E ). Autrement 
dit, 2 dS = —X E dU et on en déduit donc que 

dS 


X E = — 


dU 


2. On a: d(k B lnZ) = k ^ z d\ E = Ud\ E , d’où dS = Ud\ E - X E dU - UdX E = -X E dU. 
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Chapitre IV. Ensemble canonique 


Cette formule établit un lien entre une grandeur issue de notre calcul statistique À# et la 
fonction d’état S (U) de la thermodynamique macroscopique. Le premier principe de la 
thermodynamique (lia) a établi que dS/dU = 1/T. Un usage très répandu est de noter 
(3 = l/(kp,T) proportionnel à l’inverse de la température. Les expressions définitives de 
la distribution canonique (ou distribution de Boltzmann ) et de la fonction de partition 
canonique sont finalement données par : 


e -/3E(m ) 


Pim,= Z(l3) ’ 

Z (P) = e^ %) 

m 


À température fixée, les états microscopiques les plus probables (expo¬ 
nentiellement) sont les états de basse énergie. 

Dans le cadre classique, les formules précédemment définies (fonctions de partition, 
moyennes, ..) devront être interprétées avec la règle (1). La fonction de partition cano¬ 
nique des systèmes constitués de particules classiques indiscernables 3 s’écrit : 

1 N) - 1 [ d 3N rd 3N p 

| ' > N\ J h 3N _' 

b Distribution de Maxwell 



On peut aller un peu plus loin car les degrés de liberté d’impulsion peuvent être 
complètement intégrés dans le cas des systèmes classiques. En effet, l’énergie totale 
s’écrit : 

E(? K ,0 N ) = . 

où Un représente l’énergie potentielle d’interaction entre les particules de masse m et/ou 
avec un champ extérieur. La fonction de partition peut donc s’écrire : 


Z(f3,N) 


1 

m 



d 3N p 

JW 


,-/3£Üi 


lia il 2 

2m 


d 3N r e~PU N (r ) 


Uv 


(16) 


où V désigne volume physique où les particules sont libres de se mouvoir. 

► L’intégrale multiple sur les impulsions peut être calculée de façon exacte. Cela résulte 
des propriétés de factorisation de l’exponentielle: e -/3 £dlftll / 2m = e~P\\Pi\\ / 2m ; e t 

des intégrales gaussiennes (cf. annexe B) : 





/ r fp /2îrm\ 2 _ 

• V J i fih 2 J - 


où on a introduit la longueur d’onde thermique (ou longueur thermique de de 
Broglie) d’une particule de masse m portée à la température T par 


A = J3h 2 /2irm) 112 


h 

y/2îrmk^T 


3. La division par N\ doit être omise lorsqu’on traite des particules discernables. 


(17) 
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3. Potentiel thermodynamique 


► On peut finalement écrire la fonction de partition canonique des systèmes constitués 
de particules indiscernables obéissant à la mécanique classique sous la forme 4 : 

WW = ^3 nQnW,N)\ (18a) 

Qn{P,N) = J d 3 r N e -P u (.r N ) . ( 18b ) 

Qn s’appelle l’intégrale de configuration. 

► On notera en particulier que, pour des particules classiques à l’équilibre thermo¬ 
dynamique, la loi de distribution des impulsions, V(p N ), s’écrit, indépendamment des 
positions spatiales, sous la forme : 


N 

v(pi, ■ ■ ■ ,p N ) = n (2vr mk B T)- 3/2 . (19) 

i= 1 

Il s’agit de la loi de distribution des impulsions de Maxwell, 
c Cas de particules indépendantes 

Pour des particules ou degrés de liberté discernables tels que E( = £ i{ n i)-, alors 

Z = avec Zi = ^^e _,â£ ^ n ^ (cas discret) ou Zi = f ^ ^ e -3 h idtv) ( cas continu). 

in J 


Si, de plus, les particules ou degrés de liberté sont identiques, on a 


Z = z 


N 


si les particules sont identiques et indiscernables, 


Z = 


N\ 


(20a) 


(20b) 


3 Potentiel thermodynamique 

Pour le moment, la fonction de partition Z intervient comme une constante de 
normalisation à laquelle n’est attachée aucune signification physique. Comme on va le 
découvrir maintenant, il apparaît que ce n’est pas Z mais son logarithme qui a une 
signification physique simple. 

Pour s’en convaincre, substituons l’expression \e = —1/T dans l’entropie d’équilibre 
canonique. On trouve : 

F = — k^T In Z , 

4. Si les particules sont discernables, il faut omettre le facteur N\. 
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où on a introduit la notation F = U — T S qui représente, en physique macroscopique, 
l’énergie libre 5 du système étudié. Ainsi, dans l’ensembe canonique, le logarithme de 
la fonction de partition s’identifie avec l’énergie libre du système. Il est important de 
souligner que Z étant une fonction de /?, c’est-à-dire de la température T, on trouve 
que l’énergie libre F = F (T) est une fonction naturelle de T. Ceci découle du fait que 
F (T) est la transformée de Legendre 6 de U (S). 

Pour rendre les choses un peu plus concrètes, considérons le cas d’un fluide (gaz ou 
liquide), contenu dans un récipient de volume V, comportant N molécules, et en équi¬ 
libre avec un thermostat qui lui impose la température T. Dans cette situation simple, 
au cours d’une transformation où le volume varie de dV et le nombre de particules de 
dN, le fluide reçoit le travail ôW = —p dV + p dN, où p et p représentent la pression et 
le potentiel chimique. La différentielles dF peut être obtenue à partir de (12) et s’écrire : 

dU = TdS — pdV + pdN => dF = —S dT — pdV + pdN . (21) 

La relation (21) est très importante, car elle permet de calculer directement l’entropie, 
la pression et le potentiel chimique à partir de l’expression statistique de l’énergie libre : 


r , dF dF dF 

= ~~dT ’ P= ~~dV ’ P=+ ~dN' 


C’est pour cette raison qu’on dira que le potentiel thermodynamique dans une situation 
d’équilibre canonique est donné par l’énergie libre F du système étudié. Pour compléter 
le tableau, on se souviendra que l’énergie interne 7 a été obtenue précédemment, elle 
aussi, par dérivation de ln Z : 


dlnZ _ d(PF) 

Ô/T ~ dp ’ 

de sorte que l’obtention des grandeurs physiques au moyen de la physique statistique, 
dans le cadre canonique, résulte uniquement du calcul de l’énergie libre F = —k-gT ln Z 
(qui prend le rôle de potentiel thermodynamique). 

4 Théorème d’équipartition de l’énergie 

Les système pour lesquels on sait calculer les fonctions de partition sans approxima¬ 
tion sont évidemment exceptionnels. Presque tous les cas connus correspondent à des 
situations unidimensionnelles où les calculs, quoique complexes, peuvent parfois être 
menés à leur terme 8 . Il existe cependant un cas, celui des systèmes classiques avec in¬ 
teractions quadratiques, que nous discutons dans cette section, pour lequel un calcul 
exact est possible. Ce résultat particulier constitue le théorème dit d’équipartition de 
l ’énergie. 

5. Précisément, l’énergie libre de Helmholtz. 

6. En effet, comme T = dU/dS (par le premier principe), on a F(T) = U(S) — S 

7. On peut utiliser directement U = F + T S, une fois F et S calculées. 

8. Un contre-exemple fameux et riche d’enseignements, est le calcul exact par Onsager de l’énergie 
libre du modèle d’Ising ferromagnétique à deux dimensions d’espace. 
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5. Calcul de la pression en théorie cinétique 


Le théorème d’équipartition de l’énergie repose sur les égalités mathématiques sui¬ 
vantes : 


z = l e /3æ2 / 2 dx 

J K 



(e) = 


clin z 

dp 


1 

20 


k B T 

2 


( 22 ) 


► Dans ces expressions, z représente la fonction de partition canonique d’un système 
à un seul degré de liberté dont l’énergie s’écrit e(x) = x 2 /2 ; (e) s’identifie alors avec 
l’énergie moyenne correspondante. 

► La formule (22) se généralise au cas d’un système où l’énergie totale est une fonction 
quadratique, définie positive, de n coordonnées généralisées qu’on note génériquement q^. 
Un cas correspondant à cette situation est celui des oscillateurs linéaires en interaction. 
En utilisant les résultats donnés en annexe F et pour les mêmes raisons mathématiques 
que ci-dessus, on trouve pour ce système : 


W = è£L=i^ M ^- 


(E) 


d 1 n Z k^T 

< 9/3 ~ n 2 i 


où M est une matrice n X n réelle et symétrique. 


Pour un système classique dans l’ensemble canonique, la contribution d’un terme 
quadratique du hamiltonien à l’énergie moyenne est k B T /2. 


► Ce résultat constitue le théorème d ’équipartition de l ’énergie , valable pour tous les 
systèmes classiques, dont l’énergie est une fonction quadratique des variables de l’espace 
des phases. 

Ainsi dans le cas du gaz parfait, l’énergie est bien une fonction quadratique des 
impulsions qui constituent les seuls degrés de libertés. Pour N atomes dans le gaz, 
on aura donc U = (3N)k B T /2. Dans le cas d’un seul oscillateur harmonique à une 
dimension, il y a 2 degrés de libertés : celui associé à l’énergie potentielle kx 2 /2 , et celui 
associée à l’énergie cinétique p 2 /(2m). Pour N oscillateurs indépendants, on aura donc 
U = N k B T. Ces calculs élémentaires permettent de retrouver rapidement les différences 
entre les chaleurs spécifiques des gaz atomiques et moléculaires à haute température. 

5 Calcul de la pression en théorie cinétique 

La pression thermodynamique est définie par (13b) mais la pression a été définie 
localement par Boltzmann, on l’appelle aussi pression cinétique. C’est la valeur moyenne 
du flux sortant d’impulsion des particules à travers une surface par unité de surface et 
de temps. 9 

On calcule la pression cinétique p xx pour une surface plane S yz , orthogonale à l’axe 
x'Ox, d’aire dS : 

9. Attention, il s’agit d’un flux partiel, on sépare les particules sortantes des particules entrantes 
(les deux flux correspondants contribuent à des pressions opposées sur la surface et s’équilibrent si l’on 
ne tient pas compte des fluctuations). 
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• On calcule, en premier lieu, la pression partielle créée par les particules ayant une 
vitesse p/m déterminée à dp/m près. Pour atteindre la surface S x et intervenir 
dans le calcul de la pression, ces particules doivent être situé, au temps t, dans 
le volume dV, un cylindre d’arête vdt = p dt/m et de base dS, m étant la masse 
des particules. On vérifie sur la figure suivante que ce volume correspond aux 
particules de vitesse v = p/m. sortant à travers S x . 


Il s’écrit exactement 



m m 

et le nombre de particules qui traversent la surface S x avec une impulsion com¬ 
prise entre p et p + dp pendant l’intervalle de temps dt s’écrit dN = ndV, où 
n = N/V est la densité volumique. 

• La pression 10 est une force par unité de surface, soit 


Pxx ~ dS 


Pour ressentir la force f x , faisons l’expérience de pensée suivante: on matéria¬ 
lise la surface élémentaire S x . D’après le principe d’action et réaction, * 11 f x est 
l’opposée de la force F x ressentie par les particules. Or, d’après le théorème fon¬ 
damental de la dynamique, F x est égale à la variation d’impulsion totale par 
unité de temps et s’écrit 


F x 


E 

particules sortant 
a travers S x 


A p x 
dt 


où A p x = p' x — p x est la variation de quantité de mouvement en rebondissant sur 
la surface. Par définition, p x est orientée vers l’extérieur et p' x vers l’intérieur. 

10. Attention, ne pas confondre p x la composante sur x de l’impulsion d’une particule et p xx la 
composante de la pression. 

11. Il s’applique ici car toutes les forces sont conservatives et dérivent d’un potentiel. On rappelle 
que les contre-exemples sont les forces de frottement (non conservatives) et la force de Lorentz (créée 
par un champ magnétique) qui n’a pas de potentiel. 
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6. Stabilité thermodynamique 


• En moyenne, (p x ) = — (p' x ) car on peut supposer, dans cette expérience de 
pensée, que les collisions sont élastiques. Alors, (A p x ) = —2(p x ). En réalité, il 
n’y a pas de surface matérielle, donc pas de collisions, mais on étend ce calcul 
par analogie et on calcule F x en posant A p x = —2p x . 

• Finalement, f x s’écrit (le deux signes — se compensent) : f x = dN2p x /dt et la 
pression partielle vaut 

2dN p x _ 2 np x 
dSdt m. 

• Pour étendre le calcul à toutes les impulsions et obtenir la pression totale, il faut 
intégrer sur toutes les impulsions en pondérant par la densité de probabilité V à 
une particule dans l’espace des impulsions. 


Pxx — 




dp Z V(p x ,Py,P z ) 


2 np x 
ru. 


n 

m 


W) 


où le facteur 1/2 supplémentaire provient de la substitution j 0 °° dp x en ////. dp x 
quand on introduit la moyenne ( ). 

• Comme on considère un gaz isotrope, on a {p x ) = ( p y 2 ) = ( p z 2 ), d’où 

{Px) = l (p 2 ) • 


et 



(23a) 


• Comme le système est isotrope, la pression s’écrit pi oc = p xx . Pour le cas par¬ 
ticulier de l’ensemble canonique, on choisit logiquement la loi de Maxwell (19) 
V{\Jp x + Py + p z 2 ). La pression totale devient 


2 n 


Pxx — 


m,(2irmkBT) 3 / 2 J 0 
nmk-QT Nk^T 


dp x e 


—rEz_ 

' 2m 


dp v e 


-p p -l 

r 2m. 


dp z e 


-Pïk 


m 


V 


où on intègre seulement sur les p x > 0 (mais on a bien alors le facteur 2n/ru. 
initial). Dans ce cas particulier, on retrouve la pression thermodynamique. 


6 Stabilité thermodynamique 


a Minimisation de l’énergie libre 

L’énergie libre F(T,V,N) est minimale à l’équilibre thermodynamique. Nous ne le 
démontrerons pas. Si on perturbe le système, à T, V, N fixés, une conséquence très 
importante est que F reste constante jusqu’au retour à l’équilibre. 

Ce principe peut être utilisé uniquement quand on perturbe légèrement le système, 
par contre, si on le perturbe trop fortement, il peut choisir un nouvel état d’équilibre. 
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b Fluctuations en énergie 


On reprend le petit modèle expliqué en introduction : le système canonique corres¬ 
pond au gaz numéroté 1, et le bain thermostatique au 2. On va étudier les fluctuations 
de E\. Pour cela, on va calculer les contributions, parmi les fl micro-états du système 
total 1 + 2, de ceux d’énergie E\. 

a Développement de l’entropie 

On calcule d’abord le développement autour de l’équilibre de l’entropie totale S = 
Si + S 2 , en fonction de E\. On notera U± sa valeur d’équilibre. 

À l’équilibre, = 0, ce qui prouve que le développement n’a pas de terme du 
premier ordre. 

Pour calculer le terme du second ordre, {Ex — Ci) 2 , on part de l’expression 

générale hors équilibre de donnée en (15), que l’on dérive: 


d 2 S 
dE 2 


1 dE 1 1 dE ‘2 


T 2 dE 1 
1 dT\ 
T 2 ÜËi 


+ 


T 2 dE 1 


+ 


1 dE-Ex dT 2 


T 2 

1 2 


dE 1 dE 2 


T 2 C V 1 


T 2 C V 2 


On doit estimer cette quantité à l’équilibre thermodynamique, quand T\ = T 2 . De 
nombreux auteurs concluent, en particularisant le calcul au cas d’un gaz parfait, à 
l’aide de (14). On va plutôt profiter de l’extensivité de la capacité calorifique et définir 
la chaleur spécifique moléculaire : 


c v — C v /N , 


qui est intensive. D’où 

d 2 S _ 1 Nx c vl _ 1 

dE 2 c vl T 2 Nx N 2 c v2 ’~ c v1 T 2 Nx 

car le facteur Nx/N 2 <C 1 par définition d’un bain thermostatique, tandis que le rapport 
des chaleurs spécifiques est un nombre donné, qui ne peut être ni trop grand, ni trop 
petit. 

En définitive, on a, hors de l’équilibre, le développement au second ordre : 

S(Ex) = S(Ux) - _ * (Ex - Ux) 2 . (24) 

2c vl T 1 z N 1 

f3 Calcul de la densité dCL/dEx 

Comme l’ensemble 1 + 2 est microcanonique, on peut utiliser directement la formule 
(8) à l’envers, d’où, en substituant (24), 

SjUj) ( E1-U1 ) 2 (Ei-U !) 2 

D = e + e 2c vi k B T i 2N i = Ç} eq e 2 <t i n i ? (25) 
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6. Stabilité thermodynamique 


ggi) 

où on a posé o\ = T\\/c v ik-Q, qui est une grandeur intensive et fi eq = e fc B . 

Or, on n’a ainsi compté que les micro-états correspondant à une énergie E\ donnée, 
la bonne expression de O est, toujours dans cette approximation du second ordre, 


n = 


n 


tôt 


(Jw/2'kN] 


i Is. 


(O-O ) 2 

2 ai 2 Ni 


dE\ 


= n 


tôt 


où on a posé f2 tot = n eq (7i \/2tiN\ . 

On voit que (25) donne la densité en fonction de E\ recherchée. 

7 Calcul des fluctuations 

On peut oublier l’indice 1 dans les calculs précédents et considérer le système dans 
l’ensemble canonique seul. La répartition des micro-états en fonction de l’énergie totale 
E est une gaussienne centrée autour de U, comme représenté sur la courbe ci-dessus : 



Comme on le constate, O le nombre de micro-états qui ne sont pas dans l’intervalle 
[U — ÔE , U + SE] est négligeable devant O. On a Û <C O, quand SE = a\fN. avec a = 
T^~v l’écart type de la gaussienne, qui donne la grandeur typique des fluctuations 
de l’énergie. 

On est exactement dans la situation prédite par le théorème de la limite centrale 
(voir annexe A) : les fluctuations sont importantes en valeur absolue, puisqu’elles sont 
de l’ordre de y/N, mais elles sont négligeables en valeur relative. Si on examine, par 
exemple, les fluctuations de e = E/N, l’énergie par particule, on trouve 



qui est nul à la limite thermodynamique, 12 ce qu’on peut récrire : 


A E 1 

TE\ ** ~r= ^ 1 • 

(E) 


12. En pratique, l/y/~N <10 10 . 
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Chapitre IV. Ensemble canonique 

c Équivalence des ensembles 

L’analyse des fluctuations de l’énergie, qui viennent d’être calculées, permet de jus¬ 
tifier le principe d’équivalence des ensembles, qui a été posé précédemment. En résumé, 
on peut écrire 

À la limite thermodynamique N 3> 1, les fluctuations de l’énergie sont négli¬ 
geables et les fonctions thermodynamiques calculées dans les ensembles canonique 
et microcanonique sont identiques. 

En conséquence, on peut choisir l’ensemble avec lequel on travaille si le système est 
assez grand. 
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7 Exercices 

a Système d’oscillateurs harmoniques quantiques unidimensionnels 

On rappelle que les niveaux d’énergie d’un oscillateur harmonique à une dimension, 
calculés de façon quantique, s’écrivent e n = hüo(n + 1/2), où n est un entier positif ou 
nul, et sont non-dégénérés. 

i Calculer la fonction de partition canonique d’un seul oscillateur. 

ii En déduire la fonction de partition canonique de N oscillateurs indépendants et 
discernables (par leur position par exemple). 

iii Calculer l’énergie interne U (T) du système d’oscillateurs et la représenter. 

iv Ce modèle d’oscillateurs indépendants a été proposé par Einstein comme un 
modèle de solide permettant de calculer la chaleur spécifique C des isolants. 
L’expérience montre qu’à haute température, C = 3Nk b (loi de Dulong et Petit), 
tandis qu’à basse température C ~ T~ 3 . 

Le modèle d’Einstein est-il conforme à l’expérience? 

b Dipôles classiques dans dans un champ extérieur 

On considère un ensemble de N dipôles électriques discernables et identiques, notés 
fh et soumis à un champ électrique externe £ constant. Le système est supposé suffi¬ 
samment dilué pour qu’on puisse négliger les interactions entre dipôles et on ne tiendra 
pas compte des degrés de liberté de translation (parce que, par exemple, les dipôles sont 
situés aux nœuds d’un réseau). 

On résout ce système classiquement. L’énergie potentielle d’interaction entre un 
dipôle et le champ électrique s’écrit E ( j = —d-£. On considérera que le champ électrique 
est selon la direction Oz et on note d = ||d|| et £ = ||£||. On repère la direction 
d’un dipôle par les variables angulaires (9,(p) des coordonnées sphériques, 9 G [0,7r] et 
ip G [0, 2tt]. Le système est en contact avec un thermostat qui fixe la température T. 

Raisonnement avec la loi de Boltzmann 

On étudie l’orientation moyenne des dipôles. On la calcule pour un dipôle seul. 

i Quel est l’effet du champ électrique sur un dipôle ? Exprimer l’énergie potentielle 
d’un dipôle en fonction de d, £ et 6. 

ii Lorsque £ = 0, quelle est la probabilité d’avoir les angles ( 9 , ip) ? 

On rappelle que l’angle solide élémentaire vaut d£l = sin 9d9dp>. 

iii Utiliser la loi de Boltzmann pour déterminer la densité de probabilité p{9 ) d’avoir 
l’angle 9. 

iv En déduire (d) comme fonction de d £/k^T et représenter (d z ) sur une courbe. 
On notera C(x) = coth(x) — 1/x la fonction de Langevin. 
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Raisonnement dans l’espace des phases 

Les impulsions angulaires associées aux variables angulaires 13 (0, <p) sont notées 
( P 0 iP<p ) e t telles que l’hamiltonien d’un dipôle prend la forme 

h(0,<p,pe,p<p) = 777 ( Pe 2 + -^ 27 ; ) - d£ cos 6 , 

21 \ sin 0 J 

où I est le moment d’inertie du dipôle. Le premier terme correspond à l’énergie cinétique 
de rotation. On adoptera une description dans l’espace des phases dans laquelle le 
volume élémentaire contenant un état d’un seul dipôle est h 2 ; c’est-à-dire que l’on 
divisera par h 2 chaque volume d’intégration dOdpdpgdp ip . 

v Soit z la fonction de partition d’un seul dipôle. Exprimer la fonction de partition 
totale Z en fonction de z, puis calculer z. 

vi En déduire l’énergie libre par dipôle /. 

vii Montrer que le dipôle moyen vaut ( d ) = — ^ e 2 . 

viii Retrouver alors l’expression de (d) trouvée auparavant. 

c Système à deux niveaux dans l’ensemble canonique 

On considère N particules discernables, chacune peut se trouver dans deux niveaux 
possibles d’énergie respective 0 ou e. Le système est maintenu à une température T. 

i Comment décrit-on un micro-état ? que vaut l’énergie totale ? 

ii Ecrire la fonction de partition Z de l’ensemble des particules et montrer que l’on 
peut l’écrire Z = z N avec z la fonction de partition à une particule. Calculer z. 

iii En déduire l’énergie libre, l’énergie moyenne et l’entropie. 

iv Quelles sont les fractions /o et f £ de particules dans l’état 0 et dans l’état e? 
Montrer que l’énergie moyenne s’écrit simplement en fonction de ces fractions. 
Interpréter physiquement les différents cas limites. 

v Exprimer l’entropie S en fonction de ces fractions. Ce résultat vous rappelle-t-il 
un exercice du chapitre III ? 

vi Définir et calculer la capacité calorique C v . 

vii Calculer les fluctuations en énergie du système. Sont-elles reliées à C V 1 
Que peut-on dire de ces fluctuations à la limite thermodynamique ? 


13. Pour information, on a pe = 16 et p v = I sin 2 Otp. Ne pas utiliser ces formules dans l’exercice. 
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La distribution canonique correspondait aux systèmes pour lesquels les variables 
(N, V, T) étaient fixées, l’énergie libre F étant le potentiel thermodynamique associé. 
Dans ce chapitre, on traite un autre cas important : l’ensemble grand-canonique, où 
les variables fixées sont (//, V,T). Le système étudié est donc en équilibre à la fois 
avec un réservoir d’énergie (présence d’un thermostat, température T fixée) et avec un 
réservoir de particules (potentiel chimique qi fixé). Comme exemple, on peut citer le cas 
de la surface libre d’un liquide en contact avec sa vapeur où un échange permanent de 
molécules est possible par l’interface liquide-gaz. 

1 Entropie d’équilibre et grande fonction de partition 

Dans ce système, l’énergie et le nombre de particules peuvent fluctuer indépendam¬ 
ment. On le représente schématiquement à la figure 6. 



FIGURE 6 - Système à volume V fixé en équilibre avec un thermostat et un réservoir 
de particules. 

Afin de déterminer la distribution d’équilibre grand-canonique, on doit maximiser 
l’entropie statistique sous les trois contraintes : 

^p(m) = 1 , 

m 

\p{m)E(m ) = (E) = cte , 

m 

^ ~^p(m)N{m ) = (N) = cte . 

m 

En introduisant 3 multiplicateurs de Lagrange Àq, et A^r (puisqu’on a 3 contraintes), 
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on arrive au résultat : 


p(m) 


g(A_B-E+AArW)/fc B 

Z(lEV,T) 


où H représente le coefficient de normalisation (pour assurer la condition de normalisa¬ 
tion 5 ~2 m P ( m ) = 1) ! S est la grande fonction de partition: 


a= > e 


E' 


,(\ E E+\ N N)/k B 


Il reste à identifier les multiplicateurs de Lagrange A e et Aat, ce qui se fait, comme 
dans l’ensemble canonique, en passant par l’entropie physique S. On trouve : 

S = k B lnE-X E (E) — \ N (N) . 

Or, un calcul direct montre que : 


(E) = k B dlnE/dX E , et (N) = k B d ln E/<9A/y . 


S((E ), ( N))/k B est donc la (double) transformée de Legendre de ln E. En conséquence, 
on peut écrire 


S 


, , _ TT ds IKT . ds 

k B ln* + U — + (N) — , 


où on a encore utilisé la notation 1 (E) = U. En utilisant le premier principe, 1 2 on 
exprime les dérivées partielles de l’entropie pour obtenir : 


—A e 


—Xn 


dS_ _ 1_ 
dû ~T ’ 
dS _ n 
ÔN ~~T ■ 


Dans l’ensemble canonique, on a introduit /3 = —Xe/ k B . De plus, il découle des 
relations précédentes que A^t/Zcb = +/3/U. Finalement, la grande fonction de partition 
est donnée par : 




e -f}(E(m)-p,N(m)) 


m 


[Cif Remarques : 

• Comme le nombre de particules N n’est plus fixé, la somme sur m doit comp¬ 
tabiliser tous les micro-états, quel que soit le nombre de particules du système 
(dans l’ensemble canonique, la somme sur les micro-états s’identifiait avec une 
somme sur les états d’énergie arbitraire à nombre de particules N fixé). 

1. De la même façon que l’on a distingué la variable dynamique E et sa valeur moyenne (E), on 
devrait distinguer N et sa valeur moyenne (N). Les vraies variables conjuguées sont ((N) ,/x) et non 
(N, y), mais cet abus de notation est quasiment universel. 

2. dS = ^ + f dV - £ dN. 
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• Ainsi, pour les systèmes classiques, on doit rajouter une somme sur N dans 
la description des états microscopiques 

E(jU, V,T) = yry J — J — é~^ E ^ g ’ p ^~ pN ^ (systèmes classiques). 

On notera également le lien entre les fonctions de partition canonique et grand- 
canonique (valable pour les systèmes classiques et quantiques) : 


E(/3, M ) = ^e^ZO^TV) . 

N =0 

Cette dernière relation montre que la grande fonction de partition est la trans¬ 
formée de Laplace (discrète) de la fonction de partition canonique. 

2 Grand potentiel thermodynamique 


A la différence de la situation dans l’ensemble canonique, la variable N est une 
variable fluctuante dans l’ensemble grand-canonique. Par conséquent, en termes macro¬ 
scopiques, la bonne fonction d’état n’est plus l’énergie libre F(N,V,T ) mais une autre 
fonction J(/i,V,T). On peut passer de F(N,V,T ) à J(n,V,T ) par une transformation 
de Legendre sur le couple (N, p) : 

J(/x, V, T) = F(N , V ,T)-N^- = F-Np = U-TS-pN 


et on s’attend à ce que la fonction d’état J(/x, V, T) joue le rôle de potentiel thermo¬ 
dynamique dans l’ensemble grand-canonique, comme F(N,V,T ) l’est dans l’ensemble 
canonique. Cette fonction prend le nom de grand-potentiel. 

Vérifions maintenant, à partir du formalisme thermodynamique, que </(//, V, T) joue 
bien ce rôle de potentiel thermodynamique. Dans la section précédente, on a montré 
que 

S = k B lnZ + ^-^(N) . 

On en déduit que : 

( 26 ) 


J = U — TS — p (N) = -k B T lnS 


Pour déterminer les grandeurs physiques, il suffit d’exprimer la différentielle de 
J (fi, V, T). On trouve ainsi : 


dJ = — (N) dp-pdV - SdT 


d’où l’on déduit 


/Ar\ d J 

d J 

dJ II 

{N) - -ëï ■ 

P = ~dV ’ 

Ce 

II 

Si 


Il reste à déterminer l’énergie (interne) U. Puisque (N) et S peuvent être obtenues 
par les formules précédentes, le plus simple est d’utiliser la relation : 


U = J + TS + pL (N) . 
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3 Le gaz parfait dans le cadre grand-canonique 


Essayons à présent, en utilisant le formalisme grand-canonique, de retrouver les 
résultats du chapitre précédent sur l’expression du potentiel chimique /3p gp = InnA 3 , 
où A est donné par (17) et l’équation d’état des gaz parfaits. 

► Le nombre total N de particules n’est plus maintenant fixé de façon absolue, mais 
peut fluctuer autour d’une valeur moyenne (N) grâce au paramètre de contrôle que 
constitue le potentiel chimique p. La grande fonction de partition s’écrit 

~ 1 fe^V\ N 

N\ V A 3 ) 


Z(T,V,ii)=J2 ePtlN Z(P,N,V)=J2 


v, 
= e * 1 


.Pu 


N=o 


N=0 


soit encore, en utilisant la relation (26) et (27b), 


e^y = s(Tj y- ^ = J» ^ Pp{T, V,p) = ■ 

► Sous cette forme, on ne reconnaît guère la loi des gaz parfaits. La raison tient dans 
le fait que cette équation d’état n’est pas écrite en fonction des variables usuelles, N, V 
et T, mais en fonction de p, V ., T. Pour la comparer avec l’équation d’état sous sa forme 
habituelle, il faut donc revenir à la variable N. Dans notre cas, cela s’obtient en calculant 
le nombre (moyen) (N) de particules : 

(N) = = -^ 4» n = Pp . 

dp A 6 

On retrouve ainsi à la fois l’expression du potentiel chimique fip gp = InnA 3 et 
l’équation d’état des gaz parfaits. 

► Comme on le verra par la suite, dans des situations plus complexes, cette procédure 
est très générale quand on se place dans l’ensemble grand-canonique. Le calcul (éven¬ 
tuellement approché) de InS permet d’obtenir directement la pression en fonction du 
potentiel chimique. Il faut ensuite procéder à une inversion de la relation liant N à p 
pour retrouver la pression en fonction de la variable naturelle qu’est la densité. 


4 Étude d’un gaz diatomique 


Les molécules diatomiques se comportent typiquement comme de petits haltères liés 
par un ressort, modèle de liaison chimique simpliste mais qui nous sera suffisant. Par 
rapport au cas d’un gaz mono-atomique, cela introduit deux termes supplémentaires 
dans l’énergie : une énergie associée à la vibration de la liaison et une associée à la 
rotation de la molécule. 

a Hamiltonien et degrés de libertés 

On peut décomposer le hamiltonien selon les différents degrés de liberté qui peuvent 
stocker de l’énergie. On sépare les termes à une particule et ceux à pluieurs particules. 
Ces derniers décrivent les interactions entre molécules, via un potentiel d’interaction à 
deux corps, mais on négligera ces termes par la suite, et on se focalisera sur les premiers 
termes du hamiltonien : 
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4- Etude d’un gaz diatomique 


• énergie cinétique de translation ; 

• énergie de rotation propre ; 

• énergie de vibration de la liaison ; 

• énergie potentielle dans un champ extérieur. 

Nous allons étudier la contribution de ces différents termes à la thermodynamique 
des gaz de molécules diatomiques. 

b Capacité calorifique d’un gaz de molécules diatomiques 

Le but de cette partie est de dériver le comportement de la capacité calorifique 
C'y (T) des gaz diatomiques par les méthodes de physique statistique. On a vu dans le 
chapitre introductif que, bien que la fonction C'y (T) soit une donnée essentielle dans 
les bilans thermodynamiques, son comportement ne peut être compris dans le cadre 
de la thermodynamique. Nous allons voir que les deux contributions de rotation et de 
vibration permettent d’expliquer la dépendance en température de C'y (T) pour peu que 
l’on prenne en compte le spectre quantique. 

Les différents degrés de liberté d’une molécule donnent des contributions séparées, 
qui se factorisent généralement dans la fonction de partition 2 . Les degrés de liberté de 
rotation peuvent être traités de façon classique, comme à l’exercice 8 c. Les degrés de 
liberté de vibration sont bien décrits par ceux d’un oscillateur harmonique unidimen¬ 
sionnel (selon l’axe de la molécule), comme traitée à la suite. 

a Approximation harmonique du terme de vibration 

La liaison chimique se modélise raisonnablement bien par un potentiel d’interaction 
v{p) avec p la distance relative entre les deux atomes. Ce potentiel présente un minimum 
pour une distance po qui est la longueur de la liaison à l’équilibre. Les petits mouvements 
autour de ce minimum sont bien décrits par un oscillateur harmonique de pulsation 
oj = \M(pÔVm avec A* la masse réduite et on peut écrire v(p) ~ v(po) + " ^ (p — po) 2 - 
(3 Traitement classique avec un bain thermostatique 

Dans ces conditions, chaque molécule est indépendante et (20b) s’applique, z étant 
déterminé par 


1 r°° r°° p 2 

^vib = t dx e 2 A: b t dp e 2 mfc B T 
J — OO J —OO 



k B T 

hu 


où m est la masse réduite, x la position relative (r *2 — f\).u (u étant la direction de 
l’axe de la molécule diatomique) et p la quantité de mouvement dans le référentiel local 
dirigé par ce même axe u (p est bien le moment conjugué de x). 

7 Traitement classique pour un système isolé 

Si on étudie ce système dans l’ensemble microcanonique, il n’est plus possible de 
considérer que les particules sont indépendantes, ni de séparer les degrés de liberté 
cinétique, de vibration et de rotation. 

Il n’est normalement pas possible de négliger le mouvement de rotation des molé¬ 
cules, car cette contribution est toujours énergétiquement dominante. Toutefois, à titre 
académique, et à seule fin d’alléger le calcul, on va déterminer fî pour un ensemble de N 
molécules diatomiques dans un volume V dont on ne considère que les degrés de liberté 
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cinétique et de vibration. La démarche est similaire à celle adoptée pour le gaz parfait, 
on utilise directement la condition 'H < E, ce qui donne 


n 


v 


N 


N\h 3N hrN fn 


( 2M + 2m + 2 maJ 


d 3N Pd N pd N x 


<E 


V N {y/2MË) 3N {y/2 ^Ë) JV ( v / 27(rmÜ2)) Ar 
N\h 3N 

V N {2Ef N M 3N 


où Pi est la quantité de mouvement associé au barycentre de la molécule i, tandis que 
Pi est la quantité de mouvement relative liée à la vibration de cette même molécule i, 
m est la masse réduite et M la masse totale de la molécule. Attention que cette formule 
n’est donnée qu’à titre de modèle et doit être améliorée pour toute application réaliste. 

ô Traitement quantique dans un bain thermostatique 

On se place à nouveau dans l’ensemble canonique. Les niveaux quantiques de vibra¬ 
tion sont ceux de l’oscillateur harmonique, e n = huj(n + 1/2) avec n € N. Il vient 


OO 

2 v ib = ^ ^ 6 

n =0 


-/3hw(n- |-I) _ 


1 — e 


huj 

e 2k B T 


hui 

e k B T _ l 


On en déduit l’énergie moyenne, à la température T, 


(e> 


hu 


i 

g fihw _ ^ 



Après calcul, la contribution à la capacité calorifique du mouvement de vibration pour 
un ensemble de N molécules vaut 


C v v ïb{T) 


Nk B 



Qp-ih/T 

(e'Aib/'T _ 1)2 


et fait apparaître la température caractéristique: T v ;b = frw/k B . Le comportement ob¬ 
tenu se caractérise par le gel d’un degré de liberté. En effet, pour T 2> T/ib, on 
obtient C v vi b(T) ~ Nk B car e Tvib ^ T ~ 1 + T v ib/T ; c’est le résultat classique du théo¬ 
rème d’équipartition, car il y a exactement 2N degrés de liberté. En revanche, pour 
T < T v ib, la capacité calorifique s’annule de façon exponentielle: 

C vvib (T) ~ Nk B e~ T ^ T . 


Autrement dit, pour des températures T <C T v ib, la température n’est pas suffisante par 
rapport à la barrière hw du spectre quantique pour créer des excitations thermiques, 
et les fonctions d’onde associées au degré de liberté de vibration demeurent dans l’état 
fondamental. Les excitations thermiques étant supprimées, il n’y a pas de contribution 
à la capacité calorifique et on dit que le degré de liberté est gelé. 
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5 Densité spatiale d’un gaz parfait 
en présence d’un champ extérieur 

On veut calculer la densité spatiale moyenne ( n{r )) lorsque les particules du fluide 
sont soumises à une énergie potentielle U(r) (par exemple, l’énergie potentielle d’un 
champ de gravitation, l’énergie d’une charge électrique dans un champ électrique). On 
a N particules de masse m et on néglige les interactions entre particules. L’hamiltonien 
du système s’écrit alors : 

/ || -* ||2 

H(qi,...,p N )=J2 + 

i=l ' 

D’après la loi de Boltzmann et après avoir intégré les impulsions comme dans (16), la 
densité de probabilité de trouver la première particule en r s’écrit : 

p(f) oc e~P u ^ . 

Pour la densité de particules, on obtient 

e -PU(r) 

{n(r)) = N —7. - 

/ e-^Vg 

et on peut vérifier qu’on retrouve bien le résultat attendu dans le cas d’une densité 
uniforme égale à IV/volume lorsque U(f) = 0. L’application la plus directe est celle de 
la densité d’un gaz dans une atmosphère isotherme : ( n(z )) = riQe~ m9Z ^ kBT avec uq la 
densité au niveau du sol, z l’altitude et g l’accélération de la pesanteur. 

6 Autres ensembles 

Il est possible de considérer des situations physiques différentes de celles traitées 
jusqu’à présent. 

On a représenté sur le schéma de la Fig. 7 quelques-uns des potentiels thermody¬ 
namiques qu’on rencontre, selon le choix des paramètres macroscopiques qui sont fixés. 
Dans cette figure, U est l’énergie interne, F l’énergie libre, FL l’enthalpie, G l’enthalpie 
libre et J le grand-potentiel. 

Lorsque le système est fermé, on a dN = 0. Selon les conditions imposées au système 
étudié, ces différents potentiels thermodynamiques sont minimaux à l’équilibre. Par 
définition, un système fermé est dans la situation microcanonique quand le volume et 
l’entropie sont fixés (dV = dS = 0). Alors, U est extrémale, soit dU = 0 ; on peut 
montrer en outre que U est minimale à l’équilibre (ce qu’on appelle parfois le principe 
d’énergie minimale pour les systèmes fermés, qui est strictement équivalent au principe 
d’entropie maximale). 

De même, toujours pour les systèmes fermés, on peut montrer que les potentiels, 
resp. F, H et G, sont minimaux à l’équilibre lorsque, resp. (T, V), ( S,p ) et ( T,p ), sont 
maintenus constants. Ces résultats sont particulièrement utiles pour discuter les équi¬ 
libres chimiques. J joue le même rôle pour les systèmes ouverts où un échange de matière 
avec l’extérieur est possible. 
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(S,T) (N.m) 

U(S,V,N) —-► F(T,V,N)=U-TS > 


(V-p) 


(V-p) 


(S,T) (N.m) 

H(S,p,N)=U+pV --► G(T,p,N)=F+pV - 


j(T,V,M)=F-MN 

1™ 

—► 0 


Figure 7 - Exemples de potentiels thermodynamiques (une flèche représente une trans¬ 
formée de Legendre par rapport aux variables conjuguées écrites entre parenthèses). 

a Relation de Gibbs-Duhem 

Sur la Fig. 7, on a écrit 0 pour le potentiel thermodynamique associé à l’ensemble 
décrit par les variables non fluctuantes ( T,p,p ). À la place, on aurait dû introduire une 
énergie, notons-la provisoirement M, 

M(T,p,pi) = U-ST + pV-J2^ N °‘ • 

a 

Or, M(T,p,p) est nulle parce qu’elle est à la fois extensive et intensive. 

a Extensivité de M 

M est extensive, comme toutes les énergies. Si on regarde plus en détail, tous les 
produits TS, pV, pN sont extensifs, ainsi que U. 

(3 Intensivité de M 

M(T,p, p) est intensive car elle ne dépend que de variables intensives (T, p. p a ) par 
construction. Or, toute combinaison multiplicative ou additive de quantités intensives 
est intensive. 

7 Nullité de M 

Pour conclure, M(T,p,p ) est nulle. Faisons l’expérience de pensée suivante: on 
double la taille du système. Aucun des paramètres T, p, p n’est affecté par ce double¬ 
ment, donc M(T,p, p) est inchangée (car c’est un paramètre intensif). Mais, en tant 
que paramètre extensif, on doit changer M —>• 2 M. Donc M = 2 M M = 0. 

S Identité forte de Gibbs 

On a prouvé l’égalité 

U = ST - pV + V»N a . (27a) 

a 

En conséquence, on trouve, pour un système simple, 3 

J = -pV . (27b) 

3. C’est également vrai pour un mélange, à condition de définir un ensemble grand-canonique où 
toutes les quantités N a fluctuent, ce qui n’est pas toujours le cas. 
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b Identité faible de Gibbs 

Si on calcule la différentielle de (27a) en y substituant (12), on obtient une identité 
différentielle très intéressante : 

SdT - Vdp + Yl N a dti a = 0 . (27c) 

a 

On l’applique notamment pour l’étude des mélanges et des réactions chimiques. Il faut 
l’utiliser avec discernement mais cette relation est moins risquée à utiliser que (27a). 
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7 Exercices 

a Adsorption d’un gaz à la surface d’un solide 

Un récipient de volume V contient un gaz parfait de molécules monoatomiques, en 
contact avec un solide dont la surface peut piéger des atomes du gaz. Chaque site de cette 
surface peut adsorber un seul atome et passe alors d’un état d’énergie nulle à un état 
d’énergie —£o- Le système gaz+surface est en équilibre à la température T, tout comme 
le gaz parfait, considéré séparément. La surface constitue un système, appelé phase 
adsorbée, constitué de A sites d’adsorption ; elle est également à l’équilibre, le nombre 
N d’atomes adsorbés n’est pas fixé mais le potentiel chimique p et la température T 
sont fixés par le gaz. 

i On calcule H la grande fonction de partition à partir de la fonction de partition 
canonique. 

(a) Calculer la fonction de partition canonique Z (A, T, N) de la phase adsorbée ; 
attention, on a N < A. 

(b) En déduire la grande fonction de partition E(A, T, p) de la phase adsorbée. 

ii On va retrouver 5 par une autre approche. 

(a) Calculer la grande fonction de partition Ç(T, p) associée à un seul site. 

(b) En utilisant le fait que les sites sont indépendants, montrer que 5 se factorise 
et retrouver le résultat précédant. 

iii Calculer le nombre moyen d’atomes adsorbés (N) en fonction de £q, p. A et T. 
En déduire la probabilité d’occupation d’un site, notée 9 = (N) /A. 

iv Le potentiel chimique p est fixé par celui du gaz parfait. On notera M le nombre 
d’atomes du gaz parfait et on supposera que Al 3> N, de sorte qu’on pourra 
considérer M constant. Exprimer le nombre moyen d’atomes adsorbés (N) en 
fonction de la pression p du gaz, du nombre de sites A et de la température T. 
On posera 

MT) =kaT y v v«/*.v 

Exprimer 8 en fonction de p et de po(T). 

v Tracer les isothermes de Langmuir, c’est-à-dire les courbes 9{p). 

vi Citer des situations dans lesquelles un gaz est effectivement adsorbé sur une 
surface et donner l’ordre de grandeur typique de £q. Voyez-vous les inconvénients 
de la description simplifiée présentée dans cet exercice? 

b Contribution de la rotation à la capacité calorifique 
des gaz diatomiques 

On poursuit la démarche commencée au 4 b pour les degrés de liberté de vibration. 
On décrit ici la rotation d’une molécule diatomique qu’on suppose linéaire et dont on 
néglige le moment d’inertie autour de l’axe de la molécule. Pour des raisons qu’on ne 
justifiera pas, on traitera seulement le cas de molécules hétéroatomiques (c’est-à-dire 
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que les deux atomes sont différents) et on note I le moment d’inertie de la molécule 
par rapport à un axe passant par son centre de masse et perpendiculaire à l’axe de la 
molécule. Si on note L le moment cinétique de rotation, la contribution de la rotation 


à l’énergie individuelle s’écrit e rot 



La mécanique quantique nous enseigne que, 


pour ce problème particulier, les énergies possibles sont 


£ l ~ + 1 ) ; 


avec I > 0 un entier, et qu’il y a pour chaque valeur de I un nombre 2£ + 1 d’états 
possibles qu’on peut indicer par un entier rri£ tel que — £ < rnp < l. 

i Justifier que la fonction de partition à une particule associée au degré de liberté 
de rotation s’écrit 

OO 

Zrot = + l)e^^ £(m) 

£=0 

et déterminer la température caractéristique T ro t ■ On ne peut malheureusement 
pas calculer z rot analytiquement. 

Limite basse température T <C T rot 

ii On peut ne conserver que les deux premiers termes dans z TOt . Montrer que la 
capacité calorifique associée à ce degré de liberté se comporte alors comme 

C r v °\T) = Nk B X12(|) e -2T rot /T 


iii Tracer qualitativement ce comportement. Comment appelle-t-on ce phénomène? 

Limite haute température T 3> T rot 

iv Dans ce cas, justifier que l’on peut approximer la somme dans z IO t par une 
intégrale dans laquelle t devient une variable continue. 

T 

v En faisant le changement de variables y = £(£+1), montrer que z ro t — — — dans 

-Dot 

cette limite. 

vi En déduire la valeur de C\ , ot dans cette limite. 

Ordres de grandeur et comportement de la capacité calorifique 

vii Quel est la dimension du moment d’inertie II Pour la molécule LQ, on a T rot — 
85.3 K. Justifier cet ordre de grandeur après avoir estimé I. 

viii Justifier que, pour une molécule diatomique dont on prend en compte l’énergie 
de translation, de rotation et de vibration, la capacité calorifique totale est la 
somme des trois contributions : 

C V (T) = C„ rans + C r v °\T) + C v v lh (T ) . 

ix L’air est principalement composé de diazote N 2 dont les températures caracté¬ 
ristiques sont T ro t — 2.9 K et T v ib — 3374 K. Tracer qualitativement la courbe 
C V (T). À température ambiante, quelle valeur approchée peut-on utiliser? 
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c Équilibre chimique 

Une réaction chimique quelconque peut s’écrire sous la forme symbolique : 

^ ' v a A a — 0 . 

OL 

Dans cette expression, A\, ..., A a , ..., représentent les différents réactifs ou produits de la 
réaction, tandis que v\ , ..., v a , ..., sont les coefficients stoechiométriques correspondants 
qui peuvent être positifs ou négatifs. Par exemple, la réaction 2H 2 + 02 = 2 H 2 O s’écrit 
aussi 2 H 2 + O 2 — 2H 2 0 = 0 , et donc ur 2 = 2 , vç> 2 = 1 et vr 2 0 = — 2 . 

i On note N\, ..., N a , ..., le nombre de molécules de type A\, ..., A a , ..., comprises 
dans le volume considéré, et m, ..., pL a , • ••, les potentiels chimiques correspon¬ 
dants. Exprimer dF, la différentielle de l’énergie libre. 

ii Que devient dF lors d’une réaction chimique à volume et température constante ? 

iii On rappelle que l’énergie libre est minimale à l’équilibre thermodynamique. 
Calculer dF/dNi (toujours à volume et température constantes) et en déduire 
que les réactions chimiques à l’équilibre doivent satisfaire la relation : 

^ V Q +a = 0 . ( 28 ) 

a 


Loi d’action de masse 

On note z a , la fonction de partition canonique du réactif A a . Les N a molécules sont 
supposées sans interaction, indiscernables et telles que N a 3> 1. 

iv Déterminer la fonction de partition canonique Z a des N a molécules de type A a . 

v En déduire que l’énergie libre F a du réactif A a s’écrit : 

F a = k B TN a ( ln — - 1 

V z a 

où T est la température. 

vi En déduire le potentiel chimique fx a du réactif A a . 

vii Utilisez la réaction d’équilibre chimique (28) pour établir la relation : 




viii Chaque molécule de type A a est caractérisée par ses degrés de liberté externes 
(position et impulsion) et internes (les niveaux d’énergies e an accessibles aux 
électrons qui la constituent). Calculer explicitement la fonction de partition z a 
en considérant les degrés de libertés externes comme classiques et les degrés de 
libertés internes comme quantiques. Le résultat doit faire apparaître le volum.e, 
la température, les énergies e an et la longueur de de Broglie A. a . 

ix Justifier qu’on retienne seulement l’énergie de l’état fondamental £ a Q dans le 
calcul de z a , lorsque la réaction chimique a lieu à la température ambiante. 
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x Démontrer, avec cette approximation, la loi d’action de masse: 

l[ rtt =K(T) , (29) 

a 

où n a est la densité de molécules d’un réactif A a , et K = Y\ a e ~^ Va£o10 /une 
constante qui ne dépend que de la température k B T = 1//3 et des réactifs. 

Formule de Saha 

On se propose d’utiliser la loi d’action de masse pour déterminer les taux d’ionisation 
de l’hydrogène atomique dans divers gaz astrophysiques à l’équilibre thermodynamique. 
La principale réaction est l’équilibre ionisation-recombinaison qui s’écrit : 

H + e" = H + + 2e" 

xi Déterminer les coefficients vh, Vh+ e t v e~- 

xii On considère que la masse de l’atome H et de l’ion H + peuvent être confondues. 
Etablir la formule de Saha : 

= C (k B T) 3 / 2 e~ u ^ k b t ) , 
n H 

où C est une constante que l’on précisera, et Ui, une énergie caractéristique de 
la réaction que l’on définira. 

xiii Dans le cas d’une nébuleuse gazeuse, on a T = 10 4 K et hh = 10 12 m" 3 , tandis 
que pour la photosphère du soleil, T = 6000 K et = 10 23 m" 3 . Dans les deux 
cas, la neutralité électrique globale du gaz impose n H + ~ n e _. 

Calculer le degré d’ionisation n^+ /{nu + njj+) dans les deux cas et commentez 
les résultats obtenus. 

Pour l’hydrogène Ui = 13.6 eV. On donne également, m e = 9,1 10" 31 kg, la 
masse de l’électron, e = —1,6 10" 19 C, la charge de l’électron, h = 6.63 10" 34 Js, 
la constante de Planck, et k B = 1-38 10" 23 J/K la constante de Boltzmann. 

d Compressibilité et fluctuations du nombre de particules 

Une boîte de volume V et maintenue à la température T est en contact avec un 
réservoir de particules au potentiel chimique /x. Les particules ont une masse m et 
n’interagissent pas, elles forment un gaz parfait. 

i Calculer la grande fonction de partition S (T, V,/x). 

ii Quelle est la probabilité d’avoir N particules? De quelle loi s’agit-il? 

iii Quel est le nombre moyen (N) de particules dans le récipient ? 

iv Exprimer la pression p en fonction des variables p et T. Retrouver l’équation 
d’état. 

v Calculer les fluctuations AN du nombre de particules et les relier à (N). 
Vérifier que la relation générale AN 2 = k B T ^ N y \T est vérifiée, avec 

Xt =-f tt~ ) ) compressibilité isotherme. Commenter. 

v \ dp ) N,T 
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Chapitre VI Coefficients de transport 


L’objectif de ce chapitre est de donner une interprétation microscopique des 
phénomènes de transport qui se produisent lorsqu’un fluide est mis légèrement hors 
d’équilibre. En particulier, on verra que les coefficients de transport apparaissant 
dans les lois macroscopiques sont reliés entre eux et on discutera leur dépendance en 
fonction de la pression et de la température. 

1 Relations entre flux et gradients : lois macroscopiques 

À l’équilibre thermodynamique, les variables intensives sont uniformes, c’est-à- 
dire qu’elles ne dépendent pas de la position r au sein du fluide. Si l’on induit une 
inhomogénéité de température (resp. de densité, de vitesse), caractérisée localement 
par un gradient de ces grandeurs intensives, cela génère un flux d’énergie (resp. de 
matière, d’impulsion) qui s’oppose au gradient et tend à ramener le système à l’équilibre 
thermodynamique. On s’attend de façon générale à avoir donc une relation entre flux 
et gradient de la forme 

{ énergie ( température 

matière = — Coefficient x Gradient de < densité 
impulsion ^ vitesse 

Un coefficient, appelé coefficient de transport, apparaît naturellement dans ces re¬ 
lations pour des raisons dimensionnelles et nous allons essayer de les évaluer microsco¬ 
piquement. 

a Notion d’équilibre thermodynamique local 

En présence d’un gradient d’une quantité intensive, le système n’est plus à l’équi¬ 
libre mais hors d’équilibre. Comme en mécanique des fluides, si les gradients ne sont 
pas trop forts, c’est-à-dire que les variations des quantités intensives se produisent sur 
des échelles de longueur grandes devant les distances intermoléculaires, on peut intro¬ 
duire une échelle de longueur intermédiaire mésoscopique telle que, localement, on ait 
suffisamment de particules pour définir un système thermodynamique et lui attribuer 
une densité, une pression, une température, mais qu’aux distances grandes devant cette 
échelle, ces quantités varient avec la position r ainsi qu’avec le temps t. Au final, on 
définit à l’échelle macroscopique, des champs de densité n(f,t), de pression p(r, t). de 
température T(r, t). 

Dans la suite, on se limitera à une description unidimensionnelle pour simplifier la 
discussion et on notera x la coordonnée spatiale correspondante. 
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b Définition des coefficients de transport 

Les relations flux-gradients et les coefficients associés sont définies de la façon sui¬ 
vante : 


Transport de 

Loi de 

Relation 

Coefficient 

Unité 

chaleur 

Fourier 

i-O 

II 

1 

conductivité thermique À 

Wm -1 K -1 

matière 

Fick 

jn ~ Dq™ 

constante de diffusion D 

2 —1 
ms 

impulsion 

Newton 

jp = 

viscosité dynamique q 

Pa-s 


c Équations de transport 

On établit ici les équations d’évolution des champs en utilisant les lois de conserva¬ 
tion (ou bilan) appliquées à un élément fluide. 

a Bilan d’énergie 

On considère un élément fluide entre x et x + dx de section S. On note de plus p 
la masse volumique du fluide, supposée uniforme, et c sa capacité calorifique massique, 
également supposée uniforme. L’élément fluide reçoit un flux de chaleur j q (x, t) venant 
de la gauche et perd un flux de chaleur j q (x + dx, t ) partant à droite. En appliquant le 
premier principe de la thermodynamique à l’élément fluide entre les instants t et t + dt, 
on peut écrire 


dU = j q (x, t)Sdt — j q {x + dx, t)Sdt 


--^-Sdxdt = \-—rSdxdt 
dx dx z 


après avoir utilisé la loi de Fourier. D’autre part, la température dans l’élément fluide 
va varier de T{x,t ) à T{x,t + dt). En l’absence de travail des forces de pression, on peut 
écrire que la variation d’énergie interne sera donnée par 

dT 

dU = cpSdx(T{x, t + dt) — T(x, t)) = pc-^Sdxdt . 

masse 


En combinant les deux équations, on obtient l’équation de la chaleur pour le champ 
de température 


ÔT _ d 2 T 
dt K dx 2 


avec k = — le coefficient de diffusion thermique 
pc 


Il s’agit d’une équation de diffusion. 


(30a) 
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(3 Bilan de matière 

On effectue le même raisonnement dans la situation d’une densité de particules 
inhomogène en faisant le bilan de matière au sein de l’élément fluide entre t et t + dt. 
La variation du nombre de particules est donnée par la différence des flux entrant et 
sortant 

(J j (j- yi 

dN = j n (x,t)Sdt — j n (x + dx,t)Sdt = -—d-Sdxdt = D——Sdxdt , 

dx ox z 

en utilisant la loi de Fick. On peut également écrire cette variation comme 

dïi 

dN = Sdx(n(x, t + dt) — n(x, t )) = —Sdxdt . 

Cela donne l’équation de la diffusion de la matière sous la forme 


dn ^<9 2 n 
dt dx 2 


(30b) 


7 Solution typique de l’êquation de diffusion 

La solution de (30b) dans le cas où Nq particules sont injectées en x = 0 à t = 0, ce 
qui correspond à la condition initiale n(x,0) = ^-ô(x), est donnée par: 


n(x, t) 


NO 1 r ~x 2 /4,Dt 

S y/4^Dt 


Il s’agit d’un profil gaussien centré sur x = 0 et de largeur typique caractérisée par 
l’écart type 

a x = VADt , 

qui varie en racine carrée avec le temps. On parle de transport diffusif pour le distin¬ 
guer d’un transport balistique qui correspond à un déplacement à vitesse constante 
du type x = vt, linéaire avec le temps. On peut vérifier qu’à tout temps, le nombre de 
particules injectées reste conservé 


n(x, t)dx 


No 
S ' 


Enfin, on notera le lien avec le mouvement brownien et les marches aléatoires 
étudiés en exercices, qui permet une interprétation probabiliste et microscopique de ce 
phénomène de diffusion. 

<5 Bilan d’impulsion 

On cherche à établir une équation de diffusion similaire aux deux précédentes dans 
le cas d’une vitesse transverse inhomogène pour faire apparaître l’équation d’évolution 
gouvernant le champ de vitesse. Pour cela, on se place dans une situation d’un écou¬ 
lement fluide simplifiée de telle sorte que seul l’effet de la viscosité soit pertinent. Ce 
type d’écoulement s’appelle l’écoulement de Rayleigh. Une plaque semi-infinie en 
x = 0 est mise en mouvement à t = 0 à la vitesse vq constante et dirigée selon Oz. 
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Les éléments fluides en contact avec la paroi sont entraînés par celle-ci et on a donc la 
condition aux limites v z (x = 0, t) = vq. On considérera que loin de la plaque, à l’infini, 
la vitesse est uniformément nulle, soit v z {oo,t) = 0. Le problème est invariant selon les 
directions y et z si bien que vitesse transverse ne dépend que de x et de t. On suppose le 
fluide visqueux mais incompressible, auquel cas on utilisera l’équation de Navier-Stokes. 
Il n’y a pas initialement d’écoulement selon y et par symétrie, on aura toujours v y = 0. 
Restent les composantes v x (x,t) et v z (x,t). La conservation de la masse pour un fluide 
incompressible donne 


div v = 0 , 


soit 


dv x 

Ôx 


= 0 , 


puis v x (x,t) 


v x (0,t) = 0 , 


la dernière égalité étant due au fait que la vitesse doit être tangente à la plaque à 
l’origine. Ecrivons maintenant l’équation de Navier-Stokes pour la composante v x (x,t), 
le champ de pression étant noté p(x,t ) : 


p(dtv x + v x d x v x ') 


d x p -I- yd xx v x 


v x =0 


d x p = 0 


p = const. 


On en déduit, pour la composante v z de la vitesse, que l’équation de Navier-Stokes se 
simplifie selon 

p{d t v z + uïâjrO = ~'d^+ rjdxxVz 

et correspond alors à une équation de diffusion avec, pour coefficient de diffusion, la 

viscosité cinématique v 


dv z _ d 2 v z 
dt U dx 2 


n 

avec v = — . 

P 


On peut alors résoudre cette équation pour obtenir le profil de vitesse. On obtient alors 
l’apparition d’une couche limite au voisinage de la plaque, qui s’étend comme la racine 
carrée du temps. 

Il reste cependant un problème concernant la quantité conservée lors de la diffusion. 
En effet, l’équivalent du nombre total de particules dans la diffusion de la matière 
correspondrait à 

v z (x,t)dx 

mais cette quantité n’est pas reliée aux conditions aux limites. En revanche, si l’on 

, . . dv z 

introduit la vorticite w = ——, on peut écrire 

dx 



w(x , t)dx = v z (oo) — Vq = — vq, et 


dw 

Ht 


d 2 


w 


= V 


dx 2 


Ainsi, la quantité qui est injectée dans le système et qui va diffuser est plutôt la vorticité. 


E3?Remarques : dans les bilans précédents, on n’a pas inclus de termes de sources 
ou de termes absorbants qui peuvent apparaître dans certaines situations physiques. Il 
suffit alors de rajouter ces processus dans le raisonnement ce qui conduira à des termes 
supplémentaires dans l’équation finale. 
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2 Evaluation des coefficients de transport par une approche 
microscopique 


On cherche maintenant à décrire les relations entre gradient et flux au niveau mi¬ 
croscopique pour obtenir une expression des coefficients de transport en fonction des 
grandeurs naturelles intervenant dans les systèmes fluides. 

a Paramètres microscopiques pertinents 

On rappelle les grandeurs pertinentes pour décrire un gaz, qu’on supposera parfait 
(limite de faible densité) à l’échelle microscopique : 

• densité : n ; 

• masse d’une molécule : m ; 

• température : T ; 

• énergie moyenne d’une molécule : e = ; 

• pression : p = nk^T ; 

• vitesse typique : v = = fm ; 

• section efficace de collision : a = ix a 2 où a est le diamètre des molécules ; 

• libre parcours moyen : t ; 

• temps moyen entre deux collisions : r = i/v . 

• charge élémentaire des électrons: q e = —1,6 10~ 19 C. 

b Retour sur le libre parcours moyen 

On a rappelé, à la section I 2, l’argument permettant d’obtenir la relation entre 
libre parcours moyen, section efficace et densité : 



On notera que ce libre parcours moyen ne dépend pas de la vitesse typique v et qu’il 
varie en fonction de la pression et de la température selon 


£ ~ 


T 

P 


Remarques : dans le modèle considéré de sphères dures, on peut retrouver ce 
résultat en considérant le nombre de collisions entre molécules par unité de temps en 
utilisant la théorie cinétique (cf. le partiel de mars 2014 en L3 Mécanique). 

c Constante d’auto-diffusion 

On présente ici une méthode heuristique de détermination des coefficients de trans¬ 
port par une approche microscopique. Le but n’est pas d’avoir une approche très ri¬ 
goureuse mais de se faire une idée du mécanisme de transport pour voir comment ces 
coefficients dépendent des propriétés du gaz. On ne calculera pas les préfacteurs exacts. 

Prenons la situation la plus simple, celle de l’auto-diffusion de particules d’un gaz. 
Le courant à la position x est la résultante du flux de particules se trouvant à gauche 
et allant à droite moins celles se trouvant à droite et allant à gauche. Ces particules se 
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situent à une distance typique de l’ordre de I de chaque côté car elles ne doivent pas 
avoir subi de collisions avant de franchir la paroi fictive traversée en x. Par ailleurs, 
elles ont une vitesse typique v qu’on suppose ici indépendante de x. En revanche, on 
suppose que la densité n’est pas uniforme si bien que n(x — £) 7 ^ n(x + £) afin de faire 
apparaître une loi de Fick. Du point de vue géométrique, parmi toutes les particules se 
trouvant en x =L £. seules comptent celles qui vont selon x (facteur 1 /3) et dans un sens 
bien déterminé (facteur 1/2). Au total, on affectera un facteur l/2xl/3 = l/6 qui n’est 
pas exact mais qui rend compte du fait que toutes les particules ne contribuent pas au 
flux. On peut finalement évaluer ainsi le courant de densité 

j n = x[ün(æ -£)- vn(x + £)} . 
o 

En supposant la distance caractéristique de variation de la densité faible par rapport 
à £, on peut écrire 

dïï 

n(x ±f)~ n(x) ± —£ . 

ox 

En réinjectant cela dans la relation précédente, la loi de Fick est retrouvée 

vt dn 

Jn 7T ~n > 

3 OX 


le coefficient de diffusion vaut donc 



On verra en exercice le lien entre ce résultat et une marche aléatoire. Pour un gaz 
dilué, on obtient donc un coefficient d’auto-diffusion qui augmente avec la température 
en raison de l’augmentation de la vitesse moyenne, et qui diminue avec la pression 
D oc T 3 / 2 /p. 

d Conductivité thermique 

Le même raisonnement peut être effectué pour la conductivité thermique. Cette fois, 
on considère que la densité est à peu près homogène égale à n mais que c’est l’énergie 
des particules £ qui dépend légèrement de x en raison d’une température inhomogène. 
Le raisonnement précédent peut se résumer à: il y a un flux d’environ nv /6 particules 
venant de la gauche ayant une énergie e{x — £) et un flux d’environ nv /6 particules 
venant de la droite ayant une énergie e(x + £). La différence de ces flux conduit à écrire 


j q = -f[e(x-£)-e(x + £)]. 

Pour un gaz parfait, on peut écrire de manière générale e = c v T avec c v la chaleur 
spécifique du constituant, c v = C v /N = cm = 3&b/ 2 pour un gaz parfait monoato¬ 
mique. Si la température dépend de x, on aura donc au final 

vlnc v dT 
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ce qui fait apparaître la loi de Fourier assortie d’une conductivité thermique 

nvc v £ 

X = — ■ 

Dans notre approche, la conductivité thermique augmente comme la racine de la tem¬ 
pérature mais ne dépend pas de la pression. 

e Viscosité dynamique 


Terminons par l’étude de la viscosité : le bilan se fait maintenant sur l’impulsion et 
on suppose donc une impulsion transverse inhomogène transportée par les particules 
mv z (x ± £). Le raisonnement conduit à 


jp = —[mv z (x - £) - mv z (x + t)\ . 


soit 


j P = 


mnvi dv z 
3 dx ’ 


qui correspond à une viscosité dynamique 


V = 


mnv£ 

3 


Elle aussi ne croît a priori qu’avec la température comme y/T mais ne dépend pas de p. 

f Relations entre coefficients 


Cette approche microscopique, dont les relations sont compatibles avec un raisonne¬ 
ment d’analyse dimensionnelle, permet de voir que ces coefficients de transport ne sont 
pas indépendants les uns des autres puisque tous dépendent des paramètres microsco¬ 
piques du gaz. Ainsi, on observe que 

k DI 

— ~ c et — ~ . 

7 1 V P 

On va voir en exercice que d’autres lois remarquables peuvent être intuitées sur la base 
de ces raisonnements simples. 

g Limites de l’approche 

Cependant, l’approche n’est pas générale et reste non-quantitative. Les préfacteurs 
numériques (le 1/3) sont faux par rapport à des modèles exacts et il faut tenir compte 
d’effets physiques non pris en compte, comme le fait que les collisions dépendent de 
l’énergie et donc de la vitesse. 

De plus, ces résultats ne sont qualitativement corrects que pour des gaz dilués. Les 
phénomènes de transport et leurs équations macroscopiques sont certes valables dans les 
liquides et les solides mais les lois sont souvent différentes et plus subtiles. Par exemple, 
on sait, par expérience de la vie courante, que la viscosité d’un fluide s’abaisse avec 
la température ; au contraire, dans un solide, les processus de diffusion ont souvent 
un comportement dit activé, de la forme D(T ) = DQe~ EA ^ k ^ T où Ea est une énergie 
d’activation. 
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3. Exercices 


3 Exercices 

a Marche au hasard (équation maîtresse) 

On reprend les notations de l’exercice “Marche au hasard” du premier chapitre. On 
pose V(x,t) la probabilité d’être en x au temps t. 

i Montrer que 

V(x,t + r) = ^ V(x - a,t) + ^V(x + a,t.) . 

ii Prendre la limite continue de l’équation précédente et montrer que l’équation 
aux dérivées partielles gouvernant V(x,t) s’écrit: 

&P _ D d 2 V 
dt dx 2 

On donnera l’expression de D en fonction des paramètres microscopiques. Connaissez- 
vous le nom de cette équation et des domaines de la physique où elle apparaît ? 

iii Par des arguments dimensionnels sur l’équation, donner le comportement de la 
distance typique en fonction du temps. 

iv La solution, pour le cas où V(x,0) = ô(x), s’écrit V(x,t) = ^== e~ x2 / 4Dt . La 
comparer au résultat de l’exercice fait au premier chapitre. 

b Vidange d’un réservoir 

Un gaz monoatomique est contenu dans un récipient de volume V. Un petit trou de 
surface s est fait dans une paroi, par lequel les atomes du gaz peuvent s’échapper. Le 
trou est suffisamment petit pour qu’on puisse supposer que le gaz évolue lentement et 
reste à l’équilibre thermodynamique à chaque instant. On peut ainsi définir une pression, 
qui dépend lentement du temps. Par contre, le récipient est maintenu par un thermostat 
à une température T contante. 

i Les vitesses des atomes sont distribuées selon la loi de Maxwell. Quel est le 
nombre d’atomes sortant du récipient pendant un intervalle de temps dt ? Discu¬ 
ter la dépendance de dN/dt en fonction de la densité N/V et de la température. 

ii On suppose ici que les atomes s’échappent dans le vide. 

(a) Comment le nombre N d’atomes contenus dans le récipient varie-t-il au cours 
du temps ? 

(b) Comment la pression varie-t-elle au cours du temps ? 

(c) Soit un gaz d’hélium dans un récipient d’un litre à température ambiante. La 
surface du trou est s = 1 nnn 2 . En combien de temps le récipient perd-t-il la 
moitié de ses atomes ? 

iii Discuter le cas plus réaliste de deux enceintes, à la même température T, où 
régnent initialement des pressions différentes Pi et P 2 , en contact par un petit 
trou. 
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c Conductivité électrique a et loi de Wiedemann-Franz 

On applique un champ électrique constant £ selon la direction e x . Les porteurs de 
charge sont les électrons libres, qui se comportent comme un gaz parfait classique à 
la température T. En plus de la force électrique q£, les collisions qu’ils subissent se 
traduisent par une force de frottement —mv/r. On va évaluer la conductivité électrique 
a du milieu définie par la loi d’Ohm locale 

je = cr£ , 

où j e est la densité volumique de courant électrique. 

i Montrer qu’en régime permanent, la vitesse moyenne v x est reliée à q, m, r et £. 

2 

nq t 

ii Exprimer j e selon e x en fonction de n, q et v x . En déduire que a = -. 

m 

iii Justifier la loi de Wiedemann-Franz qui prédit le rapport universel : 

A = c 

gT 

et exprimer la constante universelle C en fonction de feg, q. 


72 


exosup.com 


page facebook 




3. Exercices 


d Suspension brownienne 

On considère N particules mésosco¬ 
piques sphériques de rayon R en suspen¬ 
sion dans un liquide dont la température 
T est fixée ; par exemple, un verre de lait, 
de l’encre,..., voir la figure ci-contre. 

Les particules, dont la densité est supé¬ 
rieure à celle du liquide, sont soumises au 
champ de gravité g = — gü z et aux chocs 
avec les molécules du liquide. On notera 
m = vq (/Opart — piiq) > 0 leur masse effec¬ 
tive dans le fluide, où on note /3p a rt et pn q 
les niasses volumiques des particules et du 
liquide et vq = AnR?/?> leur volume. On 
note enfin S la section du conteneur cylin¬ 
drique et L sa hauteur. 


Courant de diffusion 

i Expliquez pourquoi la suspension peut être décrite dans l’ensemble canonique. 
On assimilera le comportement des particules à celui d’un gaz parfait de parti¬ 
cules de masse m soumises au champ de pesanteur. Ecrire, sans la calculer, la 
fonction de partition du système. Quelle est la probabilité d’un micro-état ? 

ii Justifiez que la répartition verticale de la densité de particules est donnée par 
n(z ) = noe~ m9Z ^ kBT et déterminer no en fonction de N, S, L, m, g et k^T. 

iii Comme la densité de particules est inhomogène, un courant de particules va 
s’établir. On note Jdiff(^) le flux de particules et D le coefficient de diffusion des 
particules dans le liquide. Indiquer sur le schéma ci-dessus le sens de ce courant 
de diffusion. 

Rappelez la loi de Fick. 

iv En déduire l’expression de Jdiff(z) en fonction de D, m, g , n(z) et k^T. 

Courant de dérive 

v Ecrire la force gravitationnelle qui s’applique sur une particule. 

vi En raison de la viscosité rj du fluide, les particules subissent une force qui s’oppose 
à leur mouvement : F = —'yv avec 7 = 6ttï]R. Que vaut la vitesse limite verticale 
va, dite vitesse de dérive, atteinte par les particules ? 

vii Représentez sur le schéma précédent le sens du courant de particules associé et 
donnez son expression Jdérive(~) en fonction de m, g, n(z ) et 7 . 

viii Que se passe-t-il pour les courants de Jdiff (z) et ./dérive(-Q à l’équilibre? En 
déduire la relation d’Einstein : 


GirqR 


Z 



FIGURE 8 - Schéma d’une cartouche 
d’encre de chine. Les pigments (noir de fu¬ 
mée) sont des particules de carbone d’en¬ 
viron 50 nm de diamètre et de masse volu¬ 
mique pp art « 2 g cm" 3 . 
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Chapitre VI. Coefficients de transport 


e Transport dans un métal : coefficient Seebeck — thermocouple 



(b) 



FIGURE 9 - (a) Effet Seebeck : deux métaux de nature différente forment un circuit 
fermé. Lorsqu’on chauffe une des jonctions à l’aide d’une bougie, un courant électrique 
apparaît dans le circuit et fait dévier la boussole placée en son sein. 

(b) Schéma de principe du thermocouple. 


Les thermocouples sont à la base de nombreux thermomètres utilisés dans l’habitat 
moderne, l’industrie et la recherche. Leur fonctionnement repose sur l’effet Seebeck qui 
fut découvert en 1821 et dont on reproduit sur la figure 9(a) le montage historique. L’effet 
Seebeck correspond à l’apparition d’un courant électrique js induit par un gradient de 
température (jj-, Les porteurs de charge sont les électrons libres. On définit le coefficient 
Seebeck S par la relation : 


~dT 

js = -vS-r ■ 
ax 

i Montrer qu’une température iuhomogène entraîne un courant électrique js = 
qn 


6 


■[u(x — I) — v(x + I)\. 


ii Récrire alors js en fonction de n, q, r, m, &B et 


iii Montrer que S est une constante universelle. Donner son unité et son ordre de 
grandeur. 

En présence d’un champ électrique £ = — ^ dérivant du potentiel U, le courant 
électrique total dans le conducteur s’écrit donc j e = a£ + js- On considère un fil 
métallique, dont le coefficient S est supposé constant, ouvert à ses extrémités A et B. 
De plus, celles-ci sont maintenues à des températures différentes T ,4 et Tb- 


iv Montrer qu’il apparaît une différence de potentiel AU = Ua — Ub entre ses 
extrémités, que l’on exprimera en fonction de Ta, Tb et S. 


Le principe de base du thermocouple est de souder trois fils métalliques selon le schéma 
de la figure 9(b). Les fils CA et BD sont de même nature, avec un coefficient Seebeck 
S, tandis que le fil AB est de nature différente, de coefficient S'. 

v Sachant que les points C et D (aux bornes du voltmètre) sont à la même tem¬ 
pérature Te = Tb, montrer que le montage permet de mesurer la différence de 
température Tb — Ta- 

vi Expliquer qualitativement ce qu’il se passe dans l’expérience de la figure 9(a) ? 
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Chapitre VII Statistiques quantiques 


Dans les chapitres précédents, nous avons traité le gaz parfait classique. Par “parfait”, 
on entend un ensemble de particules dont les degrés de liberté sont indépendants les 
uns des autres. Cette indépendance se traduit par le fait que le nombre de micro-états 
accessibles du système est proportionnel à V N (plus précisément Qn(/3,N) défini dans 
Eq. (18b)), ce qui conduit à l’équation d’état pV = Nk^T. 

► Si on considère un gaz de particules identiques selon la mécanique quantique, il ne 
peut y avoir d’indépendance entre particules, même lorsqu’elles n’interagissent pas. En 
effet, la fonction d’onde du système doit être totalement antisymétrique sous l’opération 
d’échange de particules si le gaz est constitué de fermions (particules de spin demi-entier) 
et symétrique s’il s’agit de bosons (particules de spin entier). 

Cette condition, appliquée à la fonction d’onde, conduit à des relations entre les états 
que les particules peuvent occuper. Par exemple, deux fermions identiques ne peuvent 
occuper le même état quantique (principe d’exclusion de Pauli), à l’inverse du cas de 
deux bosons. 

En conséquence, un gaz de fermions ou de bosons identiques, avec ou sans interac¬ 
tion, ne peut être parfait dans le sens où, d’un point de vue quantique, il ne peut vérifier 
l’équation d’état (classique) des gaz parfaits, pV = Nk-gT. 

► Le but de ce chapitre est de se familiariser avec les méthodes générales nécessaires 
à l’étude par la physique statistique des propriétés des systèmes de particules indiscer¬ 
nables obéissant aux statistiques de Fermi-Dirac dans le cas de fermions et de Bose- 
Einstein dans le cas de bosons. De nombreux systèmes physiques font appel à ces sta¬ 
tistiques quantiques notamment en physique de la matière condensée, en électronique, 
en astrophysique, ... 

1 Systèmes quantiques sans interaction 

En mécanique quantique, les niveaux d’énergie accessibles sont les valeurs propres du 
hamiltonien qui décrit le système. Dans la situation canonique, le nombre de particules 
N est fixé, chacun des niveaux d’énergie est donc susceptible d’accueillir au plus N 
particules. Pour rappeler le fait qu’il s’agit de niveaux à N particules, on utilisera la 
notation E n jy. On peut donc écrire la fonction de partition canonique quantique Z selon 

Z (P, N) = Y e~P EnN , 

n 

la somme étant prise sur tous les niveaux d’énergie à N particules. 

Les choses ne sont pas aussi simples car les fonctions d’ondes des particules pos¬ 
sèdent des propriétés de symétrie qui dépendent du spin des particules. Nous allons 
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montrer que la factorisation de la fonction de partition n’a plus lieu dans l’espace des 
positions/impulsions, 1 comme pour les systèmes classiques, mais dans l’espace des états 
à une particule. 

En effet, lorsque les particules sont sans interaction, la fonction d’onde décrivant un 
état microscopique à N particules peut s’écrire comme un produit symétrisé (cas des 
bosons) ou antisymétrisé (cas des fermions) de fonctions d’onde à une seule particule 
(cf. le cours de mécanique quantique). En conséquence, chaque niveau n d’énergie des 
états à N particules peut s’écrire comme la somme 

En.N — j 

k 

où k énumère les différents états 2 à une particule d’énergie cp. rif. représente le nombre 
d’occupation du niveau k, nombre au plus égal à N dans le cas des bosons, et à 1 dans 
le cas des fermions, où on (re)définit N, le nombre total de particules, 

y ^n k = N . 

k 

Cette dernière contrainte est difficile à prendre en compte quand on se place dans 
l’ensemble d’équilibre canonique, il est plus facile de se placer dans l’ensemble grand- 
canonique, où fi est fixé, puis de déterminer le potentiel chimique correspondant à une 
valeur de N fixée par la relation N = — dJ/dp, . La grande fonction de partition de 
particules quantiques indépendantes s’écrit donc : 


Z(T,V,p) = J2 e ~ 0inN ~ ] = Y, e " 

N,n m,ré¬ 

introduisons la grande fonction de partition du système lorsque celui-ci est dans 
l’état quantique k. On obtient alors la forme factorisée : 

s(t,v,aO = n&. 

k 

où 4 = 

rik 

Tout se ramène donc au calcul de la grande fonction de partition de 
l’espace des états à une particule. 


1. Les opérateurs quantiques x et p ne commutent pas et il n’existe donc pas d’E.C.O.C. (Ensemble 
Complet d’Observables qui Commutent) permettant une description simultanée dans les espaces réel 
et réciproque. 

2. En général, k est un vecteur, et l’état dépend des différentes composantes ki de ce vecteur. 
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2. Statistiques classique et quantique 


2 Statistiques classique et quantique 


a Statistique de Fermi-Dirac 

Dans le cas des fermions, n k = 0 ou n k = 1, d’où = 1 + , On en déduit 

l’expression de la fonction de partition : 


S(T,V,/i) = n( 1 + =► 

-PJ(T,V,ii) = J2 ln 

l _)_ e -PFk-F) 


k 

k 



On peut obtenir toutes les autres formules thermodynamiques à partir de celle-ci. 
En particulier, par dérivation, on trouve que le nombre total de particules (N) vaut 

r) J 1 

(N) = = V —r, -r- , (32a) 

' du ^ e+P( e k~») + 1 v ' 

k 

ce qui permet, au moins en principe, de déterminer [i en fonction de (N) et T. 

Puisque (N) = Y, k ( n fc)> I e nombre moyen de particules dans le niveau k s’écrit 

(n k ) = f FD (e k ) = ^(^4) + ! ’ ( 32b ) 

que l’on appelle souvent le facteur d’occupation de Fermi. 

► La distribution de probabilité (n k ) / (N) est connue en Physique sous le nom de 
distribution de Fermi-Dirac, et donne, dans le cas de fermions libres, la probabilité 
d’occupation du niveau d’énergie e k , à température et potentiel chimique fixés. 

b Statistique de Bose-Einstein 

Pour des bosons, n k G N sans restriction. On est donc ramené à sommer la série 
géométrique 3 

c, = V e__I_ 

^ l __ ’ 

neN 

ce qui suppose, afin d’assurer la convergence de la série, que 

p dFk F) ^ 

Cela impose que le potentiel chimique des bosons sans interaction doit vérifier : 


F < eo , 


puisque e k > eo, Vk, condition qui n’était pas exigée pour les fermions. 

La fonction de partition du système de bosons identiques et sans interaction est 
donc donnée par : 


E(T,V,n) = H 


1 


1 — e~P( e k~u) 


P J {T, V, fi) = y ln fl — 


( 33 ) 


3. 1 + z + z 2 + • • • = (1 — z) 1 si \z | < 1. 
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Les nombres moyens (N) et (rik) s’obtiennent à partir de /3J(T, V, p) de Eq. (33) et 
s’écrivent comme suit : 


(N) 

= y 1 

Z-J e +(3(e k -y) _ 1 ’ 
k 

(34a) 

{rik) = /be(Q;) 

1 

e+/3(efc— y) — \ 

(34b) 


/BE(efc) représente ainsi le facteur d’occupation par des bosons identiques et libres 
d’un niveau d’énergie ep- 


► La distribution de probabilité (n?-) / (N) est connue en physique statistique sous le 
nom de distribution de Bose-Einstein, et donne, dans le cas de bosons (identiques) 
libres, la probabilité d’occupation du niveau d’énergie ep, à température et potentiel 
chimique (négatif) fixés 4 . 

3 Densité d’état des systèmes quantiques — limite continue 


Le caractère discret du spectre d’énergie d’un système quantique est une conséquence 
du volume fini V du système. L’écart entre les niveaux d’énergie varie en général comme 
E —1 , ce qui permet d’assimiler le spectre discret avec un spectre continu dans la limite 
des grands volumes. Dans ces circonstances, il est possible et souvent pratique, d’intro¬ 
duire une densité d’état d’énergie p(e) telle que p{e) de compte le nombre d’états dont 
l’énergie est comprise entre e et e + de. La densité d’état dépend du spectre d’énergie 
donc du hamiltonien étudié. 

a Densité d’états 


On définit d’abord A/"(e) le nombre de niveaux d’énergie inférieurs ou égaux à e. 

• Pour un système de gaz de particules libres dans un volume V, on va utiliser une 
condition de bord. Considérons d’abord le cas à une dimension d’espace, notons 
a la dimension du segment dans lequel les particules sont confinées. Les fonctions 
d’onde s’annulent aux bords et sont solutions de 


H 2 d 2 p _ 

2 m dx 2 ^ 


(35) 


Ces solutions sont <p(x) oc sin(kx) avec k = nir/a, n G N* et e = h 2 k 2 /(2m). À 
deux et trois dimensions, on considère un volume carré ou cubique et on note b 
et c les dimensions du volume selon les axes y'Oy et z'Oz. Les vecteurs d’onde 
sont quantifiés selon les trois axes, k x = nin/a , k y = riï'n/b et k z = n^ir/c, avec 
(ni,n-2,n3) € N* 3 , et l’énergie s’écrit e = h 2 (n 2 /a 2 + n/ 2 /b 2 +n 3 2 /c 2 )/8. 

• Pour des particules libres sur un réseau cristallin, on utilise plutôt des conditions 
périodiques. Considérons d’abord le cas à une dimension d’espace, notons a la 
période du réseau. Les fonctions d’onde vérifient <p(x + a) = (p(x) et sont toujours 

4. Le cas y = et, est singulier puisqu’alors (m) —» oo. Cette occupation macroscopique d’un seul 
état microscopique est associée à une transition, la condensation de Bose, qui est reliée au phénomène 
de superfluidité. 
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solutions de (35). Elles s’écrivent, à un facteur complexe près, p = e kx . La 
condition de continuité de la fonction d’onde et de sa dérivée impose k = 2mr/a, 
n € Z et e = h 2 n 2 /(2a 2 ). Il y a un état dans chaque volume 2tt/ a de l’espace 
des impulsions, comme on le voit sur la figure ci-dessous : 


2'ïïja 2n/a 



FIGURE 10 Impulsions d’une particule dans une boîte de taille L unidimensionelle. 


De même, à deux et trois dimensions, si on note b et c les périodes selon les axes 
y'Oy et z'Oz, les vecteurs d’onde sont quantifiés, k x = 2n x 7r/a, k y = 2n y ir/b et 
k~ = 2n z n /c, avec ( n x ,n y ,n z ) € Z 3 et l’énergie s’écrit e = h 2 (n 2 /a 2 + n 2 /b 2 + 
n 2 /c 2 )/2. 

• Quelque soit la condition de bords (libres ou périodiques), on trouve le même 
nombre d’états (chaque facteur 2 dans la quantification des impulsions donne un 
facteur 1/2 dans JV, qui est exactement compensé par le nombre N* —>• Z d’états 
associés, qui est doublé). C’est la limite thermodynamique classique. 

• Il est plus facile de visualiser le calcul de J\f (s) pour les conditions périodiques. 
À une dimension, comme k peut être positif ou négatif, le nombre d’états A/"(e) 
s’écrit 


A f{e) 


2k 

2tt/ a 


a 

TT 


( 2m \ 1/2 

w 




À trois dimensions, il s’écrit approximativement 


vol. esp. récip. 47 r(v / 2 me) 3 / 3 ( 2 s + 1) _ 8nVV2(m£) 3 ^ 2 (2s + 1) 

vol.cell. base 8h 3 n 3 /(abc ) 3h 3 


où on a utilisé V = abc et les volumes ont été calculés dans l’espace réciproque, 
qui est bien adapté parce que p est un E.C.O.C. 5 pour les particules libres. On 
a ajouté le facteur 2s + 1 de dégénérescence de spin. Le cas à deux dimensions 
est illustré sur la figure ci-dessous. On voit que le nombre d’états est exactement 
le nombre de sommets du réseau inclus dans le cercle, et donc, à quelques unités 
près, au nombre de cellules de base inclues dans la sphère car il y a quatre 
sommets de poids 1/4, soit 4/4 = 1 sommets par cellule. 



5. Voir la note 1 page 76. 
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• Pour finir, on calcule la densité d’état à une particule par la formule 

. . dAf, . . . 

p(e) = -^-(e) • (36) 

Par exemple, à une dimension, la densité d’état vaut 

, , dN{e) L f 2m \ 1/2 1 

l,{c) = — = 2i ffirj Te' 

Il ne faut pas oublier de prendre en compte la dégénérescence des états due au 
spin des particules. Le résultat précédent doit être multiplié par 2s + 1 pour des 
particules de spin s. 

• On reprend le calcul pour une particule classique libre, dont l’énergie s’écrit : 

e _Müî _ 

2 m 

Le niveau d’énergie fondamental est £q = 0. La fonction de partition canonique 
à une particule s’écrit dans ce cas simple : 


/ 

iM 3 xI 3 


a HpII 2 d 3 rd 3 p 

=> C* 2m. - 


h 3 


V 

h? 



d 3 p . 


Essayons d’exprimer cette même intégrale dans l’espace des énergies e, et ainsi 
de définir une densité d’état dans cet espace. Les deux seules relations utiles ici 
sont les suivantes : 


d 3 p = dpxdpydpz = 47 rp 2 dp ; 

de = — dp ; 

m 

on a introduit la notation p = ||p|| = (p x 2 +Py ’VVz'Ÿ^- combinant ces 
deux équations avec la relation de définition de l’énergie (e = p 2 /(2m)), on 
trouve facilement que d 3 p = Amm\/2me de. La fonction de partition s’écrit donc 
également 6 : 


2nV{2 m) 3 / 2 

h 3 


yfee~P e de 


e~P e V(e) de . 


► La densité d’états des particules classiques libres dans l’espace des énergie 
s’écrit donc sous la forme : 

, 2nV (2m) 3 / 2 r 
p( E ) = — 

(36) établie à partir de l’équation de Schrôdinger donne le même résultat au 
facteur 2s + 1 près. 

6. En utilisant J 0 °° de = 2 ^/ 2 , on retrouve bien le résultat attendu z = L/A 3 , où A est 

donné par (17). 
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3. Densité d’état des systèmes quantiques - limite continue 


b Expressions intégrales des fonctions thermodynamiques 

Une fois la densité d’état obtenue, on calcule les grandeurs thermodynamiques en 
remplaçant les sommes par des intégrales. On trouve : 


F = =F^ AfeTln(l ±e-P ek ) 

k 

poo 

-> TftBT / 

J 0 

p(uj) ln(l ± e l3hui )duj ; 

(37a) 

(N) 

II 

g 

W 

II 

poo 

/ f( £ )p( £ ) de ; 

(37b) 


k 

J eç 


U = (E) 

= ^2 e k (n k ) 

k 

p oo 

= / £ f( £ )p( £ ) de , 

J £0 

(37c) 


où /(e) représente le facteur d’occupation -la probabilité d’occupation- du niveau 
d’énergie e (selon une des statistiques mentionnées plus haut) et =t = + les fermions, 
± = — pour les bosons. Ne pas confondre p{e) de, le nombre d’états accessibles dans 
[e, e + de], et f(e)p(e) de le nombre d’états effectivement occupés à la température T, 
compte tenu de la statistique à laquelle obéissent les particules. 

c Statistique de Planck 

Les photons sont des bosons sans masse et sans interaction entre eux (gaz parfait) 
dont le nombre n’est pas fixé (ils peuvent être continuellement créés ou absorbés). En 
conséquence, les photons obéissent à la statistique de Bose-Einstein avec un potentiel 
chimique p identiquement nul. Les niveaux d’énergie accessibles aux photons vérifient : 

£k = hv = hck . 


Chacun de ces états est deux fois dégénéré (cette dégénérescence correspond aux deux 
états de polarisation du photon). Dans cette expression, h est la constante de Planck, 
h = h/ 2 tt), c est la vitesse de la lumière dans le vide, v la fréquence du rayonnement, 
et k = k ( k étant le vecteur d’onde). 

En utilisant les résultats de la section précédente, on en déduit aussitôt que le nombre 
moyen de photons (n^) dans le niveau k s’écrit : 

1 


(™k> = 


Q(3hck _ ^ 

La distribution associée est la distribution de Planck. L’énergie vaut alors 

Huj 


(£w> = 


q/3Huj _ -y 


a Loi de Planck 

On définit tt(w; T) la densité spectrale d’énergie volumique de sorte que V u{u\T)du> 
soit l’énergie électromagnétique contenue des les modes de pulsations entre uj et ui + du 
dans une boîte de volume V, à une température T. On trouve 


u(uj; T) 


h u 3 
7 r 2 c 3 — 1 


(38a) 
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Chapitre VII. Statistiques quantiques 


/3 Loi de Wien 

Le maximum d’émission d’un corps noir correspond au maximum de u et est atteint 
pour une pulsation w max proportionnelle à la température : 

u; max = 2.821^ . (38b) 

On observe que w max augmente avec la température. 

7 Énergie moyenne — loi de Stefan-Boltzmann 

L’énergie contenue dans un corps noir idéal suit la loi 


(E) = VaT 4 


avec 


k B 

15 ( hc ) 3 


une constante universelle. 


(39) 


Elle est reliée avec le fait que le nombre total moyen de pilotons varie comme T 3 
comme l’indiquent les lois de Planck et de Wien conjuguées. 

ô Pression et entropie 


On trouve que la pression du rayonnement suit la loi : p 

T , . „ 4 (E) 

L entropie vaut b = ——. 

3 T 

4 Gaz parfait de fermions 


1(E}_ 
3 V ' 


a Surface de Fermi 


On s’intéresse ici uniquement à un système de fermions. La fonction /fd devient 
une marche discontinue, à température nulle ; plus précisément, quand T = 0, on trouve 
/fd(s) = 1 pour e < //(O) et /fd(^) = 0 pour e > /r(0), où on a précisé la dépendance 
en température du potentiel chimique. Quand T est non nulle, la distribution de Fermi- 
Dirac décroît continûment, mais reste proche du cas T = 0 tant que k^T -C /i(0), 
comme on l’observe sur la figure suivante : 



- T=0 

- T=Ti>0 

—. T=T2>T, 

— t=t 3 >t 2 


a Énergie de Fermi 

On appelle énergie de Fermi l’énergie £f = p(T = 0). À température nulle, l’état 
quantique à N particules correspond à la situation suivante : tous les niveaux d’énergie 
e < Ep sont remplis et tous les niveaux d’énergie £ > ep sont vides. L’énergie de Fermi 
correspond donc au dernier niveau rempli, et comme ce sont ceux qu’on peut exciter le 
plus facilement et qui sont donc les plus fragiles, l’essentiel des mécanismes physiques 
d’un système de fermions concernent les états dont l’énergie est proche de Ef- 
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4- Gaz parfait de fermions 


(3 Température de Fermi 

On appelle température de Fermi Tp = sp/kp. C’est la température en-dessous (ou 
autour) de laquelle les effets quantiques sont significatifs. 

b Résultats à température nulle 
a Niveau et température de Fermi 

Au nombre de particules N correspondant la densité n = N/V. On trouve, pour des 
fermions libres, 

£F = (è n ) 7 = L ^ 2 ' 3 ■ 

(3 Énergie et pression quantique 

Toujours pour des fermions libres, on trouve l’énergie totale : 

(E) = ^Nep 

et la pression du gaz de fermions libres : 

2 2 

p = -ns f = —nkpflp . 

5 5 

c Développement à basse température 

a Développement de Sommerfeld 

Pour une fonction /(e) régulière autour de s = p, on peut écrire: 

J^ f(s)(n(s))ds = jT f{s)ds + y (k B T) 2 f (p) + O 

/3 Fonctions thermodynamiques 

On peut ainsi développer le potentiel chimique p(N) en fonction de la température, 
on obtient 

Mr,A0 = EF (i-î- (A) +...j . 

On en déduit pour l’énergie 

U(T, N) = 17(0, + ^ (Jr) 2 + • • • 

ce qui donne une chaleur spécifique linéaire en température si bien qu’il en sera de même 
pour l’entropie. 

Enfin, on obtient rapidement le comportement de l’équation d’état à basse tempé¬ 
rature en vertu de la relation pression-énergie 

p(T,n) =p(0,N) fl + ^- (^rj +••• 
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5. Exercices 


5 Exercices 

a Densités d’état des électrons libres 

Montrer que la densité d’état des électrons libres contenus dans une boîte de volume 
L d , où d est la dimension d’espace (d = 1,2 ou 3), s’écrit : 

V(e) = K d ë , 

où K ( i est une constante qui ne dépend que de d, L, m et h. 

b Gaz d’électrons soumis à un champ magnétique 

On considère un gaz parfait d’électrons plongé dans un champ magnétique uniforme 
et constant, B = Bü z . On note fl le moment magnétique d’un électron. La relation entre 
le moment magnétique et le spin S de l’électron s’écrit fl = gp^/hS où /tb = q e h/2m 
est le magnéton de Bohr et g ~ —2 le facteur de Landé d’un électron libre. 

i Montrer que la densité d’états p±(e) des électrons de spin f (up, s z = h/2) ou 
de spin J. (down s z = —h/2) est 

p±(e) = ^P(e=FMB B) , 

où T>(e) est la densité d’états calculée à l’exercice a à 3 dimensions pour s = 1/2. 

ii Tracer les courbes p+(e) et /0-(e) ensemble, et y placer le niveau de fermi. 

iii Calculer N + et 1V_ les nombres d’électrons de spin f et h à température nulle. 
En déduire l’expression de £p, le niveau de Fermi, en fonction de N = N + + IV_ 
le nombre total d’électrons, quand /ie B <C £f(B = 0). Vérifier que, dans ces 
hypothèses, £f est indépendant de B. 

iv Sachant que /xb = 5,8 X 10 -5 eV/Tesla, montrer que la condtion pbB <C £f est 
très bien vérifiée expérimentalement. 

v Montrer que l’aimantation M à température nulle se met sous la forme 



2 £f 

vi Comparer le résultat ci-dessus avec le résultat équivalent pour des particules 
discernables. Expliquer l’origine de l’énorme réduction de l’aimantation. 
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Chapitre VII. Statistiques quantiques 

c Naine blanche et limite de Chandrasekhar 


Généralités 

On considère un ensemble de N fermions, de masse m, de spin 1/2, sans interactions, 
compris dans un volume V = L 3 . Les niveaux d’énergie accessibles, e, ne dépendent que 
du module p = ||p|| de l’impulsion qui est quantifiée : 


( Px 
P= I Py 


2irh 2n h 2-Kh ■> 

avec p x = n x—j~, Py = %-j—, P- = Hz ^~ et \ n xi n yi n z) G £ 


Dans toute cette partie les fermions sont considérés à température nulle. 

i Donner l’expression du facteur d’occupation de Fermi-Dirac. Le représenter sur 
un schéma en fonction de p/pFf, où p-p est l’impulsion au niveau de Fermi. 

ii Montrer que la densité d’états V(p) dans l’espace des impulsions s’écrit : 

V(V) = 5 

on n’oubliera pas la dégénérescence due au spin des particules. 

iii Montrer que la densité n est donnée par l’expression suivante : 


3vr 2 h 3 ' 

i V 2 

iv L’énergie de particules non relativistes s’écrit e = j—. 

Démontrer l’expression de la densité d’énergie interne u = U/V : 


10ir 2 h 3 m . 

Exprimer u en fonction de n, m et h. 

v L’énergie de particules ultra-relativistes s’écrit e = pc, où c est la vitesse de la 
lumière dans le vide. Démontrer l’expression de la densité d’énergie interne : 


47 T 2 H 3 


Exprimer u en fonction de n, c et h. 

Données numériques utiles pour la suite 

- masse du soleil Mq = 2 X 10 30 kg ; 

- masse d’un nucléon tun = 1,67 X 10~ 2 ' kg; 

- masse de l’électron m e = 9,11 x 10 -31 kg; 

- constante de Boltzmann kp = 1,38 X 1CD 23 J K -1 ; 

- constante de Planck h = 1,05 x 10 -34 J s ; 

- vitesse de la lumière c = 3 x 10 8 m s^ 1 ; 

- constante universelle de la gravitation G = 6,67 x 10~ n Nm 2 kg -2 
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Limite de Chandrasekhar 

On considère une naine blanche. Son rayon est proche de celui de la terre, sa masse 
proche de celle du soleil. Contrairement au Soleil, les naines blanches ont brûlé tout 
leur combustible nucléaire et produisent peu d’énergie. On les assimilera ici à un milieu 
totalement ionisé comportant un nombre égal N de protons, de neutrons et d’électrons. 

vi On considère d’abord que les électrons forment un gaz de fermions libres non 
relativistes à température nulle. 


(a) En utilisant les résultats de la première partie, montrer que l’énergie cinétique 
moyenne E c des N électrons de l’étoile peut s’écrire : 


E r = C 


iV 5 / 3 

f/2/3 


où C est une constante qu’on précisera et V le volume de l’étoile, 

(b) On assimile l’étoile à une sphère de rayon R. Etablir le résultat : 


E c = C l 


N 5 / 3 

R 2 


avec C\ = 


3h 2 f 9tt\ 2/3 


10m e \ 4 / 


où m e est la masse des électrons. 

(c) L’énergie potentielle gravitationnelle de l’étoile est dominée par les contribu¬ 
tions des N neutrons et N protons. On écrit m pro ton = m neu tron = mjv la 
masse d’un nucléon. L’énergie potentielle gravitationnelle s’écrit, en suppo¬ 
sant une densité de masse uniforme au sein de l’étoile, 

N 2 12 

e p = -C 2 ^, avec C 2 = y C{m N f , 
où G est la constante de gravitation universelle. 

Déterminer le rayon pour lequel l’énergie totale E c + E p est minimale. 

(d) On néglige les énergies de cohésion, M = 2Nm^- Déterminer numériquement 
N. En déduire une valeur numérique du rayon d’équilibre. 


vii Lorsque la densité est élevée, l’analyse précédente doit être reprise en considérant 
les électrons comme ultra-relativistes. 


(a) En utilisant les résultats de la première partie, et en assimilant encore l’étoile 
avec une sphère de rayon R, montrer que l’énergie cinétique moyenne E c des 
N électrons peut s’écrire : 


N 4 / 3 

E c = C 3 —— avec 

SX 


<*=¥{%) 


1/3 


4 \4 


(b) Comparer la dépendance en R des énergies cinétique potentielle. 

En déduire qu’il existe une valeur critique de N, noté N c , au-delà de laquelle 
l’étoile s’effondre sur elle-même. Exprimer N c en fonction de h, c, m .jy et G. 

(c) Calculer la masse critique correspondant à N c et la comparer à celle du soleil. 
Cette masse théorique maximale des naines blanches est appelée la limite de 
Chandrasekhar (physicien indien, 1910-1995). 
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d Semiconducteur intrinsèque 


Un semiconducteur est un matériau 
dont la bande de valence est remplie et la 
bande de conduction vide, à température 
nulle. On représente les bandes sur la fi¬ 
gure ci-contre. Elles sont supposées semi- 
infinies. L’origine des énergies est choisie 
au sommet de la bande de valence. L’éner¬ 
gie minimale de la bande de conduction est 
notée e g et appelée énergie du gap. 

À température non nulle, l’agitation 
thermique permet à des électrons de la 
bande de valence de franchir le gap sépa¬ 
rant les bandes. Le matériau devient alors 
faiblement conducteur. 


bande dé conduction 

_ i — 



I 


bande d$ valence 

I 

| 

| 

| 

_ 1 _ 


Figure 11 - Bandes de valence et de 
conduction séparée par la bande interdite. 


Un électron migrant de la bande de valence à celle de conduction est assimilé à un 
fermion libre de spin 1/2 et de masse effective m* e . De plus, le trou laissé dans la bande 
de valence est assimilé à un fermion libre de spin 1/2, de masse effective m / et de charge 
opposée à celle de l’électron. 

i Donner un encadrement rapide de l’énergie de Fermi £p. 

ii Ecrire la densité d’états en énergie p c on d( £ ) pour les électrons de la bande de 
conduction, et son analogue p v a i(e) pour les trous de la bande de valence. 

iii Montrer que le nombre moyen de trous occupant un état d’énergie e est 

^ n ^ Q-PC-C) + i 

iv Donner les expressions générales de N e , le nombre d’électrons dans la bande de 
conduction, et de Nt, celui de trous dans la bande de valence. 

v Déterminer explicitement N e et A/, quand /3p 1 et /3(s g — p) 1. 

Vous utiliserez le résultat d’intégration J 0 °° y/xe~ x = 

vi On a N e = Nt pour un semiconducteur intrinsèque. En déduire la densité n = 
N e /V des porteurs de charge dans la bande de conduction, puis le potentiel 
chimique p(T). Situer p(T) pour vérifier les hypothèses de v. 

vii Le tableau ci-contre donne des valeurs 
caractéristiques pour le Silicium et le 
Germanium, à T = 300 K. Calculer 
p(T) puis n pour chacun de ces maté¬ 
riaux. Comparer n Sl et n Ge entre eux et 
à la densité d’électrons dans un métal 
tel le cuivre. Vérifier que les hypothèses 
du v sont bien vérifiées. 



£g(éV) 

ml/m. 

771./ /771. 

Si 

1,12 

1,13 

0,55 

Ge 

0,67 

0,55 

0,29 











chapitre vin Transitions de phase 


1 Introduction 

L’un des plus grands sujets de la physique statistique est l’étude des transitions 
de phases. On observe des phases thermodynamiquement stables, bien séparées selon 
les paramètres d’état (par exemple, T et p, mais on peut faire un diagramme à trois 
paramètres, pour un corps simple). L’existence de phases différentes a été établie en 
1942 par Onsager sur le modèle d’Ising, qui décrit un système de fermions de spin 1/2. 
Toutefois, le mécanisme de condensation de Bose-Einstein, qui implique une transition 
entre un état condensé à basse température et un état normal à plus haute température, 
est antérieur à ces travaux. 

Pour un système aussi simple que l’eau, qui existe sous forme gazeuse (vapeur), 
liquide ou solide (glace), il n’existe aucune démonstration de l’existence de phases dif¬ 
férentes, ni donc de transitions entre ces phases. Pour illustrer la difficulté, on peut 
mentionner qu’il existe des chemins, quand on se déplace dans le diagramme de phase 
( p, T ), qui joignent continûment sans transition brutale la phase liquide et la phase 
gazeuse ; il suffit de passer au-dessus du point critique, comme on le comprend sur le 
diagramme suivant : 


Pression (bar) 



La caractéristique d’une transition de phases, qui ressort de ces exemples, est donc 
l’existence de la discontinuité d’une quantité physique quand on varie un paramètre 
thermodynamique. Plus précisément, c’est le potentiel thermodynamique de l’ensemble 
d’équilibre dont on cherchera à caractériser une discontinuité, puisque toutes les autres 
quantités physiques en découlent. 
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Chapitre VIII. Transitions de phase 


Selon la terminologie inventée par Landau, on dit qu’il y a une transition du premier 
ordre, lorsque ce potentiel est discontinu, et du second ordre, lorsque une dérivée est 
discontinue. 

Il existe des systèmes, pour lesquels il n’existe même aucune transition discontinue. 
On peut citer notamment les systèmes à une seule dimension à température non nulle. 
Pour décrire les changements de phases, les théoriciens utilisent alors le terme anglais de 
cross-over. Ces transitions douces sont également intéressantes mais ne présentent pas 
les mêmes caractéristiques et ne doivent pas être assimilées à des transitions de phases. 

2 La condensation de Bose-Einstein 

Certains phénomènes physiques qui défient les lois habituelles, comme les fontaines 
spontanées d’hélium He A+ ou la supraconductivité du mercure à T < 4, 2 K, sont liés 
au mécanisme de condensation de Bose-Einstein, qu’on étudie maintenant. 

a Potentiel chimique d’un système de bosons 

La fonction /be définie en (34b) diverge quand e — >• pi. Elle n’est définie que, pour 
e > pi, sinon (n(e) < 0 a une valeur non physique. On en déduit immédiatement que 
pi < £q, où £q est l’énergie minimale, qui apparaît dans les formules (37b) ou (37c). 

De fait, on a souvent sup (pi) = £o car le potentiel chimique ne possède aucune autre 
contrainte qui le borne. Le paramètre qui permet de varier le potentiel chimique est la 
température, et, généralement, pi (T) est une fonction décroissante, ce qui implique alors 


lim u(T) = £n 

T —>0 


b Divergence du nombre de particules dans l’état fondamental 

Quand pi —>• £q, le nombre de particules dans l’état fondamental /be(^o) diverge 
tandis que les autres remplissages /be(e) pour £ > £o restent petits. 

Cela signifie que toutes les particules se condensent dans l’état fondamental, tandis 
que les autres niveaux se dépeuplent simultanément. C’est le mécanisme de condensation 
de Bose-Einstein. 

c Mécanisme de condensation 

Quand la température T devient très petite, de sorte que k B T <C £i — £o, où £i est 
l’énergie du premier état excité, alors /be(^i) ~ 0. De même, /be(^) pour toute £ > £i, 
et la somme f £> fBE{^)p(^)ds converge et reste également négligeable. 

À l’inverse, /be(^o) ~ N, ce qui permet de déterminer exactement pi. On peut, en 
particulier, faire un développement limité, pour estimer l’écart £q — pi: 

k B T 
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2. La condensation de Bose-Einstein 


d Comportement à haute température 

a Comportement du potentiel chimique dans l’approximation continue 

On suppose des bosons de spin nul pour simplifier la discussion. On montre qu’au 
dessus d’une certaine température, le nombre de particules vérifie : 


N oc T a J a avec = 


f 


e x /<p — 


i dæ = E 


n =1 


La valeur de a et le coefficient de proportionalité dépendent du modèle choisi. Par 
exemple, pour des bosons libres dans une boîte à trois dimensions, en utilisant les 
résultats sur la densité d’état et (17), on trouve a = 3/2 et 

" = • 

(3 Étude graphique et température de Bose-Einstein 

On trouve un comportement typique d’une transition du premier ordre. On re¬ 
présente l’exemple de l’oscillateur harmonique bidimensionnel sur la figure suivante, à 
laquelle on a joint le cas d’un simple cross-over, tiré de l’oscillateur harmonique unidi¬ 
mensionnel. 




FIGURE 12 - Courbes du nombre d’occupation de l’état fondamental d’un système 
de bosons en fonction de la température. La première courbe concerne l’oscillateur 
harmonique bidimensionnel : on y distingue la transition à la température critique T c ~ 
0,78 K, entre l’état normal et l’état condensé. La seconde courbe concerne l’oscillateur 
harmonique unidimensionnel : on y constate que le passage entre l’état normal et l’état 
condensé est doux, il ne s’agit pas d’une transition mais d’un cross-over. 


La température critique de transition s’écrit 


T c oc 


l/a 

C(a) 1 /» 


avec Ç (a) 


JM ), 


où a a la valeur définie à la section précédente. Pour les bosons libres dans une boîte, 
on trouve 


k B T c 


Airti 2 ( n \ 2 ^ 3 

~ïrn U(3/2)/ 
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e Comportement à basse température 


a Fraction condensée et approximations 

On sépare l’occupation du fondamental avec £o = 0 dans la somme sur les états et on 
la note IVo = (no). Il y a condensation lorsque IVo oc N, c’est-à-dire IVo est une fraction 
macroscopique. Comme /i ~ — k^T/No, on pourra prendre /i = 0 dans les équations 
avec une bonne approximation. 

On peut enfin négliger les occupations des états excités par rapport à Nq et les 
traiter dans l’approximation continue. On en déduit finalement 

w »< T >="O-(£)“)- 


où a a la valeur définie précédemment. 

(3 Comportement thermodynamique en dessous de T c 

L’énergie moyenne pour T <T C vaut 


U ={E) oc 


a+l 


où a est définie comme précédemment. Pour des bosons libres dans une boîte, on trouve 

si 


rp \ 5/2 

(E) ~ 0,77 Nk B T c ( - 
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3. Exercices 


3 Exercices 


a Gaz de bosons et condensation de Bose 

Un gaz quantique constitué de N bosons sans interaction, ponctuels, de masse m et 
de spin 5 = 0, est contenu dans un cube d’arête L et de volume V = U 3 . On note Nq 
le nombre de bosons dans l’état fondamental (p = 0) et N' le nombre de bosons dans 

/ n x \ 

les états excités ( p fi 0). On suppose que: p = ; n = \ n y et (n x ,n y ,n z ) € Z 3 . 

\nj 

On notera e(p) = ||p|| 2 /(2m). On rappelle que T(|) = fifr/2 et T(|) = et on 

définit, pour tout z € [0,1], 


9u(z) = 


1 


roo U u-1 


i» 


'0 


e 

-1 

z 


^ z n Q x. 

du = 53(1) = C(x) ~ 2,612 et gs (1) = C(x) ~ 1,314 

77/2 Z 2 Z 

n= 1 


Préliminaires 

i Rappeler l’expression du nombre moyen d’occupation (n$) de l’état p en fonction 
de e{p). En déduire que le potentiel chimique doit vérifier p < 0. 

ii Exprimer IVo et N f en fonction des 

Calcul explicite de N' 

iii On se place à la limite thermodynamique. Justifier la transformation : 


y 


d 3 p 


iv Dans cette limite, expliciter le calcul de N' sans l’effectuer. Montrer, en particu¬ 
lier, qu’il est licite de choisir la borne inférieure d’intégration p =|| p ||= 0. 

Ceci justifie a posteriori le traitement fait à part du niveau fondamental. 

v Donner l’expression de N 1 à l’aide de g 3 , k-g, T et p. 

2 

Condition de condensation 

vi Etudier la variation de g% (z) dans son domaine de définition. 

vii En déduire que N' est majoré par A^^ ax (U, A), où on explicitera la fonction A^^ ax 
en fonction de V et de la longueur d’onde thermique A = ~^== = . 

viii Montrer qu’il se produit une transition quand T < T c et préciser T c . 

ix Pour T > T Cl montrer que N 1 ~ N et p < 0. 

Étude pour T <T C 

x Montrer que l’énergie moyenne du gaz de bosons est donnée par l’expression : 

3Vk B T 5 
' E ' A 3 ^2 ' 

L’exprimer en fonction de N, T et T c . 

xi Calculer la capacité calorifique à volume constant du gaz en fonction de N, T et 
T c . Comparer ce résultat à la capacité calorifique d’un gaz classique. 
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Chapitre VIII. Transitions de phase 
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Annexe A Probabilités 


1 Combinatoire 


La combinatoire est l’ensemble des méthodes qui permettent de dénombrer de ma¬ 
nière efficace des ensembles lorsqu’il y a un grand nombre de possibilités et que l’énu¬ 
mération directe est si fastidieuse qu’elle devient impossible. Ces méthodes sont donc 
particulièrement utiles pour les méthodes probabilistes et statistiques. Elles permettent 
de déterminer le nombre total de possibilités et celui de certains types de configurations. 

a Choix multiplicatifs 

Le nombre de manière de choisir n objets pour chacun desquels on a q choix s’écrit : 

q x q X • • • X q = q n . 

Si on a deux ensembles A et B ayant un nombre d’éléments fl a et fl B, le nombre de 
couples (a, b) de A U B est x fl g. Cela sera utile pour la notion d’entropie. 

b Factorielle 


Le nombre de listes ordonnées de n objets numérotés (c’est aussi le nombre de 
permutations de ces n objets) s’écrit : 

n\ = nx (n — 1) x ■ • • x 2 x 1 . 

Ce nombre croît extraordinairement vite avec n. 

c Formule de Stirling 


Elle est valable lorsque n 3> 1 mais donne des approximations remarquablement 
bonnes pour n petit également (le premier terme correctif dans la parenthèse est l/12n) : 



iïJÊt Pour le logarithme, on utilisera souvent l’approximation suivante qui néglige les 
termes d’ordre ln n : 

lnn! ~ nlnn — n + 0(\nn) . 

d Arrangements 


Le nombre de listes ordonnées (ou arrangements) de k objets parmi n (0 < k < n) 
s’écrit : 


A 1 ^ = n x (n — 1) x • • • x (n — k + 1) = 


n: 


(n — k)\ 
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Annexe A. Probabilités 


e Combinaisons 


Le nombre de listes non ordonnées (ou combinaisons) de k objets parmi n (0 < k < n) 
s’écrit : 


ni 


k\(n — k)\ 


Une autre manière de voir les choses est de dire qu’on ne distingue pas les objets qui 
sont tirés. Ce sera le plus souvent le cas en physique. 


f Coefficients binomiaux 


Quelques relations utiles : 


V ( n \ = ‘2 n • ( n \ = ( U V (n\ = n(n- 1\ 

^ W ’ W \n ~ k) ’ \k) k\k-lj- 

Pour les coefficients binomiaux, on suppose n, k 2> 1 et n — k 3> 1, on pose k = nx si 
bien que 0 < x < 1. Alors, on peut écrire: 

~ n [—xhvx — (1 — x) ln(l — x)] . 

Une autre formule utile : si on imagine que n> 1 mais k est fixé, on peut écrire 


k\ 


_ k -1 

quand k <C y/n , en utilisant ^ i = /c(A: — l)/2. 

i=0 

On rappelle la formule du binôme : 

{p + q) n = {p + q) x (p + q ) x • • • x (p + q) 



n 


n 

ÂT 



+ o 


rr 


n 

ÂJ 




2 Quelques rappels sur les probabilités et les distributions 


a Définitions 
a Issue 

Soit une expérience. Son résultat est dit aléatoire s’il y a plusieurs résultats possibles. 
Chaque réalisation, que l’on notera i, est appelée une issue et on lui attribue une 
probabilité P*. 

Les probabilités vérifient, par définition, les deux propriétés suivantes : 

0 < Pi < 1, V i et ^ Pi = 1. 
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2. Quelques rappels sur les probabilités et les distributions 


/3 Évènement 

On appelle évènement un ensemble de résultats. Par exemple, si on lance un dé à six 
faces numérotées de 1 à 6, l’évènement “le résultat est pair” correspond à l’ensemble des 
réalisations i = 2, i = 4, i = 6. Une issue est également appelée évènement élémentaire. 

Un événement est dit impossible si sa probabilité est nulle, certain si sa probabilité 
est 1. 

Soit un évènement A. on associe à cet événement sa probabilité de réalisation P a 
définie par la somme des probabilités Pj des issues favorables à la réalisation de A, 

pa= Y p ; 

i favorable à A 


« la probabilité d’un événement est la somme des probabilités des issues favorables. » 
Si chaque issue est équiprobable, on trouve 


P = 


1 N a 

N Ct Pa = ~N ' 
où Na est le nombre d’issues favorables à A et N le nombre total de réalisation. 


7 Probabilités composées 

Dans l’espace des évènements possibles £1, on peut également définir des évènements 
composés qui correspondent à un sous-ensemble de £. Un exemple très simple est celui 
d’un jeu de dés : les évènements élémentaires correspondent à chacune des faces possibles 
tandis que l’évènement “un chiffre pair est sorti” est un évènement composé. La théorie 
des ensembles et ses notations sont très utiles pour les probabilités. L’évènement U A et 
B ” se note Ad B, U A ou B ” se note AU B. On trouve le résultat fondamental suivant : 


P(A U B) = P{A) + P(B) - P(A n B) 


(Al) 


Le cas où P(A 0 B) = 0 correspond à la situation où les événements A et B sont 

incompatibles : 


P(AuB) = P(A) + P(B) 


(A et B incompatibles) 


ô Probabilités conditionnelles 

On peut considérer des évènements tels qu’ils soient conditionnés par d’autres évè¬ 
nements : quelle est la probabilité d’avoir A sachant que B s’est produit? On note A\B 
('“A si B ”) cet évènement (attention: A\B B\A). En revenant à la définition des 
probabilités, on trouve la loi de composition : 


P{A\B) 


P{AdB ) 
P(B) 


[ïrrOn en déduit pour les évènements composés : 


P(A 0 B) = P(B)P{A\B) = P(A)P(B\A) . 


La dernière égalité correspond à la formule de Bayes : 


P(B\A) 


P(B) 

P(A) 


P(A\B) . 
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Annexe A. Probabilités 


e Indépendance des événements 

Il se peut que la probabilité de A ne dépende pas de la réalisation de B. On dit alors 
que les événements A et B sont indépendants. Dans ce cas, on a : 

P(Afl B) = P(A) P (B) (A et B indépendants) . (A 2) 

Dans ce cas, P(A\B) = P (A), et donc P(B\A) = P (B). Il n’est pas toujours facile de 
savoir a priori si deux événements sont indépendants. Pour trancher, on peut éventuel¬ 
lement calculer séparément P (An B) et P (A) P (B), et les comparer. 

Pour résumer les différents cas, il est utile d’un point de vue mnémotechnique d’uti- 



lüT Attention : B peut être indépendant de A et C sans que A et C soient indépen¬ 
dants entre eux. 

b Lois de probabilités 

Souvent, on n’a pas accès directement à une description simple des évènements, 
mais on effectue des mesures dont le résultat dépend de l’évènement qui s’est produit. 
On parle de variable aléatoire. Ces variables peuvent prendre des valeurs discrètes ou 
continues. 

a Lois de probabilités discrètes 

Soit X une variable pouvant prendre un nombre de valeurs discrètes (i.e. on peut les 
mettre en bijection avec N), aq, ■ ■ ■, Xk, ayant comme probabilités, respectivement, 
Pi, ■ ■ ■, Pjç, ■ ■ ■, et telles que Pi + • • • + P& + • • • = 1. 

► Exemple : Soit X la variable aléatoire correspondant au nombre de garçons dans 
les familles de 3 enfants. Les valeurs (discrètes) pouvant être prises par X sont xq = 
0, x\ = 1, X 2 = 2 ou x% = 3. Le nombre total d’issues associées à cette variable aléatoire 
est 8 = 2 3 (deux possibilités pour la première naissance, deux pour la seconde et deux 
pour la troisième; symboliquement GGG, GFF, FGF, FFG, GGF, GFG, FGG, FFF). 
En supposant égale la probabilité de naissance des garçons et des hiles, les valeurs de 
la distribution de probabilité sont : 




2. Quelques rappels sur les probabilités et les distributions 


(3 Lois de probabilité continues Densité de probabilité 

On rencontre également fréquemment le cas où la variable X peut prendre ses valeurs 
sur un ensemble continu I. On parle alors de variable aléatoire continue. 

Par exemple, la vitesse d’une particule selon un axe v x . Pour un évènement donné, 
la variable aléatoire X va prendre une valeur x avec une probabilité qu’on appelle loi de 
probabilité. On parle également de densité de probabilité , de distribution statistique ou 
plus simplement de distribution. C’est une fonction de la variable x. La loi de probabilité 
est obtenue comme la probabilité de l’union de tous les évènements qui donnent la 
valeur x. Dans le cas discret où X peut prendre les valeurs x n . on la note P(x n ) pour 
“la probabilité d’obtenir x' n ”. On a bien sûr P(x n ) = 1. Lorsque la variable aléatoire 
X est continue, elle prend typiquement des valeurs dans un intervalle [a, 6] (les bornes 
pouvant être ±oo). Si on s’intéresse à la probabilité que la variable aléatoire tombe dans 
l’intervalle élémentaire entre x et x + dx, on introduit la densité de probabilité V{x) telle 
que : 

P(x < X < x + dx) = V(x)dx . 


La fonction V(x) est toujours positive. La densité de probabilité est normée: 

r*b 


[ V{x)dx = 1 
J a 


On a également la relation : 


V(x) = 4- P (a < X < x) . 
dx 

En physique, on utilise également le mot distribution pour V(x) et elle aura la dimension 
de l’inverse de celle de x. 

► Exemple : A titre d’exemple, considérons une variable aléatoire continue X qui 
ne prend de valeur que dans l’intervalle [0,4], et soit V(x) = ax + 1/2 sa densité de 
probabilité, où a est une constante. On demande de calculer la probabilité que X prenne 
une valeur comprise entre 1 et 2. 

7 Fonction de répartition 

C’est la fonction associée à la variable aléatoire X et notée Fx, telle que Fx(x) = 
P(X < x). Sauf subtilité mathématique, c’est la primitive de la densité de probabilité 
associée telle que 

Fx(x) = ( V(y)dy . 

On a immédiatement quelques propriétés : 

* Fx est une fonction croissante ; 

’ 0 < Fx( x ) < 1. 

Si X est discrète, la fonction de répartition sera en escaliers. En pratique, il est parfois 
plus facile de calculer Fx(x) et de la dériver pour obtenir la densité de probabilité. 
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Annexe A. Probabilités 


ô Valeur la plus probable, moyenne, variance 

La valeur la plus probable x* maximise la loi de probabilité P(x) ou la densité 
V(x). Il peut y en avoir plusieurs. La m.oyenne d’une variable aléatoire discrète s’écrit 
naturellement : 

(x) = y^x n P(x n ) 

n 

et en version continue 

[ b 

(x) = / xV(x)dx . 

J a 

On peut aussi prendre la valeur moyenne d’une fonction /(x) qui se calcule selon: 

(/) = f f{x)V{x)dx , 

J a 

et vérifier sans peine la linéarité de cette opération : 

(A/ + pg) = A (/) + p (g) . 

Par exemple, la variance o\ est définie par : 

°x = d x - (*)) 2 ) = ( x2 ) - ( x ÿ ■ 

Cette quantité caractérise l’écart à la moyenne : ax, qui a la même dimension que x, est 
appelé écart type. La variance, parfois appelé aussi écart quadratique m.oyen, est le carré 
de l’écart type. En statistique, on parle de dispersion des résultats autour de la valeur 
moyenne. En physique, l’écart type caractérise les fluctuations de la grandeur étudiée. 

Attention : il peut arriver que la variance soit infinie, comme pour la loi de Cauchy 
Vix) = définie sur [-oo,+oo]. 

c Moments d’une loi de probabilité 


On peut définir, de façon analogue, un moment d’ordre n, qu’on notera /x n , de la loi 
de probabilité continue V est le nombre, lorsqu’il existe, définie par la relation : 

Pn = [ x n V(x)dx = (x n ) v , 

J K 


où la dernière expression se lit “moyenne de x n avec la loi P". Ainsi : 

■ le moment d’ordre 0 vaut 1 par la propriété de normalisation et n’apporte pas 
d’information ; 

* le moment d’ordre 1 se confond avec la moyenne de la variable aléatoire x de la 
loi P : 

(x ) p = pi ; 

* le moment d’ordre 2 est associé à la variance , qui peut s’écrire définie par la 
relation : 

ox = H2~ Pi ■ 

Les moments d’une loi de probabilité discrète sont définis par les relations : 


l^n 


E 


T n p 
x m r rr 


m 
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2. Quelques rappels sur les probabilités et les distributions 


d Transformations sur les lois de probabilité 


a Changement de variable 

Supposons que Y soit une variable aléatoire, fonction de la variable aléatoire X 
qu’on notera y = f(x). Connaissant la loi de probabilité Vx{x) pour X, on cherche 
à déterminer celle de Y notée V Y (y). Lorsque y est une fonction croissante de x, on 
montre la relation V Y (y) = Vx(x)/f'(x) dans laquelle il faut remplacer x par / -1 (y). 
On la mémorise souvent sous la forme : 


V Y (y)dy = V x (x)dx . 


Il faut bien noter que la valeur de x dans le membre de droite est x 
N.B. : dans le cas général (pas souvent utile en pratique), on a: 


T’y (y) 


£ 

Xi/y=f{xi) 


1 

\f'(Xi)\ 


Vx{xi) . 


f 1 (y)- 


(3 Loi de probabilité jointe à plusieurs variables 

Lorsqu’on s’intéresse à plusieurs variables aléatoires (on les prendra continues, les 
définitions dans le cas discret étant similaires), on peut introduire la loi de probabi¬ 
lité jointe à ces variables. Pour deux variables aléatoires X et Y, on a la densité de 
probabilité jointe à XY : 

P(x < X < x + dx et y < Y < y + dy) = Vxy(x, y)dxdy , 


avec 


Xb Vb 


Vxv{x,y)dxdy = 1 . 


%a Va 


On peut définir une loi marginale, qui est la densité de probabilité d’avoir X = x, 
indépendamment de la valeur de Y. On intègre donc sur toutes les valeurs possibles que 
peut prendre la variable Y : 


Vb 


Vx(x) = / V X y{x, y)dy . 


Va 


Par exemple, pour calculer la valeur moyenne de x, on a : 

Xb Vb x b 

(x) = (f xVx Y {x,y)dxdy = [ xVx{x)dx . 


%a Va 


Un résultat similaire existe pour les écarts types ux et <jy. On a également pour 
une somme de variables aléatoires : 


(x + y) = (x) + (y) . 

On peut définir l ’indépendance d’une variable par rapport à une autre selon: 

Vxy{x, y) = Vx(x)V Y (y) , 

ce qui entraîne (xy) = (x) (y). On peut généraliser cela au cas de N variables aléatoires. 
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Annexe A. Probabilités 


7 Corrélations 

Deux variables aléatoires X et Y sont corrélées si leur coefficient de corrélation 


Cxy 


(xy) - (x) (y) ^ Q 

a x crv 


Xb Vb 

où, par définition, (xy) = f f xyVxY (x, y)dxdy. On remarquera que deux variables 

*£<z Va 

indépendantes sont décorrélées (le coefficient de corrélation est nul) puisque (xy) = 
(x) (y) ; ceci est intuitif et justifie l’expression du numérateur (on notera également que 
(xy) — (x) (y) = ((x — (x))(y — (y))))- Attention, la réciproque est fausse. Le dénomina¬ 
teur est choisi pour que le coefficient de corrélation soit compris entre —1 et 1 (il peut 
en effet être négatif, on parle d’anti-corrélation). Ce résultat découle de l’inégalité de 
Schwartz | ((x — (x))(y — (y))) | < \J((x — (x)) 2 ) ((y — (y)) 2 )- En physique, on utilise 
souvent la dénomination “les corrélations” pour désigner ce type de quantité et on ne 
met pas toujours le dénominateur. On parle par exemple des “corrélations de densité” 
pour la quantité (n(r)n(r')) — (n(r)) (n(r')) qui sont utiles pour caractériser la structure 
du fluide ou du solide. 


e Variance d’une somme de variables aléatoires indépendantes 

On sait déjà que (x + y) = (x) + (y). On calcule (a x+y ) 2 ■ 

( a x +YŸ = ((x + y) 2 ) - (x + y ) 2 

= (x 2 ) + (y 2 ) + 2 (xy) - [(x) 2 + (y ) 2 + 2 (x) (y)} 

= ox+ °Y + 2 [(xy) - (x) (y)} , 


d’où on tire que, si les variables X et Y sont indépendantes, 

{(y x+y ) 2 = <yx +oy , 

et par récurrence pour une somme E = YliLi de N variables aléatoires Xi indé¬ 
pendantes (autrement dit V(x\, X 2 , ■ ■ ■, xjv) = 'Pi(xi)T > 2(x2) ■ ■ - Vn(xn))- La valeur 
moyenne vaut (S) = (xf) et la variance : 

N 

„ 2 _ \ " ^ 2 

• 

i =1 


f Loi des grands nombres et théorème de la limite centrale 

On considère que les X{ indépendantes ont toutes même valeur moyenne (x) et même 
variance a y et que celles-ci sont finies (on peut imaginer par exemple pour simplifier 
que toutes les X, t suivent la même loi de probabilité), on a d’après le résultat précédent : 

(S) = N (x) et as? = Nax , 

où E = 'ffiLi Xi est la somme des variables. Sa variance vaut, en valeur relative 

_ 1 <?x 

(S) ~ y/N (x) ' 
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3. Exemples importants de distributions statistiques 


Cette remarque est très importante pour les expériences : supposons qu’on mesure une 
grandeur X qui fluctue un peu en raison d’un processus aléatoire (par exemple une 
source de bruit) et que les mesures sont indépendantes. Alors, plus on fait de mesures, 
plus on améliore la barre d ’erreur ou encore l ’incertitude relative sur la détermination 
de la valeur moyenne de X. 

Le théorème de la limite centrale va plus loin en déterminant vers quelle loi de 
probabilité converge la distribution de S. On introduit la variable aléatoire 

y _ X-N(x) 

VNa x ’ 


alors, selon ce théorème, dans la limite N —>• oo, la distribution de la variable Y est une 
gaussienne centrée sur 0 et de variance unité: Vy{y) — -X=eT v I 2 . On peut également 
écrire : 


V*(S) ~ 


(S-N(x )) 2 

2N °X 


2irNcr 


2 

X 


Dès qu’on échantillonne un ensemble de variables aléatoires indépendantes qui sont 
gouvernées par la même loi, la valeur moyenne obtenue sera aura une loi qui tend (assez 
vite en général) vers une gaussienne. C’est la raison pour laquelle les distributions statis¬ 
tiques expérimentalement observées sont très souvent gaussiennes. En mathématiques 
également, de nombreux processus convergent vers la distribution gaussienne. 

Un exemple de loi de probabilité à laquelle on ne peut pas appliquer ce théorème 
est la loi de Cauchy évoquée plus haut. Ces situations sont détaillées dans les exemples 
qui suivent. 


3 Exemples importants de distributions statistiques 


a La loi binomiale 


Si un évènement a une probabilité p de se produire (et donc 1 —p qu’il ne se produise 
pas), après n épreuves, on a la probabilité 

Pn(k)= (fyp k (l-p) n ~ k 

que l’évènement se soit produit k fois. Le facteur ()(j tient compte du fait qu’on ne s’in¬ 
téresse pas à l’ordre dans lequel les évènements se sont produits. C’est la loi binomiale. 
On voit qu’elle est bien normalisée 


Y Pn(k) = {p + 1 - p) n = 1 , 

k =0 

d’après le développement du binôme éponyme. Pour calculer sa valeur moyenne et son 
écart type, on peut utiliser les relations sur les coefficients binomiaux. On note également 
q = 1 — p pour alléger les notations. Par définition, 


(k) = YkP n (k) = Yk 


fc =0 


k=\ 


U X p k q n ~ k 


= Y n 

k= 1 


n — 1 
k — 1 


n— 1 

p k q n ~ k = npY 
k '=0 


n — 1 
k' 


pk'qu-l-k' 


= np . 
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Une autre manière, plus astucieuse, est d’introduire une fonction 

f(x) = (px + q) n = J2 (fyx k p k q n ~ k , 

pour laquelle on a /( 1) = 1. On remarque alors que la moyenne est reliée à la dérivée 
première : 


(k) = ^ kp k q n k = f'(x = 1) = np{p x 1 + q) n 1 = np . 
k =1 ' ' 

On s’intéresse maintenant à l’écart type. On devine qu’il faut calculer la dérivée seconde : 
f"(x = !) = è Q k ( k ~ l )p k q n ~ k = ( k2 ) ~ (k) = n{n - l)p 2 , 


d’où on tire 


a 2 = ( k 2 ) — (k) 2 = n(n — l)p 2 + np — (np) 2 = np( 1 — p) . 

b Distribution normale ou Gaussienne 


Cette loi est la principale loi à connaître car elle intervient comme limite de nom¬ 
breuses autres lois et se retrouve par conséquent dans de nombreux domaines de la 
physique et des statistiques. Elle s’écrit en fonction d’une variable aléatoire continue 
x G [— 00 , 00 ] sous la forme : 


V(x) 



(x - a? ) 

2 b J ' 


On vérifiera que V est bien normalisée. La nouvelle variable aléatoire Z définie à partir 
de X par la définition : 

<7 

est associée à la distribution normale, dite centrée réduite, définie par la relation : 


V(z) 



que l’on a représentée sur la figure A 2. 

À l’aide des intégrales gaussiennes (voir annexe B), on trouve que: 


(x) = a, a 2 = b , 


si bien qu’il est utile de la retenir sous la forme 


V(x) 


1 

y/2i rcr 2 



{x-(x)) 2 \ 
2 ^ / 
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3. Exemples importants de distributions statistiques 



Figure A 2 - Loi de distribution normale. 


Un calcul numérique permet de montrer que P {—1 < z < 1) ~ 68,27% et que 
P (—3 < z < 3) ~ 99,73%. Les valeurs un tant soit peu écartées de l’origine sont donc 
hautement improbables dans le cas de la loi normale. 

c La gaussienne comme une limite de la loi binomiale 

On s’intéresse au comportement de la loi binomiale dans la limite d’un grand nombre 
de réalisations n lorsque p est fixée. On va montrer que la loi est alors très proche 
d’une gaussienne. On suppose que ?j > 1 et que (k) = np > 1 et on cherche le 
maximum du logarithme de la distribution lnP(fe). Ce résultat peut également être 
trouvé en considérant k comme une somme de n variables aléatoires indépendantes X{ 
qui valent 1 si l’évènement est réalisé avec p et 0 sinon. On a k = Y27=i x i- On trouve 
(xi) = p X 1 + (1 — p) x 0 et ax = p(l — p). D’après la loi des grands nombres et le 
théorème de la limite centrale, la loi binomiale doit bien tendre vers une gaussienne de 
moyenne np et de variance np( 1 — p). 

d La loi de Poisson 


C’est une loi paramétrée par un réel À > 0 telle que la probabilité d’avoir l’entier 
k > 0 suive la loi 

On voit qu’elle est bien normalisée : 


OO OO 

E P W=e-^- 


e A e A = 1 . 


Si A < 1, la probabilité est maximum pour k = 0 puis décroît rapidement avec k. Sinon, 
la loi passe par un maximum pour un certain k* (voir ci-après le lien avec la gaussienne). 
Le calcul montre que sa moyenne est égale à sa variance et vaut 


(k) = er 2 = A . 


Cette loi intervient pour décrire les évènements rares et qui se produisent de façon 
indépendante. 
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Annexe A. Probabilités 


e La loi de Poisson comme une limite de la loi binomiale 

On va voir que la loi binomiale, dans la limite où n > 1 mais np fixé (et donc p —>• 0, 
“évènement rare”), tend vers une loi de Poisson. Ecrivons pour k fixé: 

fn\ fr/n nn-fr n(n — 1) • • • (n — k + 1) fc/1 , n _ k 

LJp fc (i-p) n k = - -p k (i-p) n 

(j P n ) k n(n - 1) • • • (n - k + 1) _ k ^ nln(l _ p) 

k\ n k 

kl n n 

~ k\ 

Le comportement est donc très proche d’une loi de Poisson de paramètre À = np = (k) 
(nombre fixé). On remarque d’ailleurs que dans la limite considérée, la variance de la 
loi binomiale a 2 = np( 1 — p) ~ np = A qui est bien le résultat sur la loi de Poisson. Là 
encore, c’est une loi un peu plus simple que la loi binomiale. 

f La gaussienne comme limite de la loi de Poisson 

On peut voir que dans la limite où A > 1 (et on considérera donc k comme une 
variable continue), la loi de Poisson ressemble elle aussi à une gaussienne. Pour cela, on 
procède comme pour la loi binomiale en utilisant la formule de Stirling : 

ln P{k) ~ k ln A — A — k ln k + k . 


On trouve que le maximum est atteint pour k = A et la dérivée seconde donne un facteur 
— 1/À si bien que dans cette limite : 


V(k) ~ - exp 

V2ttÂ 


(k - A) 2 ) 
2A J ’ 


c’est-à-dire une gaussienne avec les mêmes valeur moyenne et variance que la loi de 
Poisson. 
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Annexe B 


Intégrales gaussiennes 


Elles interviennent très souvent en physique, particulièrement en physique statis¬ 
tique comme on va le voir. Les relations dont dérivent les autres sont : 


- a (x- b ? dx = 


et 



On donne également le résultat sur les intégrales du type : 

+ OO 


I, 


.(«> = / 
o 


= / x m e~ ax dx . 


Par dérivation par rapport à a et intégration par partie, on a I m . 1-2 (a) = Aa' A;.( a )- 
On distingue les résultats sur les cas pair et impair. On n’aura besoin que des premiers 
termes : 


1 TT 


/o(a) “2Vâ ! /2(a) “ 


1 TT 


4 a\ a 


, h{a) = 


8 a 2 V a ’ 


h(a) = Y a ■ J 3(a) = ^ ’ A(®) = ^ 3 a . 
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Annexe B. Intégrales gaussiennes 


108 




Annexe C 


Volume d'une hyperboule et fonction 
gamma 


On cherche à calculer Vd(R), le volume de l’hyperboule de rayon R dans un espace de 
dimension D. Par homogénéité, on doit avoir Vd(R) = CdR D où Cd est une constante 
numérique que l’on cherche à déterminer et qui dépend de la dimension D. On utilise 
l’astuce suivante 1 : on calcule de deux manières l’intégrale 

/ oo /*oo 

dx i • • • / dx D e-(*i 2+ - +a -d) 

-OO J —OO 


La première manière, c’est de voir que e ( x i+---+ x d ) = Y\9 = 1 e x i et donc que 


Id = 


( r-\-oo \ D 

J oc dxe ~ X J = {V * )D 


La deuxième manière est de remarquer que, en passant en coordonnée sphérique à 
D dimensions avec le rayon r 2 = x' 2 + ... + x£, on peut écrire 

/»oo 

Id= e~ r ~S D (r)dr , 

J o 

avec Sf)(r) la surface de l’hyperboule de rayon r. Comme Sr>(r) = —^ - = D Cd r D ~ l 


on peut écrire 


dr 


I d = DC d r D ~ 1 e~ 7 ' 2 dr . 


On fait alors le changement de variable u = r 2 (du = 2 rdr) pour avoir 

D /*oo 

I D = —C D / u^~ 1 e~ u du . 

1 J o 

L’intégrale restante est la définition de la fonction T(x = D/2). En effet, on définit 
la fonction gamma par 

roo 

r(x)= / u x ^ 1 e~ u du . 

J o 

En intégrant par partie T(x), on montre la relation fondamentale (généralisation de 
la factorielle) 

r(x + 1) = xT(x) . 


1. Qui en fait vient de l’isotropie de l’intégrale choisie et de la factorisation de l’exponentielle. 
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Annexe C. Volume d’une hyperboule et fonction gamma 


Pour un entier n, on a en particulier T(n + 1) = ni ce qui donne ici : 


Id 


C D y r 



c D r (| + 1 ) = ^ 


On obtient finalement : 


Vd(R) 


ir D / 2 

r(f+ 1) 


r d 


Vérifions qu’on retrouve bien la longueur d’une ligne, la surface du disque et le volume 
de la sphère pour D = 1,2,3. Pour cela, on a besoin de savoir que r(l/2) = sJtt (en 
refaisant le changement de variable inverse à toute à l’heure, r = y/u, on se ramène à 
une intégrale gaussienne bien connue). 

* Pour D = 1, on a T(3/2) = \/tt/2 qui donne V\(R) = 2 R. C’est bien la longueur 
d’une ligne de rayon R. 

■ Pour D = 2, on utilise r(l) = 1! = 1 et on retrouve bien la surface du disque 
V 2 (R) = n R 2 . 

* Pour la sphère, avec r(5/2) = |T(3/2) = |\/vr, on a bien V 3 (R) = ^-kR?. 

’ On peut calculer également le volume de la boule de dimension 4, on trouve 

V 4 (R) = ^-R 4 . 

Enfin, à l’instar de la factorielle, on a pour les x grands la formule de Stirling suivante : 


r(x + 
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Annexe D Multiplicateurs de Lagrange 


Considérons une fonction / de n variables indépendantes: f(x\, ■ ■ ■, x n ). Sa diffé¬ 
rentielle au point x = (aq, • • •, x n ) s’écrit 


df 


df , 

tt- dx 1 + 

ÔX\ 


+ 


df_ 

dx n 


dx r , 


Une condition nécessaire pour que la fonction / possède un extremum (maximum ou 
minimum) au point x est que df = 0, soit : 


df_ 

dx\ 


df_ 

dXr, 


= 0 


(A3) 


Supposons maintenant qu’il existe une contrainte entre les variables aq, ■ ■ ■, x n ; au¬ 
trement dit, il existe une certaine relation fonctionnelle 


g{x 1 , • ■ ■ ,x n ) = 0 

entre ces variables. Cette relation montre qu’une des variables, disons x n dépend des 
n—1 autres. Il n’y a donc plus n variables indépendantes, mais seulement n—1. Ainsi, les 
conditions d’extrémalisation données par (A 3) ne sont pas valables en présence d’une 
contrainte. 

Une première solution consiste, lorsque cela est possible, à exprimer une des variables 
en fonctions des n — 1 autres, à substituer son expression dans /, puis à écrire les 
conditions d’extrémalisation (A3) sur les n — 1 variables restantes. 

► Exemple : Soit à déterminer le rectangle d’aire maximale pour un périmètre donné. 
Avec des notations évidentes, on a L(x,y) = 2(x + y) et S(x,y ) = xy. En éliminant y 
entre ces 2 équations, on obtient S(x ) = x(L/2 — x), x étant maintenant une variable 
non contrainte. La condition S x = 0 conduit à la relation L/2 — 2x = 0, et donc à la 
solution x = y = L/4, c’est-à-dire, à un carré. 

Un autre procédé plus général a été proposé par Lagrange. L’idée consiste à absorber 
la contrainte dans la définition d’une nouvelle fonction. Soit en effet la fonction / ± Xg, 
où A est un paramètre appelé multiplicateur de Lagrange. Il s’agit d’une fonction des n 
variables indépendantes (aq, • • • ,x n ) puisqu’aucune contrainte extérieure ne s’applique 
désormais sur les variables. On peut donc écrire directement : 

d f dq 

d(f±Xg) = 0 => — ±A— = 0, i = l,---,n. 

OXi OXi 

Ces n équations et la contrainte g = 0 déterminent donc les n+1 inconnues x\, ■ ■ ■, x n , A. 
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Annexe D. Multiplicateurs de Lagrange 


► Exemple : Reprenons l’exemple précédent. On a f(x, y) = xy et g(x, y) = L— 2(x + 
y ) = 0. L’extremum de / est déterminé par les 2 relations f x + Lg x = 0 et f y + \g y = 0, 
soit y — 2A = x — 2A = 0, qui conduit bien à la solution x = y. 


Théorème : Soient f et g deux fonctions C 1 définies sur M n . On suppose en outre 
qu’il existe i tel que dg/dxi 0. 

Un extrem.um (maximum ou minimum) de la fonction f soumise à la contrainte 
g(x i, • • •, x n ) = 0 est déterminé (s’il existe) par les n + 1 équations : 


ÊL 

dx i 


= o 

ÔX-i 


i = 1, ■ ■ ■, n 


g{x i,-’ • ,x n ) = 0 , 


où A est un paramètre réel appelé multiplicateur de Lagrange. 


En introduisant autant de multiplicateurs qu’il y a de contraintes, on généralise ce 
procédé aux cas des fonctions soumises à plusieurs contraintes. 
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Annexe E Transformation de Legendre 


Il est bien connu qu’une courbe dans le plan y = f(x) est parfaitement définie dès lors 
que sa fonction dérivée f est connue en tout point x. On rencontre un certain nombre de 
situations en physique et en mathématique où la donnée naturelle du problème étudié 
est la dérivée f'(x ) plutôt que la variable x elle-même. La transformation de Legendre 
est un moyen systématique qui permet de définir, sans perte d’information, une nouvelle 
fonction dont la variable indépendante sera la pente de la fonction / au point x, c’est- 
à-dire f'(x). 

Définition : Soit f : M i —> R, une fonction bornée à variation bornée. 

La transformée de Legendre de f est la fonction g définie par l ’une des deux 
relations : 

9 = xf(x) - f(x) ou g (f'(x)) = f(x) - xf'(x) . 

La première convention de signe est celle utilisée en Mathématique, la seconde 
est celle utilisée en Thermodynamique. 

Montrons que cette définition est équivalente à la définition 

gip) = max(pæ- f{x)) , 

qui a un contenu géométrique clair : g(p) est la distance maximale entre la droite de 
pente p qui passe par l’origine et la fonction f. La solution est obtenue pour x p tel que 
p = f'(xp), c’est-à-dire que la distance entre la droite et la fonction est maximale au 
point x p où la tangente à f est égale à la pente de la droite (voir Fig. A3). On a donc 
aussi g(p) = px p — f(x p ), soit encore g(f(x p )) = f'(x p )x p — f(x p ), ce qui correspond 
à la définition initiale. Les variables x et f'(x ) = p sont dites variables conjuguées par 
rapport au couple de fonctions / et g. 



FIGURE A3 - Interprétation géométrique de la transformée de Legendre. 
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Annexe E. Transformation de Legendre 


► Exemple: Montrons que la transformée de Legendre de la fonction f(x) = mi 2 /2 
est la fonction gfp) = p 2 /2m. En effet, comme p = f'(x) = mx. on a par définition 
g(f'(x )) = x(mx) — mx 2 /2 = mx 2 / 2, donc g(p) = p 2 /2m. Cet exemple simple montre 
que la transformée de Legendre peut être utilisée pour obtenir la mécanique hamilto¬ 
nienne à partir du formalisme lagrangien. 

Voici quelques exemples supplémentaires de transformées de Legendre que l’on éta¬ 
blira à titre d’exercices : _ 


/(*) 

g(p ) 

x n 

n 

e x 

p ln p — p 

lnx 

Inep 


On vérifiera sur ces exemples que les dérivées de / et g sont des fonctions inverses l’une 
de l’autre, ce qui, du reste, est une autre définition de la transformée de Legendre d’une 
fonction. 

Le passage de f(x) à g(f'(x)) peut d’ailleurs être retrouvé directement à partir de 
la définition. Soit en effet une fonction g définie à partir de x. f et f : 

d(x, f(x), f'(x)) = f{x) - f\x) x . 

En utilisant le fait que df{x) = f{x)dx (/ est différentiable en x), on voit que la 
différentielle totale de la fonction g au point x est telle que 


dg = df(x) — f'(x ) dx — x df'(x) = — x df'(x) . 


Il en résulte que d x g = df^g = 0 de sorte que g est bien une fonction différentiable de 
la seule variable indépendante f'(x). 

Cette transformation se généralise aux cas des fonctions de plusieurs variables où il est pos¬ 
sible de faire une transformation de Legendre associée à une partie (ou à la totalité) des variables. 
Considérons ainsi une fonction différentiable / des n variables x\, ■ ■ ■ , x n et supposons que l’on 
cherche à définir une nouvelle fonction g qui dépende des variables x\, ■ ■■*, x r , u r +i, où les 

Ui = df /dxi (i = r + 1, • • •, n) sont les variables conjuguées des variables initiales • • •, x n . 
En généralisant la procédure à une dimension, il suffit de définir 


g = f - {u r+ ix r+ i H-h u n x n ) 


pour constater que 


dg = u\ dx i + • • • + u r dx r — &>+i du r +1 — ■ ■ ■ x n du n , 


où on a utilisé le fait que df = i u i ^ x i ■ 9 est donc une fonction différentiable des n 
variables Xi, ■ ■ ■, x r , u r + 1 , 

On notera en outre une importante propriété de la transformation de Legendre : 
même de dg montre qu’on a 


df_ 

dxj 


dg 

dui 


pour i = r + 1, • • •, n . 


nouvelles 

l’écriture 
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et Xi = 





► Remarque : En mécanique analytique ou en mathématiques, on prend plutôt comme 
définition G (y) = yx — F{x), ce qui change quelques signes et la convexité de F dans 
la discussion ci-dessus (cf la définition du hamiltonien). 

► Exemple : En thermodynamique, on peut écrire l’énergie interne du système E, ou 
U, en fonction des variables S, V, N sous la forme E(S, V, N). La variable conjuguée de 
l’entropie S par rapport à E définit la température T : 



Si on veut passer du jeu de variables (S, V, N) au jeu de variables (T, V, N), on effectue 
une transformation de Legendre sur la variable S en introduisant une nouvelle fonction 
d’état F(T,V,N), appelée énergie libre, et définie donc par 


F = E — T S . 


En thermodynamique, on raisonne plutôt à l’aide des différentielles. On a l’identité 
thermodynamique : 

dE = TdS - pdV + pdN , 

si bien que 

dF = TdS - pdV + pdN - TdS - SdT = -SdT - pdV + pdN , 


où l’on voit bien que F dépend naturellement de la variable T (puisqu’on a dT et non 
plus dS dans la différentielle), tandis qu’on retrouve l’entropie en utilisant 
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Annexe E. Transformation de Legendre 
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Annexe F Description lagrangienne et hamiltonienne 
d'un système classique 


On résume ici le cadre formel de la description d’un système classique, dans le 
formalisme lagrangien d’abord, puis le formalisme hamiltonien ensuite. Ce dernier se 
révèle mieux adapté à la description dans l’espace des phases. 

1 Equations de Lagrange 


On décrit géométriquement la dynamique d’un système par le choix de N coordon¬ 
nées généralisées souvent notées qi(t), i = 1,..., N, paramétrée par le temps t. Il s’agit 
souvent des positions en coordonnées cartésiennes mais on connaît également les angles 
qui décrivent l’orientation d’un solide ou d’autres formes de coordonnées. 

La dynamique étudie l’évolution dans le temps du système si bien que l’on s’intéresse 
à la manière dont ces coordonnées évoluent en introduisant les vitesses associées q t ( t ) = 

d’une part, et une fonction £(qi, qi, ..., çn, Ç[n, t) appelée lagrangien et dont on 
peut tirer les équations du mouvement. On a donc 2N variables indépendantes en tout 
et par suite, on condensera la notation (qi(t), qi(t),.. . , gjv(t), gjv(A) en pour 

alléger les équations. 

Il s’agit d’une reformulation des lois de Newton, équivalente, mais qui fournit une 
description et des méthodes plus génériques, et parfois plus simples. Toute l’information 
sur la dynamique du système est alors contenue dans la définition du lagrangien. Dans le 
cas relativement général où les forces dérivent d’un potentiel U(qi), en notant K(qi,qi ) 
l’énergie cinétique, le lagrangien s’écrit C(qi,qi) = K(qi,qi ) — U(qi). Les N équations 
du mouvement gouvernant les (qt, q t ) s’obtiennent à partir des équations de Lagrange : 


A ( 9C \ 9C n 

dt \dçu J dqi 


avec la dérivée totale 


d_ 

dt 


d_ dçu _Ô_ dqj d 

dt, dt dqi dt dqi 

1=1 


(A4) 


Elles sont du second ordre par rapport à t. Il faut bien entendu rajouter les conditions 
initiales pour obtenir la trajectoire effectuée par le système. 


2 Formalisme d’Hamilton 


On peut récrire les équations précédentes avec un nouveau jeu de variables, les ( qi,Pi ) 
où les moments généralisés sont définis par 

dC 
dqi ' 
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Annexe F. Description lagrangienne et hamiltonienne d’un système classique 


Le passage à ces nouvelles variables introduit V hamiltonien du système par trans¬ 
formée de Legendre du lagrangien : 


(A 5) 


N 

U(qi,Vi) = 5 >« - £(<?*, çn) 

i =1 


En utilisant (A4) et (A5), les équations du mouvement prennent alors une forme simple 
sous la donnée de 2N équations du premier ordre par rapport au temps : 



(A 6) 


3 Évolution dans l’espace des phases (hors programme) 

L’espace des phases est l’espace à 2N dimensions dont les axes représentent les va¬ 
riables • • •, çn^Pn)- La trajectoire du système y est alors une courbe paramétrée 

par le temps t. 

Théorème de Liouville : le volum.e élémentaire de l’espace des phases 
N 

dr = dqidpi autour du point ( qi(t),Pi(t )) est conservé au cours de 
i =1 

l’évolution temporelle, le long de la trajectoire. 

DUP Le montrer en exercice. Indications : il faut considérer de passage d ’un point de 
la trajectoire de t à t + dt et voir comment évolue le volume élémentaire lors de ce 
changement de base. 


Évolution d’une densité de probabilité sur l’espace des phases : On peut 
considérer qu’un état macroscopique est déterminé à un instant t par une densité de 
probabilité dans l’espace des phases notée u(qi,pi ,..., qig,pN',t). En utilisant le théo¬ 
rème de Liouville et la conservation de la probabilité lors d’une évolution entre t et 
t + dt, on obtient l’équation de Liouville : 



(A 7) 


en utilisant la notation de la dérivée totale de l’équation (A4). La densité de probabilité 
s’écoule donc de manière analogue à un fluide incompressible dans l’espace des phases. 
03? Montrer l’équation (A 7) en exercice. 

En mécanique statistique classique, on travaille avec le formalisme hamiltonien parce 
que l’espace des phases se prête bien au passage d’une moyenne sur le temps d’une 
observable à une moyenne d’ensemble, en partie en raison du théorème de Liouville. 
Une autre raison est que, pour un système présentant la symétrie d’invariance dans 
le temps, l’hamiltonien correspond à la quantité conservée qui est l’énergie totale. Les 
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j. Exemples d’hamiltonien élémentaire 


trajectoires possibles à énergie constante sont donc confinées dans un hyperplan. On se 
convaincra également que deux trajectoires distinctes ne peuvent se couper dans l’espace 
des phases. On notera de plus que, parmi les densités de probabilité stationnaires (telles 
que = 0) correspondant à un hamiltonien indépendant du temps 1 , les densités qui 
s’écrivent comme fonction des constantes du mouvement (une constante du mouvement 
G est telle ^ = 0) satisfont naturellement à (A 7). En particulier, si on prend co(qi,pi ) = 
f('H{qi,Pi)), avec / une fonction quelconque, cette distribution est bien stationnaire 
et on verra que les densités de probabilité qu’on utilise en physique statistique pour 
décrire les états d’équilibre, qui sont ici stationnaires, s’expriment effectivement à l’aide 
du hamiltonien. 


4 Exemples d’hamiltonien élémentaire 


• Cas général. Une particule de masse m dans un potentiel extérieur V(q): 

U(q,p) = ^ + V(q) . 

2m 


On retrouve alors les équations du mouvement de Newton en appliquant (A 6) : 


dPq p _ d 3 p 
dt m dt 


-^V(q) . 


• Oscillateur harmonique. Dans ce cas, V(q) 


1 


mu 2 ||g || 2 


avec U3 la fréquence 


de l’oscillateur. 

• Particule dans un potentiel central V(\\q ||). En passant des coordonnées car¬ 
tésiennes (x, x, y , ÿ, z, z ) aux coordonnées sphériques (r, r, 6 , 9 , ip, <p), on trouve : 


'H{r,p r ,0,po,ip,p v ) 


1 

2m 


Pr 




r sin 9 2 


+ V(r) . 


C’est un exemple où le terme cinétique est fonction non seulement des moments 
généralisés mais également des coordonnées. 

• Deux particules en interaction via un potentiel U. 

U(qi,pi,q 2 ,P 2 ) = MIL + U|k!L -I- U(\\qi - q 2 \\) ■ 

ISP Retrouver le résultat sur le potentiel central en appliquant la définition des mo¬ 
ments généralisés. 

11!? Tracer la trajectoire de la particule dans l’espace des phases dans les cas suivant : 

* oscillateur harmonique à une dimension. 

' particule dans une boîte à une dimension. 

' particule qui rebondit élastiquement sur le sol, soumise à la gravité. 


1. Si H ne dépend pas explicitement du temps, on parle de système conservatif, c’est en particulier 
le cas pour un système isolé. 
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